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STRESSES LAYERED ELASTIC SOLIDS 


SYNOPSIS 


this paper, the general expressions for stress are derived for multi- 
layer systems which satisfy the conditions plane strain. illustrate the 
procedures involved, expressions for stress are derived for two-layer sys- 
tem, subjected surface line strip load. Each layer considered 
elastic solid. The upper layers are considered finite depth, but 
infinite horizontal extent, and the lowest layer infinite depth and horizontal 
extent. The general equations developed herein may used calculate the 
stresses supporting layered media for embankments, continuous footings, 
pavement systems, compressor-station foundation mats, canals, etc., the 
supporting media canbe approximately plane strain 
system. Expressions for stress are developed for conditions the interface 
perfect continuity and zero friction. Such system permits analytical 
solution many problems engineering importance. The procedures em- 
ployed are means limited two-layer systems. The three and four 
layer cases may solved with little additional difficulty with the aid ofa digital 
computer. 


The numerical solutions for the two-layer cases were carried out for 
value Poisson’s ratio equal one-quarter, value which more closely ap- 
proximates that most engineering materials than the limiting values zero 
and one-half used primarily for convenience many the earlier investi- 
gators determining the expressions for stresses and displacements 

radially symmetric problems. Any arbitrary values 1/2) Poisson’s 
ratio could have peen used view the generality the methods which have 
been developed. 

For given value Poisson’s ratio, and given condition the interface, 
the important parameters the two-layer system are most conveniently ex- 
pressed terms the non-dimensional ratios a/h and Numerical 


Note: Discussion open until January extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Engineering Mechanics Division, Proceedings the 
American Society Civil Engineers, Vol. 86, No. August, 1960. 
Mgr. Strength Analysis Sec., Southwest Research Inst., San Antonio, Texas. 
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solutions were obtained for and for 0.1 for the 
two interface conditions. However, the interest brevity, numerical values 
for stress are tabulated this paper for 0.02 only. Many the pave- 
ment and foundation systems engineering importance fall within this range. 
Stresses have been evaluated along the upper surface, interface, and other 
pertinent locations. 


Notation.—The letter symbols adopted for use this paper are defined 
where they first appear, the illustrations the text, and are arranged 
alphabetically, for convenience reference, Appendix 


METHOD ANALYSIS FOR TWO-LAYER PLANE 
STRAIN PROBLEMS 


Introduction.—The statical problem concerning infinite elastic solid 
bounded plane that subjected given distribution surface forces 
(or deformation) has attracted the attention anumber eminent elasticians. 
The first solution for the problem involving purely normal load was worked 
out Lame and Clapeyron(1)2 and was published 1831. Then, 
1884, solved the problem loaded slab liquid subgrade. 

The next important contribution the subject was made 1885 
(3) who solved number the important cases, including the one 
normal concentrated load acting the surface semi-infinite solid. 
1902, Lamb’s (4) solution for the case symmetrically distributed 
loading the surface semi-infinite elastic solid was published. 

plane strain solid revolution, making use stress function, and 
listed general solutions the basic differential equations this function. His 
development the basis most present-day work for layered systems. 

Pickett, (8) and developed formulas for stress and displacement 
for the one-layer system, resting rigid base, for line axial loadings. 

problem thin slab infinite size supported semi-infinite elastic 
solid. Their solutions were restricted symmetrically loaded areas. The 
solution for unsymmetrical loading was not worked out until 1947 
Volterra. (12) 

Burmister (13) published paper 1943 establishing the equations 
stress for the two-layer system subjected radially symmetric loading, 
treating each layer elastic solid having value Poisson’s ratio 
one-half, for conditions perfect continuity and zero friction the interface, 
the upper layer infinite horizontal extent, the lower layer being infinite 
depth. 1945, extended his theory the three-layer The 


Numerals parenthesis—thus: (1)—refer corresponding items the Bibliogra- 
phy. 
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computation stresses from his formulas, for loadings practical impor- 
tance, means simple. The determination the actual numerical val- 
ues for stress displacements points off the axis loading involves con- 
siderable additional computations. 

1948, Fox(15) developed explicit solutions for the two-layer case 
the form integrals for the computation stresses the axis loading, 
utilizing Burmister’s general equations for the two-layer system. The inte- 
grals were evaluated mechanical quadrature. Using 
relaxation method, solved the differential equations such manner 
also obtain the components stress off the axis loading, again using ap- 
proximate numerical methods. 

1951, Acum and extended this work the three- 
layer elastic solid system, again starting with Burmister’s equations, ob- 
tain formulas for the horizontal and vertical stresses along the axis load- 
ing. This time they found the application relaxation methods completely 
impractical terms the labor required for even reasonable accuracy. 

The solutions developed Pickett, Biot, and Marguerre, for the plane 
strain problem elastic solid subjected this type loading, previ- 
ously mentioned, were restricted the one-layer system supported 
finitely rigid base, and such are mere special cases the general two- 
layer solutions developed herein. should also noted that the equations for 
stress developed Burmister and later modified Fox are not applicable 
the physical systems herein investigated, even though their methods can, 
use the Newmark(18) charts, extended include non-circular 
loadings. 

general solution the two-layer problem derived here- 
embodied the usual assumptions the theory elasticity and further re- 
quired that certain conditions the boundaries and interface satisfied, 
follows: 


The material composing each layer was assumed homogeneous, 
isotropic, and elastic, obeying Hooke’s law. 

The modulus elasticity, and Poisson’s ratio, were assumed 
different for each layer, but for convenience performing the numerical 
computations, Poisson’s ratio was assumed the same both layers. 

The boundary conditions required that the top Layer No. free 
shearing stresses, and free normal stresses everywhere except where the 
strip line load was applied, and that the stresses and displacements in- 
finite depth zero. 

The two continuity conditions investigated the interface required that: 
(a) perfect continuity exist, that is, the two layers act together with slip- 
page the interface, and (b) zero friction exist, that is, the layersare con- 
tinuous contact but perfectly free move horizontally relative each other. 
reality, the actual conditions the interface are, course, between these 
two extremes. For brevity, the conditions (a) and (b) are also referred 
the and “smooth” cases, respectively. 

Theory. 

Development the Basic Equations Stress for the Two-Layer System 
with Perfect Continuity the Interface.—The two-layer system investigated 
shown schematically Figs. Positive directions for the and 
coordinates and the respective displacements and are indicated. The 
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surface loading compressive. The solution carried out for sinusoidal 
stress the y-direction and the parameter loading. will shown 
later that the stresses resulting from strip load intensity and width 
line load equal total force, may obtained from the solution for the 
sinusoidal loading application the Fourier integral. 

Equilibrium Equations.—The basic problem solved represents 
problem plane strain. The basic equations for equilibrium, terms the 
displacements are (19) 


which the displacement function for the horizontal displacement and 
the displacement function for the vertical displacement (Hereinafter, 
for convenience, the subscript dropped.) can that the follow- 
ing displacement functions satisfy the equilibrium and boundary equations. 


C3, C4, are zero Eqs. 3(c) and 3(d), order satisfy the condition 

Upon rewriting the equations equilibrium (1) terms the partial 
derivatives the displacements, obtain for layer 


eu, au, a2y, 
and, for layer 
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Performing the required partial differentiations the displacement func- 
tions, inserting them Eqs. 5(a) and 5(b), and equating the the 
each separately equal zero, obtain 


Continuity Equations.—(1) RoughInterface Perfect Continuity Displace- 
ment the Interface. 


For the case perfect continuity the interface, have 


From and Hooke’s law, €y, obtain 
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ELASTIC SOLIDS 


(2) Smooth Interface (Zero friction) 
For the case zero friction the interface, have 


From Eqs. and Eqs. 17, then obtain 


Boundary Equations.—On the top surfaces, have COS ax, 
Using these boundary conditions (and the expressions for and each 
layer), obtain the remaining equations 


A; - Ag - Ag - Aq + By + Bg = 0 


which are applicable both the rough and smooth interface cases. 

Eqs. through 11, through 16, and constitute independent equa- 
tions for solution the rough interface problem terms the unknown 
displacement function constants Ag, By, Bg, Bg, Cy, 
terms the loading parameter Similarly, Eqs. through and 
through constitute independent equations for the smooth case. For con- 
venience, these equations are summarized Tables and for which the 
parameters are expressed non-dimensional terms. 

The solutions the two layer systems derived herein, require solving the 
simultaneous equations for the unknown constants the displacement 
ing the unknown constants, terms these parameters, formal algebraic 
solution impractical. Instead, numerical values, covering broad range 
may assigned these parameters. Numerical values the constants 
were computed for values ranging from 0.05 10, 0.02 and 0.1, 
a/h and 1/4 for both the rough and smooth cases. These 
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TABLE 1.—SUMMARY THE BASIC EQUATIONS FOR THE TWO-LAYER 
SYSTEM ROUGH INTERFACE 


Cy + Ko Uo Cy + ky Uo Dy 2(ky Up + 


-ky My Ag + 2 Aq - ky My Bg - ky Hy Bg = 


My Ay + ky My Ag + ky My By My (5) 
-ky Hy Ag - ky My A4 - ky Hy B3 - 2(ky My + 1)B4 


" 


A, - Ap - Ag - Aq + B, + By = 0 (11) 


constants were computed means digital computer eleven digits 
after the decimal point. Nevertheless, the high values (greater than 
seven), lack precision was noted the evaluation Dy, due 
increased round-off” error. 

increase the precision would have required computing the residuals for 
each solution and adjusting the roots accordingly. Fortunately, however, the 
methods developed herein for the determination stresses are such that the 
contributions the higher a’s the total stress become progressive- 
smaller with the increased Thus, given error the constants, due 
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TABLE THE BASIC EQUATIONS FOR THE TWO-LAYER 
SYSTEM SMOOTH INTERFACE 


Ug Cy + Ky Ug Co + Ky Ug Dy - + = 0 


-k, My Ay + ky My Ag + ky Hy - 2(ky 1)By 
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the inherent round-off” the higher values has little effect, percent- 
age-wise, the accuracy the final stresses, even though the percentage 
error for the constant itself might considered significant. 
Determination Stresses for Strip Line Loadings Using the Fourier 
Integral 
Each solution the simultaneous equations obtained from the digital 
computer was for sinusoidal loading specific frequency and specific 
value the amplitude which function applied the top surface 
the upper layer. Having obtained the numerical values for the constants the 
displacement functions from these solutions, the stress any point due the 
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sinusoidal loading may obtained from the following general expressions for 


stress: 


(24a) 


(24b) 


(24e) 


obtain the stresses for strip line loading from these solutions, was 


then necessary resort procedure which would permit integration over 
This was accomplished expressing the strip loading terms Fourier 


integral 


2 


where the primes indicate the variable integration and p(x) the loading 
function represented. [For detailed discussion the use the Fourier 


integral the reader referred the work Churchill (22)]. The numerical 
values the integrals were obtained plotting the integrand and integrating 


the area under the curve, using Simpson’s rule. 


Figs. and show thecharacteristic shapes typical curves the incre- 


mental stress versus Since the stresses dy, 


for each layer are ob- 


tained integrating these curves from zero infinity, evident that the 


expression for total stress. 


incremental stresses physically represent the changes stress per unit in- 
crease or, mathematically stated, represent the integrands the final 


t 


q 
& 
| 
a 
(25) 
: 
= 
ae 
al 


zt 
‘ 


ELASTIC SOLIDS 


ysnoy 


zt 
‘ 


fe | 
|_| 
|_| 
x 
N 


Summaries for the stresses for each layer for six problems 
are found Tables through negative sign these tables indicates 
compression. 

Stresses Elastic Solid Resting Rough Infinitely Rigid 
Numerous investigators, including Marguerre,(6) Biot,(7) and Pickett(8) have 
considered the problem stresses anelastic solid resting rigid layer. 

The first correct solution this problem, for plane strain, was 
given Marguerre.(6) Pickett(8) also obtained solutions this problem for 
the stresses along the top surface and the interface, for various values 
Poisson’s ratio. 

order provide independent check the compactness the basic 
equations and methods solution for the general two-layer system developed 
herein, special case was considered, with equal infinity. The solution 
this problem should, course, the same Marguerre’s solution the 
problem the semi-infinite elastic solid resting rigid layer. 

The comments herein are restricted the case which perfect continuity 
exists the interface, although the equations derived earlier for the two- 
layer system, with zero friction the interface, are completely applicable 
the solution this problem for the case zero friction. 

For convenience, problem was considered for which =0, 
The vertical displacements the upper layer, the axis 
loading, may expressed 


Using Hooke’s law and Eqs. 25, the vertical stress the axis loading then 
becomes 


noted that represents the increase vertical stress due 
unit increase spectral frequency. Thus, the total vertical stress over the 


The B’s course, are functions and were obtained solving the 
general equations for the two-layer system for the arbitrarily assigned nu- 
merical values 

Fig. plot the integrand the vertical stress Integrating the 


area under this curve, using Simpson’s Rule, value for -0.0697 was 
obtained. 


Marguerre’s (6) solution for vertical stress the interface the axis 
loading 
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Author's curve 


———-— Marguerre’s curve 


Value 


0 0.08 0.16 0.24 0.32 0.40 0.48 0.56 0.64 0.72 0.80 088 0.96 1.04 
Value of @ 


5.—COMPARISON AUTHOR’S INTEGRAND FOR VERTICAL STRESS, RIGID 
BASE PROBLEM, WITH MARGUERRE’S INTEGRAND 


—-—-— Slab on elastic solid 
| ——-—— Slab on liquid subgrade 


+ 


Value 


Value of x 


SOLID, LAAB ELASTIC SOLID, AND SLAB LIQUID SUBGRADE 
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TABLE 3.—SUMMARY STRESSES 


a 

0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
3.2814 0.0000 
2.9769 -0.2137 
2.4390 
1.2242 -0.0994 
-0.2177 0.0000 
-0.0070 -0.0284 
-0.1942 
+0.0040 -0.0514 0.0000 
+0.0040 
+0.0038 
+0.0035 


0.0000 
-4,3034 0.0000 
-3.6033 -1.0000 0.0000 
0.0000 0.0000 
2.8802 -0.2048 0.0000 
2.7259 -0.2013 -0.0175 
2.2978 
-0.0078 0.0000 
-0.0104 -0.2013 -0.0175 
-0.0161 -0.0311 
0.0040 0.0000 
0.0040 
0.0038 
0.0035 -0.0507 


-7.4399 0.0000 
0.0000 
0.0000 
-5.2910 0.0000 0.0000 
4.5790 0.0000 
4.4883 
2.8789 -0.3259 -0.0823 
-0.0308 -0.3800 0.0000 
-0.3715 
-0.3259 -0.0823 
-0.1025 0.0000 
-0.1025 -0.0022 
-0.1022 
-0.1011 -0.0088 
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TABLE 3.—SUMMARY STRESSES (Continued) 


a 

0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
4.3291 0.0000 
4.1033 -0.1999 0.0000 
3.3842 0.0000 
0.0000 
-0.2062 0.0000 
0.0000 
0.0000 
-0.1615 0.0000 
0.0000 
-0.0015 
-0.0029 


Problem Smooth Interface, 0.02, a/h 


0.0000 
0.0000 
-3.7608 0.0000 
0.0000 0.0000 
4.0068 0.0000 
3.8422 
3.4739 0.0000 
2.2690 0.0000 
0.0000 
-0.1929 0.0000 
-0.1849 0.1849 0.0000 
-0.1620 


0.0000 
-0.0059 


Problem Smooth Interface, 0.02, a/h 


0.0000 
0.0000 
-0.3522 0.0000 


-1.0000 0.0000 
0.0000 
0.0000 
0.0000 0.0000 
6.8189 0.0000 
6.6878 0.0000 
6.3001 -0.3522 0.0000 
-0.3239 0.0000 
e 


-0.3239 -0.3239 0.0000 


0.0000 
-0.0030 
-0.1259 


4 
4 


The integrand the above expression also plotted Fig. for pur- 
poses comparison. Performing the integration before, value for 
-0.0702 was obtained. should noted that this value agrees with the 
writer’s values within 0.56%. which excellent considering the precision 
with which the numerical integration was performed. 


DISCUSSION RESULTS AND CONCLUSIONS 


Twelve physically distinct problems were solved, covering fairly broad 
range the parameters, order study systematically the effect stress 
prescribed changes values these parameters. and illustrate the 
methods and techniques. However, the interest brevity, numerical val- 
ues stress are tabulated this paper for only (Table 3). 
Nevertheless, the observations which follow are valid for values 
ranging from 0.02 1.0. 

For convenience performing the numerical calculations was assumed 
since the two-layer system essentially two-parameter problem, expressi- 
ble terms the non-dimensional ratios a/h and evident that 
the stresses obtained from the numerical examples and summarized Table 
are directly applicable the solution problems for which and are 
other than unity, provided that a/h and not change. For problems 
involving distributed loading. and not appear explicitly the final 
expressions for stress. but appear non-dimensional ratios x/h, y/h, and 

For problems involving concentrated line loads (a/h which and 
are other than the stresses found the above-mentioned tables still 
represent the stresses for provided that they are multiplied the 
factor 1/h. Also. making plots the data these tables. one may obtain 
stresses for intermediate values and a/h for either line loads 
strip loads. 

Under direction, Westergaard’s equations for stress 
thin slab resting liquid subgrade, and those Hogg for thin slab rest- 
ing elastic solid, were solved obtain influence charts and tables 
from which numerical values for stress may readily determined. Such 
charts and tables are found the above reference, and from them stresses 
for any general distribution axial loading may found. 

Fig. are plotted the horizontal stresses, ox, obtained from equations 
for the layered elastic solid developed herein (with rough interface), along 
with the flexural stresses obtained from Pickett’s numerical solutions for 
thin slab resting either liquid elastic solid subgrade. should 
noted that Pickett’s solutions for the slab resting elastic solid were for 
value Poisson’s ratio equal 0.4 only. The pertinent data for each sys- 
tem summarized follows: 


Two-layer elastic solid: 
P(concentrated line load) 
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ELASTIC SOLIDS 


Slab elastic solid: 


P(concentrated axial load) 
Slab liquid subgrade:3 


P(concentrated axial load) 


noted that the horizontal stresses, for the two-layer elastic solid 
and the slab resting liquid subgrade follow the same general trend, and 
are fairly close agreement 0.5, except that the origin, where 
the slab solution gives infinite stress. Beyond this point, the two-layer 
solid values are closer agreement with those for the slab resting the 
elastic solid. should also noted that the values flexural stress obtained 
for the slab resting liquid subgrade are tensile and consistently higher 
than those obtained from either the other solutions, least 
Since the factor safety for concrete pavement usually based the ratio 
the modulus rupture the concrete the maximum computed tensile 
stress, appears that pavement design based liquid subgrade would 
have somewhat higher factor safety against cracking than indicated the 
computation. 

the basis the numerical results this investigation, the following 
observations regarding the effects the pertinent parameters the stresses 
may made: 


Effects a/h stress: ox, oy, and increase steadily with increasing 
a/h all points when the surface loading uniformly distributed and con- 
stant When the loading line loading, that is, a/h 0,the 
horizontal stresses, and the interface are slightly higher 
general (except, course, the origin where they are infinite) than the cor- 
responding stresses for a/h=1. This might expected, because the total 
loads for a/h and are equal, and concentrated loads cause higher 
stresses the vicinity the load than distributed loads the same total 
value. also interesting note that the stresses for a/h and 
are identical four decimal places 10h, whichis fairly remote from 
the surface load. This accordance with St. Venant’s principle. 


Effects friction the interface stress: changing the conditions the 
interface from rough smooth, there usually slight decrease the 
stress) the upper layer, and considerable increase (compressive 
stress) the lower layer. points remote from the surface loading, say 
there only slight change ox, generally from slight tension 
slight compression. There only very slight change all points 
usually small decrease stress, and little change the value 
Txy, except the interface, where they all, course, become zero. This 
clearly indicates the pronounced effects friction has the horizontal stresses 


average value E/36 was used for the modulus subgrade reaction, val- 
commonly used the Texas State Highway Dept. The values for range from ap- 
proximately E/27 E/45, depending upon the type soil and other factors. theo- 
retical determination from the elastic constants, alone, not possible. 
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both layers, and hence the necessity taking friction into account the 
rational analysis design any slab pavement system. 

Finally, should mentioned that most interesting characteristic was 
observed for the smooth cases almost perfect flexural distribution the 
stress the upper layers the system. Ifa section the upper layer 
were considered beam depth point (Fig. would correspond 
point the upper surface the beam, and point would represent the cor- 
responding point the lower surface the beam, each point distance 
h/2 from the neutral axis. For pure bending, the flexural stresses these 
two points should, course, equal and opposite. Inspection Table 
Problem indicates values -6.8201 and +6.8189, respectively, points 


Value of O, 


O2h 04h 06h 10h 
Value of x 


INTERFACE TW-LAYER SYSTEM SUB- 
JECTED CONCENTRATED LOAD, FOR 
VARIOUS STIFFNESS RATIOS 


and which excellent agreement with this requirement. This flexural 
behavior most pronounced for low values where the underlying 
layer relatively flexible comparison the upper layer. 


Effects stiffness ratio: stress: For the smooth cases was 
found that increasing the stiffness ratio, from 0.02 0.1 nearly 
doubles the interface. the same time, however, (compressive) 
the surface loading virtually reduced factor two, are the hori- 
zontal stresses (tensile) the interface the upper layer. This would seem 
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indicate that increasing the stiffness ratio results less load being car- 
ried flexure, and more direct bearing the supporting layer, noted 
previously. 

For the rough cases, increasing the stiffness ratiofrom 0.02 0.1 increases 
approximately one-half, while, the same time, decreasing bya 
factor slightly less than three. The shearing stresses, the interface 
are, general, increased factor nearly three. 

Fig. shows the effect increasing the stiffness ratio, 
should noted that the stresses for the stresses 
base. noted that there less and less spreading the load with increas- 
mately one-quarter greater than the corresponding stresses for the homogene- 
ous case, and nearly three times greater than the stresses corre- 
sponding 0.02. This clearly points the pronounced effect pave- 
ment stiffness has subgrade vertical pressures. 


APPENDIX I.—NOTATION 
Ag, Ag, the displacement functions for horizontal dis- 
placement the upper layer 


the displacement functions for vertical displace- 
ment the upper layer 


constants the displacement functions for horizontal dis- 
placement the second layer 

constants the displacement functions for vertical displace- 
ment the second layer 

Young’s modulus elasticity 

modulus rigidity 

magnitude the concentrated line load, 

amplitude the uniformly distributed strip loading 

displacement function for the horizontal displacement 

displacement function for the vertical displacement 

width the loaded area 


depth Layer 
kj 1/1-2y,; 
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amplitude the sinusoidal loading 
spectively 
rectangular coordinates 
incremental normal stresses and directions, respec- 
tively 
Txy) incremental shearing stress the plane 
normal stresses and directions, respectively 
shear stress the plane 
Poisson’s ratio 


NOTE: Subscripts are used the text identify the layer question; thus 


Uj, Kj, refer the ith layer. 
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LATERAL STABILITY FRAMES ENERGY METHOD 


Donald Johnson,! A.M. ASCE 


SYNOPSIS 


method analyzing the elastic stability rigid frame structures sub- 
ject sidesway presented. The energy method applied assumed 
frame deformation compatible with its restraints. The assumed deformations 
are chosen that compatibility may achieved one standard moment 
distribution. Formulas describing the energies involved are given and pre- 
sented graph form. Two examples are given and compared results ob- 
tained more rigorous methods. 


INTRODUCTION 


well known that multi-story, rigid frame structures may buckle 
laterally sidesway under loads significantly less than those predicted in- 
dividual column-buckling formulas. Structures this category are usually 
multi-story frames without vertical truss systems. The methods available for 


rigorous solution the lateral stability such frames involve the task 


solving large numbers transcendental equations all but the simplest 
cases. Thepurpose this paper present simplified method that can 
easily used find accurate estimates the loads under which given frame 
will become unstable. 


Note: Discussion open until January extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Engineering Mechanics Division, Proceedings the 
American Society Civil Engineers, Vol. 86, No. August, 1960. 

Graduate Student, Mechanics Department, College Engineering, Cornell Univer- 
sity, Ithaca, N.Y. 
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Notation.— The letter symbols adopted for use this paper are defined 
where they first appear, the illustration the text, and are arranged 
alphabetically, for convenience reference, the Appendix. The symbols 
used herein (except Po, Lo, Ip, apply particular member which 
under consideration. member being considered, then the subscript 
implied but omitted simplify the equations. When two more members 
are under consideration simultaneously, subscripts are used avoid am- 
biguity. 


GENERAL ANALYSIS 


Assumptions .—Although the method presented could extended other 
cases, the problem limits itself follows: 


Frames considered consist horizontal and vertical members 
only. The configuration limited frames the beam-column slab- 
column type such shown Individual column heights need only 
identical for given floor and individual bay spacings need only identical 
when located aboveor below one another. Moments inertia need not iden- 
tical. Bottoms first floor columns will beconsidered clamped, although they 
can modified easily allow for various degrees restraint. All members 
are assumed uniform cross section. 

Only loads associated with lateral buckling sidesway, such shown 


Fig. are considered. Thus, joints each progressively higher floor 


undergo progressively greater displacements. Loads any member are as- 
sumed axial but need not identical. Only frame deformations due pure 
bending are considered. assumed that direct shortening lengthening 
any column due axial loads alone may neglected describing the buckled 
shape and computing the critical loads. Although this same assumption 
made the more rigorous methods analyzing the lateral stability multi- 
story frames, may result overestimating the buckling loads very tall 
and narrow frames that buckle modes shown Fig. Such frames are 
more accurately analyzed taking into account the direct shortening 
lengthening columns that occurs during buckling. 

Instability assumed occur the elastic range and consist 
small deflections. All assumed slender enough allow applica- 
tion the basic assumptions simple beam All joints are 
considered rigid. 


applied the entire structure, yielding the equation AT. the 
bending energy absorbed the structure and the total work done the 
external loads. Both quantities are due assumed buckling deformation 
the entire structure. 

The deformation each member assumed third degree poly- 
nomial variable which locates all points along the member’s length 
shown Fig. coefficients defining the polynomial for each mem- 
ber are expressed termsof the enddisplacements andslopes. The moments 
and energies related each member are, related the coefficients 


“Theory Elastic Stability” Timoshenko, McGraw-Hill Book Co., 1936. 
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ENERGY METHOD 


FIG, 1.—TYPE BUCKLING CONSIDERED 


FIG. 
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the polynomial. substitution, general formulas for end moments and 
energies for any member are derived. 

Derivation Deformation Equations.—The top end all columns will 
considered the right end. Hence thedesired deformation member (column 
beam) representable, general, shown Fig. Note that the total 
Since the energies and are related only the first and second deriva- 
tives the displacement, may neglected and only the relative displace- 
ment w(x) need considered. The assumed deformation w(x) chosen 
third degree polynomial 


wa) 


From Fig. can seen that differs from the amounte. Because 


Deformed 


Initial Position 


FIG. 3.—BUCKLING DEFORMATION TYPICAL 
MEMBER 


very small deflections are assumed, negligible relative tol Hence 
Introducing y(x) into Eq. yields 


differentiating Eq. gives 


d2w(x) 
dx 


similarly 


f 


= 


w(x) 


ENERGY METHOD 


Now have four unknown coefficients (A, and four boundary condi- 
tions: 


dw(x) 


dw(x) 


substituting the boundary conditions Eqs. into Eqs. and and solving for 


Derivation Moment and Energy have for the moment 


d2w(x) 
m(x) 
but the dimensionless measure m(x) M(x) Therefore 


Substituting Eqs. 6(c) and 6(d) into Eq. 7(b) and evaluating for +1, 
yields the following: 


(Note that the sign has been adapted moment distribution convention; 
clockwise rotation end member considered positive.) 


For horizontal beam members the relative displacements are neglected 

computing end moments from 8(a) and 8(b). This completely consistent 
with the assumptions made and explained the following manner. Consider 
frame its buckling deformation, such shown Fig. having lateral 
joint displacements each floor and slopes for each 
also increased Associated with the initial displacements 
are also displacements for each mem- 


(Sa) 


re: 
a 
a 
j 
43 


proportion virtue the equation 


(9) 


that applicable small deflections. Thus reduce the proportion 
slopes The relative displacement for horizontal 
beam member composed only the displacements the two adjacent 
rows columns below that beam member. Other displacements contributing 
are neglected virtue assumption number two. Hence 
computing end moments beam members. This often good 
approximation even the very small, because 
most frames the joints oor relatively level under buckling 
deformations. 

For beams 


= 

> 
+ 
D 


expressed terms the dimensionless quantities and 
m(x) 
(11b) 


substituting into values and from 6(c) and 6(d), perform- 
ing the integration and introducing 


(12) 
yields 


\ 


Eq. 12(b) only applicable for that is, the member contains 
inflection point. Values versus are presented graphically Fig. 
for values are required beyond this range, Eqs. and 13(b) 
and Fig. can applied facing the member from the opposite direction. 
The value represents the portion energy member being summed. 
formula found: 


This expression much more cumbersome than Eq. 13(b) but applicable 
the member contains inflection point. 

The term the energy equal the quantity for given member. 
The sum over all members the energy released all the external 
loads. may expressed 


for small displacements. 
Expressing its dimensionless form 


Substituting values from Eq. and performing the integration 
yields 


Values for various values and are presented graphically 
Fig. Because quite constant for ratios and used 


sidesway, this graph can used with sufficient accuracy. 
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FIG, 6.—EXAMPLE 


(2) (5) (8) 
(1) (9) (10) 


General Choose convenient arbitrary lateral 
displacements Wr]. Wr2. Wr3. ... for each floor that could most probably 
describe the true buckled shape. Although the quantities 
must actually much less than described the energy 
and moment formula applicable small deflection theory, should noted 
that any convenient series numbers can used and interpreted one 
magnification factor times series small fractions. can shown that 
both and would then proportional the square this factor, which 
would cancel the equating the two. These numbers should selected 
with good however. any reasonable guess will give excellent re- 
sults. suggested that floor displacements chosen simplify the 
minimum. 

Select, for each member end slopes and end moments 
that are compatible with the assumed lateral deflections Wr3, 

all joints equilibrium. For simple cases these can 
found solving equations representing moment summations each joint 
based Eqs. The functions w(x) have been chosen that they are 
identical functions that describe beams bending under the influence end 
moments only. the desired compatibility can obtained use one 
application moment distribution. The dimensionless fixed end moment in- 
duced each column then 


After the moments are completely distributed, slopes joints may found 
proceeding from known slope (usually bottom) using the moment 
area relation which reduced 


or. more simply 


for any member containing inflection point. This formula correct sign 
for use proceeding from the left the right member. Itis recom- 
mended that values all slopes checked against final end moments use 
Eqs. and 10. 


Compute for every member and and for eachcolumn. 
Using graphs for and the values and for each member. 
yields: 
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Hence the buckling coefficient 


Either Eq. may used determine the cricial load. The energy 
summation extended over the entire structure over any represen- 
tative portion for which known that the ratio identical that 
the entire structure. 


EXAMPLES 


Example I.—This problem involves the frame shown Figs. 6(a), the as- 


sumed buckling shape, and 6(b) the givenframe. The solution obtained 
follows 


Lo, and are chosen according Fig. 6(b). 

The assumed deformation chosen shown Fig. 6(a). Notice that 
all joint rotations each floor are the same, and @9, respectively. ad- 
dition, the columns above each floor have identical deformed shapes that 
all loads each floor descend the same distance. other words, 
and =eg. Each beam deformation antisymmetric and its in- 
flection point shown hinge Fig. 6(a). can shown applying 
Eqs. 13(c) and 15(a), similar intuitive reasoning, that, for the deforma- 
tion assumed, the portion the frame the left section Z-Z has the same 
ratio energies the entire structure. 

reasonable assume greater than because the bottom the 
structure clamped. 

The joint equilibrium the portion the frame the left section 
Z-Z can analyzed contained the hinges shown Fig. 6(a) because 
the antisymmetry. Hence modified stiffness factors are used and the mo- 
ment distribution performed only the portion shown below (see also 
Table 2): 


For the computation factors use the general formula 


Table obtained. 
For the computation fixed end moments Eq. applied obtain 


6(2) (8) 


+96 


4 
and the formula for member with hinged end 
41, 
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The deflected slopes the two joints and are determined proceeding 
from zero slope the bottom using Eq. 18(b): 


and 


Checking the beam moments from Eq. 10(a), noting that the deflection anti- 
symmetric about Z-Z, yields 


91.2 vrs 90.9 from moment distribution 


and 


6(2.32) 41.7 vrs 40.8 from moment distribution 


This accuracy sufficient because the general method gives results only 
within 2%. 


TABLE 


Member 


Joint 


Joint 


Computations energies using graphs Figs. and are presented 
Tables and 


This identical problem has been previously computed two rigorous 
methods, one using the method four moment equations and the other em- 
ploying the method slope deflection equations, Mr. John 
These methods both yielded 7.09. Hence the error induced using the 


Elastic Stability Rigid Frame Structures Subject Sidesway” 
Master’s Thesis, Lehigh University, 1958. 
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ENERGY METHOD 


TABLE 2.— TABLE BENDING ENERGY 
MOMENT DISTRI- 
BUTION FOR 
EXAMPLE ONE 


Only energy the left section Z-Z summed. 


TABLE 4.—SUMMATION LOAD ENERGY 


60.0 
-10.9 -49.1 (2) 394 
-20.8 |-41.6 |-93.6 
-5.4 
\Fo 
96.0 
75.9 


energy method this case very small. Results obtained from the energy 
method for assumed lateral joint displacements and 
energy method lateral displacement ratios varying much 40% 
from each other, yield critical loads varying less than from each other. 

Since these computations involve chart and slide rule accuracies common- 
related errors the order 1%, should noted that any reasonable 
selection and would have given satisfactory results. 

Example example involves multi-story frame similar Fig. 
with the following features: 


Total load all columns 

The frame consists large number stories and bays. Thus the 
problem analogous that formulated for analysis tall wide 
frames. 


Solution.—1.In order avoid the unnecessary work distributing mo- 
ments and summing energies over many members, minor modifications 
boundary conditions similar those shown Bleich will made. Fig. 
shows the modifications boundary conditions made the present analysis. 

The solution presented here not limited the size frame shown 
Fig. but applies any large system with this type boundary modification. 
The reductions moments inertia and loadings the top and side edges 
the frame, shown Fig. can seen have little influence the total 
critical load large sized structure. further simplify the problem, 
however, necessary reduce the bottom restraints that given set 
beams, also shown Fig. that renders the loaded structure symmetrical 
two directions. With these modifications, each column and beam are as- 
sumed deflect the same exact shape, can shown that moments all 
joints will balance. Furthermore, can seen that the energy ratio 
for, any representative area, shown dotted lines Fig. equal 
the same ratio for the entire structure. The deformation this representa- 
tive area shown Fig. Because the repetition deformed shapes 
throughout the frame, the energy summation need only summed over mem- 
bers and BE. Thus thecritical load can determined considering only 
the energies one interior beam and column together. 

The moments any internal joint shown Figs. and9 are balanced 
(Note that since all interior columns and beams have same 
deflected shape and their moments are the same throughout the structure.) 


yields: 


“Buckling Strength Metal Structures” Bleich and Bleich, McGraw Hill 
Book Co., 1952. 
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FIG. BOUNDARY MODIFICATIONS FOR EXAMPLE 


FIG. 8.—ASSUMED BUCKLING DEFORMATION 
REPRESENTATIVE AREA ADCE 


x 


uo 


Since can any convenient number, let then 
From 13(c) 


=6a+6+6 (1 + a@)2 - 12 (1 + @) 
30 


Lo Ip 


Table shows the comparison between values computed from Eq. 
and those given Bleich for 10-story building for various values 


2 
energy method From Bleich 
1/4 2.39 2.32 
1/3 2.97 2.91 
5.71 5.69 


10.00 9.87 


Bleich’s values have been converted from effective lengths buckling coeffi- 
cients. Numbers were scaled from graph. 


CONC LUSIONS 


addition the assumptions described the section general analysis, 
the following limitations, and/or extensions, should considered: 
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ENERGY METHOD 


Pin restraints bottom columns can described letting 
Other boundary conditions can represented use Eqs. and 10. 

The method can extended describe failure the plastic zone 
trial anderror. First estimate Then determine the tangent modulus 
for each column for the estimated use stress-strain diagram. 
The formulas presented herein are then applicable each column re- 


duced Apply the energy method using the reduced values yield 


accurate Continue the procedure using the more accurate 
and on. 

Inspection Eq. and the preceding derivation reveals that for giv- 
the values the factor are proportional Thus see from 
Fig. that the value for any column relatively independent its end 
slopes and almost entirely dependent Although may seen from 
Fig. that variations end slopes may cause variations high 20%, 
this extreme not realized many frames. frames having the most pro- 
nounced variation member sizes the vertical values the 
joint rotations along any particular floor will quite similar. Consequently 
the row columns above any particular floor will have very similar end 
slopes and very similar values the total load energy 
would approximately the same regardless the way which the total load 
each floor divided among the columns. This fact might used allow 
investigator move loads from column column along each floor 
simplify the mathematics the more rigorous methods without great 
loss accuracy. 

metric buckling mode deformation that compatible with the joint equi- 
librium requirements, even the loading isnot symmetric. Such modes allow 
the investigator confine the problem smaller “representative portions” 
similar those used the given examples and thereby reduce the office 
work required distribute moments and sum energies. 

The method could adapted frames, other than those shown 


that contain only vertical and horizontal members. such adaptations 
are beyond the scope this paper. 


This method intended for use checking the lateral stability frames 
already designed. Low degree polynomials have been used the past solve 
with accuracies few percent many stability problems concerning one 
member. The method has been developed that slide rule and well known 
moment distribution methods can achieve results quickly and easily within the 
usual required accuracies. This method requires little additional work par- 
ticularly designs which lateral moment distributions have already been 
made during design for wind, and forth. 

The method can used structures with varying column and beam sizes 
and spacings, varying types stiffeners, and on; cases which more 
rigorous approaches are often unfeasible. The designer’s intuitive knowledge 
sidesway buckling modes can utilized yet need not exact order 


yield good results. 
ADDITIONAL READING REFERENCES 
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APPENDIX. NOTATION 


length member 

L/2 half length member 

reference length, arbitrarily chosen 

half reference length 

moment inertia member 

reference moment inertia, arbitrarily chosen 

axial load member 

reference axial load, arbitrarily chosen 

variable locating points along member (see Fig. 
function describing the assumed buckled shape (see Fig. 
coefficients polynomial used describe w(x) 
value w(x) right end member (see Fig. 

w,/L dimensionless measure 


relative displacement between ends member axial direction 
(see Fig. 


slope buckled shape right end member 
slope buckled shape left end member 
moment right end member 


moment left end member 


energy due bending member 


Lo2 
buckling coefficient the frame 


absolute value 
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parameter locating inflection point 


proportion energy member being summed 
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SEEPAGE LOSSES FROM PARALLEL CANAL SYSTEMS 


Hammad! 


SYNOPSIS 


This paper deals with the two-dimensional problem steady seepage flow 
under gravity from system parallel, identical, and equally spaced canals 
into semi-pervious clay layer finite thickness underlain freely per- 
meable layer sand and which the piezometric headis very near the 
canal water level. Two steps conformal mapping are used and approxi- 
the canal profile adopted. this approximation the specified canal 


width and depth are not changed. 


INTRODUCTION 


most agricultural zones, especially river valleys, the stratified soil 
formed top clay layer finite thickness and medium low permea- 
bility, underlain permeable layers sand and gravel. The piezo- 
metric head the lower, freely permeable layers controlled the river 


running the valley. 


If, for any reason, the canals running the top clay cap are run higher 
water level than the piezometric level the lower sand and gravel layers, 
steady seepage flow will take place from the canals through the clay layer 


the sand and gravel layers beneath. 


Note.—Discussion open until January 1961. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Engineering Mechanics Division, Proceedings the 


American Society Civil Engineers, Vol. 86, No. August, 1960. 


Asst. Prof. Applied Mech,, Faculty Engrg., Alexandria Univ., Alexandria, 
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the present treatment, discharge formula obtained for the seepage 
loss from system parallel, identical, and equally spaced canals, which 
the water level slightly higher than the piezometric head the sand and 
gravel layers beneath. This last condition contrast with Wedernikow’s 
problem2 deep-lying watertable the lower sand and gravel aquifer. The 
same problem, limited single canal, has been attempted the writer.3 
simple seepage discharge obtained for the present problem, from 
which the single canal problem derives special case when the spacing 
the canals assumed tend infinity. 


GOVERNING EQUATIONS AND BOUNDARY CONDITIONS 


may found texts flow through porous the governing dif- 
ferential equation the flow that Laplace: 


given 


which the permeability the porous medium, the pressure any 
point, the fluid density and the acceleration due gravity. 
According these assumptions, the boundary conditions are: 


separation between the semi-pervious clay the lower 
freely permeable layer sand and gravel equipotential line. 

The canal profiles are also equipotential lines. 

The vertical lines symmetry are streamlines, because there can 
flow across them. These lines may, fact, replaced solid bounda- 
ries. 


Thepressure vanishes water surface, which also astream- 
line. 


METHOD SOLUTION 


Because the piezometric head the underlying sand assumed very near 
the canal water level, the free water surface the clay layer will, under steady 
conditions, also very near the canal water-level. Consequently, the free 
water surface may approximated substituting horizontal line coinciding 
with the canal water level. This extra the problem will cause 
slight exaggeration the seepage discharge obtained this solution. 
may noted, however, that the deviations the true free surface from the 
substituted horizontal line are greatest points A.B.C..... (Fig. which are 

“Flow Homogeneous Fluids through Porous Muskat, Edwards, 
1946, 


Losses from Irrigation Canals,” Hammad, Proceedings Paper 
1991, ASCE, April, 1959. 
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SEEPAGE LOSSES 


FIG, CONFORMAL MAPPING 


| 


stagnation points the flow field and, consequently, the approximation willre- 
sult serious error. 

first order approximate solution the problem series 
parallel shallow canals (depth semi-breadth and spacing center 
center, may obtained considering the seepage under gravity from 
infinite number breadth andspacing inthe solid boundary 
shown Fig. 2(a). Because the lines *L, are lines 
symmetry, they may considered solid boundaries that will divide the 
field flow into series identical panels. will then sufficient con- 
sider the flow one these panels, Fig. 2(a). 

streamlines, whereas the line equipotential. 

has already that water surface the seepage layer 
high and very near the canal water level. Consequently, the original as- 
<L) and(y have been taken the top streamlines instead 
the actual free surface that very near them. 

conformal transformation using the Schwarz-Christoffel theorem each 
step. The term the velocity-potential, refers the stream-function, and 
the usual notation. The two steps are follows: 

Step 1.—Mapping conformally the rectangular panel bounded 
and Lin the Z-plane onto the lower half the -plane, Fig. obtain, 
terms complex functions, 


The points having the same numbers both planes correspond one another 
and constant. Integrating Fig. produces 


fore 


(5) 


terms Jacobi’s elliptic functions 


the first kind corresponding the complementary modulus 
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The result the first step conformal mapping then, 


Step 2.—Following the same steps mapping, conformally, the rectangle 
(Fig. derived terms Jacobi’s elliptic functions, 


plete elliptic integral the first kind corresponding the complementary 


modulus The result the second stepof conformal mapping 


The two relations, Figs. and give the relation the flow. 

evaluate the constants and terms the dimensions the 
flow panel and (Fig. 2), proceed follows. 

point (Fig. Z=b, w=a, andat point Z=L+iD, 
Substituting these two points Eqs. and produces 


Using the properties elliptic functions Eq. set equations may 
reduced the following useful equations: 
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Because practice, the clay thickness, small compared with the 
canal spacing follows from 11(c) that small compared with 
Consequently, the modulus small. The complete elliptic 
functions and expanded convergent series are4 


first order approximation may put 


From Eq. 
and, 


The resulting values and may compared with their more exact 
values the mathematical tables the complete elliptic functions. Compu- 
tations show that this approximation quite sufficient for all practical values 
D/L. 

introduce the dimensions canal section represented the breadth 

Introducing these values into Eqs. and 


For small values and Eq. further reduced, applying the: 


properties elliptic functions Eq. the following simple relations: 


“Function Complex Variable,” Th. Macrobert, Macmillan and Co., Lon- 
don, 1925. 


4 2 


tan 
> 1 tan (34) = (16b) 


tanh 


which 


SEEPAGE DISCHARGE FORMULA 


Having evaluated the unknowns terms the known values 


the flow panel and (Eqs. 11, and 16) one may proceed obtain 
general formula for the seepage loss per unit length each canal. 

From argument given Terzaghi,° the seepage flux out each unit 
length canal, normal the plane motion, (Fig. is, accordingly, given 


which the head loss between the canal boundary and the line 
the difference between the canal water level and the piezometric head the 
freely permeable sandy layer, and the permeability the clay cap. 


Substituting the values for and obtained from Eqs. and into Eq. 


yields 


which and are the complete elliptic integrals the first kind whose 
moduli are and kg' respectively given Eqs. 17. 

may worthwhile mentioning that the case drain which the 
water level lower amount than the piezometric head the sand and 
gravel aquifer, the seepage flux will reversed and the drain will gaining 
rather than losing the rate given Eq. 19. 


“Theoretical Soil Mechanics,” Terzaghi, John Wiley and Sons, New 
pp. 246, 


and its complementary 
2 
D'-c 
att 


SPECIAL CASE 


The special case single canal attempted the writer separate 
paper, may derived from the present result making the canal spacing 
tend infinity. that case the moduli the complete elliptic functions 
given Eq. will reduce 


the discharge formula given Eq. same form but the moduli 
its complete elliptic functions must taken from Eq. instead Eq. 17. 


NUMERICAL EXAMPLE 


For parallel canal system spacing depth each canal 
head the sand and gravel aquifer, the following relation would true: 


tan 


0.07830 


From Fig. 17, 


51 


Thus, 0.00613 and the square its complementary modulus 0.99387. 
From tables elliptic functions the corresponding complete elliptic inte- 
grals the first kind are 1.5733 and 3.9385. 


The seepage loss per unit length each canal may then computed from 
Eq. 


1.5733 


0.0147 gpm per length canal. 
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VIBRATIONS STRUCTURAL SYSTEMS 
COMPONENT MODE SYNTHESIS 


Walter Hurty! 


SYNOPSIS 


Natural modes and frequencies structural systems are determined 
energy method using mode functions applicable the complete system, 
subsystem thereof, that are synthesized from admissible se- 
lected for the component members the system. The synthesis accom- 
plished using equations constraint that follow from conditions force 
equilibrium and deflection compatibility the junctions. 


INTRODUCTION 


The analysis vibration dynamical behavior complex struc- 
tures has increasingly occupied the attention analysts during recent years. 
Methods analysis applicable simple structural elements have been known 
and used for wellover century and these methods have been applied ideal- 
ized models complex structures. Examples are found analyses multi- 
story building structures and high aspect-ratio airplane wings which these 
structures are dealt with ideally beams. Such analyses may quite accu- 
rate and useful long the bases for the idealizations are valid. 

more recent years methods have been devised which the behavior 
structure may predicted interms the properties its elements. These 
methods analysis involve, explicitly, every element the structure, hence 
they enable one determine the effects changes these elements the 
behavior the entire structure. The method analysis presented this 
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paper this category and the words “Structural Systems” the title are 
used imply systems analysis this sense. 

order clarify the position this study relative other structural 
systems studies, pertinent refer some those that have appeared 
the recent literature. method analysis considerable accuracy and 
versatility has grown out the well-known methods Holzer,2 Myklestad, 
and Thomson.4 has been generalized Pestel and and 
described method transfer matrices. applied the determination 
natural modes and frequencies, characterized the setting 
frequency determinant that follows consequence satisfying boundary 
conditions. The computations, which are conveniently carried out matrix 
operations, require the insertion trial values the frequency into the 
transfer matrices. Since many such trial values may necessary order 
identify the frequency spectrum, great amount computation generally 
involved and itis conceded that high-speed electronic digital computing facili- 
ties are required order that the method practicable. The versatility 
the method lies inthe fact that response static loads well steady-state 
forced vibrations may handled readily. 

Another powerful method analysis has been applied frames 
Described method receptances, deals with the forced 
vibrations systems interms the receptances its component parts. The 
convenience this method lies the fact that these individual receptances 
can catalogued for ready reference. For determining natural modes and 
frequencies the structure, those frequencies are sought for which the over- 
all receptance vanishes. This procedure also requires trial-and-error so- 
lution. 

Veletsos and have devised method described exten- 
sion Holzer’s method for continuous beams. also applicable frames 
without sidesway, that is, frames which the junctions corners not 
translate. For natural frequency determination, this method seeks those fre- 
quencies for which exciting couple applied the beam vanishes. This, 
the methods already mentioned, requires repeated trial computations. 


Die Berechung der Drehschwingungen, Holzer, Springer-Verlag, Berlin, 1921. 
Vibration Analysis, Myklestad, McGraw-Hill Book Co., 1944, 


“Matrix Solution for the Vibration Nonuniform Beams,” Thomson, Jour- 


nal Applied Mechanics Trans. ASME, Vol. 17, No. Sept. 1950, pp. 337-339. 
“Ein allgemeines Verfahren zur Berechnung frier underzqungener Schwingungen 


von Stabwerken,” Eduard Pestel, Abhandlungen der Braunschweigischen Wissen- 
chaftlichen Gesellschaft., Vol. VI, 1954, pp. 227-242. 


“Katalog von Ubertragungsmatrizen zur Berechnung Technische Schwingungspro- 
bleme,” Pestel, Schumpich, and Spierig, Vol. 35, 1959, pp. 
11-28, 


«The Analysis and Synthesis Vibrating Systems,” Bishop, Journal 

the Royal Aeronautical Society, 58, Oct. 
“The Analysis Vibrating Systems which Embody Beams Flexure,” 
Proceedings the Mechanical Engineers, Vol. 169, 1031, 1955. 


Mechanical Engineers, Vol. 170. 955, 1956. 


“Natural Frequencies Continuous Flexural Members,” A.S. Veletsos and 


Newmark, Proceedings American Society Civil Engineers, Vol. 81,1955, Paper 
735, 


Determination the Natural Frequencies Continuous Beams Flexible Sup- 
ports, Veletosos and Newmark, (to published). 
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contrast with the above methods, the procedure described this paper 
leads the formulation eigenvalue problem terms generalized co- 
ordinates and generalized inertia and stiffness coeificients that are applicable 
the complete structure. Once this problem formulated, direct solution 
obtained one several possible techniques that are well-known. The 
generalized properties are obtained through the application energy methods 
using sets admissible mode functions for the various elements the sys- 
tem. Best approximations the true mode shapes are obtained superposi- 
tion finite number these functions. This method, applied simple 
structures idealized models complex structures, has been used widely 
and often called the method. applied structural sys- 
tems, has often been regarded inaccurate because the difficulty de- 
fining suitable mode However, the present analysis, mode shapes 
applicable tothe entire system are not defined the outset. Instead, they are 
synthesized from mode shapes that are selected for the separate elements 
the system. Their synthesis accomplished technique that results from 
the application equations deflection compatibility and force equilibrium 
the junctions these elements. The delineation this technique con- 
sidered the principal contribution this paper. addition, the tech- 
nique extended enable one synthesize modes and generalized proper- 
ties for subsystems that may then put together form the complete sys- 
tem. The computations required for each the subsystems may proceed inde- 
pendently. Points similarity exist between this analysis and that Saibel 
and D’Appolonial3 their treatment continuous beams. that analysis 
mode functions are synthesized from the eigen functions simple beam ob- 
tained removing the intermediate supports. 

Notation.—The letter symbols adopted for use this paper are defined 
where they first appear, the illustrations the text, and are arranged 
alphabetically, for convenience reference, the Appendix. 


ENERGY ANALYSIS THE SYSTEM 


For this analysis shall consider system structural members such 
beams and rods connected their ends. These connections may pins 
rigid joints combinations these, and the various members may sub- 
jected bending, torsion, axial loading combinations these loadings. 
shallrestrict the motions small amplitude sothat free vibrations 
the separate elements longitudinal, transverse and twisting motion are 
uncoupled. principle, the geometry the structure immaterial and the 
individual members need not uniform. The complete structure will sub- 
ject boundary conditions that may defined one more the junctions. 

The first task the analysis select mode functions for each member 
and superpose these functions define longitudinal, transverse and tor- 
sional vibrations throughout the system. write any one these modes 


Aeroelasticity, Bisplinghoff, Holt Ashley, and Halfman, Addison- 


Wesley Publishing Co., Inc., 1955. 
“Forced Vibrations Continuous Beams,” Edward Saibel and Elio D’Appolonia, 
Proceedings the ASCE, Vol. 77, No. 100, Nov. 1951. 
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vibration interms the selected functions and set coordinates follows: 


which the longitudinal displacement, reférs transverse displace- 
ment, denotes torsional displacement, the jth mode function, indi- 
cates the independent variable measured along length member, and the 
jth coordinate. Although not explicitly noted, are functions time 
well and the p’s are functions time. indicated functions and 
coordinates are chosen for the structure and these are distributed among the 
various members. Some the functions represent longitudinal motion and 
others transverse and torsional motion. 

Although wide choice exists the selection functions several condi- 
tions must exist. For example, they must satisfy the boundary conditions 
the structure whole, that is, given element the structure subject 
such boundary conditions then the functions chosen for that element must 
satisfy them. Conditions constraint the junctions two more mem- 
bers are not considered boundary conditions but will termed junction 
conditions. general, functions must chosen which admit displacements 
the junctions and also admit the existence forces the junctions. For 
example, member enters junction and that member may have 
transverse displacement that point, then functions must chosen for 
which transverse shear force may also exist that point, then 
functions must included for which where the prime notation 
used denote differentiation with respect Functions which satisfy the 
necessary conditions are termed admissible functions this study. 

Limited experience problem solutions gained thus far indicates that 
simple functions may chosen. the examples included this paper inte- 


gral powers are used with good results. Thus, displacements are 
represented the series 


the case transverse displacement, the functions 1(x) and represent 


rigid-body modes translation and rotation, respectively, and the remaining 
functions represent deformation modes. longitudinal displacement the 
function 1(x) represents rigid-body translation. The use such rigid-body 
modes permits the inclusion joint translation the vibration modes the 
complete structure. For example, sidesway building frames may handled 
this technique. 

Consideration the junction conditions leads equations deflection 
compatibility and force equilibrium which serve equations constraint 
among the coordinates Therefore, these are not general- 
ized coordinates the Lagrangian sense. Application Hamilton’s principle 
this problem leads the equation!4 
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which the kinetic energy the system, refers the strain energy 
the system, indicates the variation the coordinates pj, and the 
time. The dot notation used denote differentiation with respect time. 
Since the p’s are not independent, the variations 5p; cannot arbitrary, 
hence, the expressions the braces Eq. not vanish separately. The 
equations constraint permit express the coordinates terms 
set independent generalized coordinates thus, 


which the are coefficients the linear transformation. Note that 
equations constraint exist, then 


n 


Substituting this into Eq. and taking account the independence the vari- 
ations obtain 


For small amplitude vibrations conservative system, the terms Eq. are 
given the following equations. 


m 
— = 0 
m 


these equations the are inertia coefficients and the co- 
efficients associated with the coordinates For translation, either the lon- 
gitudinal transverse direction 


where the mass the element per unit length point and the integra- 
tion taken over the length the member. For transverse bending 


vers 

‘ 
k=1 

i=1 
ie 
2 
| 

ag 

| 
= 

a4 


where the modulus rigidity flexure. For longitudinal strain and 
torsional strain, the stiffnesses are given, respectively, 


and 


The quantities and are moduli rigidity tension-compression and 
torsion, respectively. 

may substitute Eq. into Eq. and write the set equations ma- 
trix form. 


Here, the transpose the transformation matrix [C]. now sub- 
stitute Eq. which, matrix form 


and note that for harmonic vibrations 


where the angular frequency the vibration. then obtain the equation 


This equation expresses the desired eigenvalue problem whose solution yields 


should noted that the matrices [m] and [k] are square symmetric ma- 
trices order having elements the inertia and stiffness coefficients 
such obtained from Eqs. through 12. From Eq. seen that the 

passing, should noted that the effects rotatory inertia and trans- 
verse shear strain may included the analysis introducing the appro- 
priate terms the computation the mjj and the kjj. 


EQUATIONS CONSTRAINT 


this section shall concerned with the equations constraint among 
the coordinates that lead the transformation Eq. 14. These equations ex- 
press the compatibility displacements and force equilibrium the junctions. 
The terms displacement and force are usedinthe general sense include ro- 
tation well translation and moment well force. 

clarify the development these equations let consider simple junc- 
tion two perpendicular members shown Fig. 1(a) with forces and displace- 
ments only the plane the members. 

Fig. 1(b) shows section the junction with positive axial forces, posi- 
tive shears, and positive bending moments, Fig. 1(c) positive dis- 
placements and are shown relative the positive direction for 
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(a) FIRST MODE (b) SECOND MODE 


FIG. 1.—A SIMPLE JUNCTION 


SUB-SYSTEM 


FIG, 3.—BENT CANTILEVER BEAM 


(a) (b) (c) 
FIG, TWO NATURAL MODES BENT CANTILEVER BEAM 
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either member. The angle represents rotation due bending and, the 
absence shear deflection, equal 
The displacement equations are 


Wy (11) ug(l2) = 68 . (17) 

The force equilibrium equations are 
Ny + Va2(l2) = 0 


clear that comparable equations can written for members intersecting 
other angles. plane system with forces and displacements its plane 
six equations exist for two members. For each additional member the junc- 
tion three additional equations displacement compatibility exist. 
general, for members junction, equations constraint exist that 
junction, which 


For junction 3-space, can shown that 


Upon introducing the equations constraints become linear homogene- 
ous, algebraic equations the p’s which may expressed 


where the matrix [A] order Now, the p’s, may se- 
coordinates interms which the remaining maybe expressed. 
These latter shall call constrained coordinates, Then Eq. 
written 


The matrix [A]; order and the free coordinates must chosen 


leads 


This equation may used construct matrix which column 
taining all the coordinates may expressed terms the free co- 
ordinates Thus, 


The matrix [C] order and the constraint transformation ma- 
trix Eq. 14. the construction carried out this manner, the col- 
umn {q} Eq. identified with column Matrix [C] used Eq. 
set the eigenvalue problem for the complete structure. 
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STRUCTURAL VIBRATIONS 
PARTITIONING INTO SUB-SYSTEMS 


Eqs. and show that for systems with more than few junctions the 
number equations constraint may very large and the computations re- 
quired multiply and invert large matrices tend become tedious. There- 

fore, variation the foregoing procedure offered for constructing 
which permits dealing with subsystems. 

Consider portions two subsystems and that connect corner 
shown Fig. For clarity shall consider these subsystems lying 

common plane and shall consider only forces and displacements that plane 
Interacting forces exist between these subsystems. The forces Mj, Xj, 
shown Fig. are those imposed subsystem subsystem The equa- 
tions force equilibrium among the members subsystem this junction 
are not now homogeneous since these force terms must included the 
right-hand side the equations. this case Eq. takes the form 


which the partitioned column the right consists zeroes and 


three force terms (Mj, Xj) the plane case. The subcolumn may take 
either two forms 


\ 


where and may any characteristic values. The purpose intro- 

ducing these quantities this manner give the elements {F} the di- 

mensions length that the elements matrix [A] may dimensionless. 


The partitioned matrix order The upper part null matrix 


order and the lower part, the identity matrix. The right side 
Eq. may arranged give the following form: 


The orders the submatrices are given. before, matrix may in- 
verted and may express the constrained p’s terms the free p’s 
and the forces acting the subsystem. 
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now possible construct the transformation matrix using the follow- 
ing arbitrary, but convenient, procedure. First set the top element 


equal unity and all other elements zero. Using Eq. compute the cor- 
responding elements Thus, obtain column containing all 
the constrained p’s and chosen free p’s. Next, set the second 


element unity and allothers zero. Thus, second column 


obtained. When carry this completion clear that may obtain 
such columns. Matrix for subsystem then constructed from 
these columns follows: 


procedure arbitrary since any set combinations the elements 


may used construct the columns [C];. However, must noted that 
any such selection elements must linear combinations those 
chosen the suggested procedure, more than independent columns 
may had any event. 

The procedure for constructing [C]; may thought alternately, 
procedure for synthesizing admissible functions for the connected subsystem. 
These functions are synthesized from the functions chosen for the indi- 
vidual unconnected elements. Eq. substituted into Eq. obtain dis- 
placements the connected system terms the free coordinates, Thus 


j=1 
which 


one the synthesized admissible functions. the case the subsystem, 
with three interacting forces one junction these functions are 
obtained, three which serve admit the existence those forces. the 
case the complete system, only admissible functions are obtained, and, 
thus, only degrees freedom exist for the complete system. Thus, may 
seem that gained three degrees freedom the subsystem. Such not 
the case, however, since these apparent degrees freedom are lost com- 
bining subsystems into the complete system. 


COMBINING SUBSYSTEMS 


Rather than attempting treat this problem complete generality, 
shall only consider the problem combining two subsystems one junction 
such illustrated Fig. The plane system again considered. be- 
lieved that this illustrated will serve indicate procedures followed 
more general cases. 
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STRUCTURAL VIBRATIONS 


Again, equations constraint must derived compatibility 
displacement and force equilibrium the junction the two subsystems and 
From the results the last section, may determine the forces and dis- 
placements terms the coordinates thus 


and 


which column matrix with six elements the order 


For the system shown, the two matrices [B] are order (nj and 
(nj 3), respectively. Equating the two deflection-force vectors the junc- 


tion are led constraining relationship among the elements the two 
coordinate vectors 


From this construct new constraint matrix that relates the sets 
q’s pertaining both subsystems minimum set free q’s. Since the 
two sets, have totalof coordinates and since have six equa- 
tions constraint follows that the complete system composed the two 
subsystems has free generalized coordinates and, hence, that num- 
ber degrees freedom. Generalized inertia and stiffness matrices for the 
complete system are found operating the corresponding matrices for 
the unconnected subsystems bythe constraint matrix shown Eq. 16. 


EXAMPLE 


Bent Cantilever Beam.— The first example cantilever beam length 
bent through right angle its center shown Fig. wish de- 
termine the first two natural modes and frequencies vibration the plane 
the beam. This simple example chosen illustrate the method ap- 
plied complete system. Bishop worked9 this problem the exact method 
receptances, hence shall have standard against which check the 
present method. The boundary conditions are: member (1); 


order compare with Bishop’s solution, shall postulate infinite axial 
rigidity both members. Therefore, two equations force equilibrium 


the corner are not included the equations constraint. The four required 
equations are 
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The functions selected follow. For member (1) the transverse displace- 
ment functions are and For member (2) the trans- 


and the longitudinal displacement function may verified 
readily that these functions satisfy the boundary conditions and that they are 
admissible functions. The uncoupled inertia and stiffness matrices are obtained 
using Eqs. and 


The four equations constraint are 


select and thefree coordinates and determine the matrix [C] tobe 
follows: 


4.50 
0.33333 0.5 
Following the operations Eq. the eigenvalue problem becomes 


frequencies inthis problem than Bishop’s 0.5% and 1.0% respec- 
tively. The two natural modes are shown Fig. 
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Two-Story Building Frame—Anti-Symmetrical Modes.—This example 
added illustrate the procedure for handling subsystems. The two-story 
building frame shown Fig. 5(a) divided into two parts. The first story 
(subsystem shown Fig. 5(b) and the second story (subsystem Fig. 
5(c). The structure has vertical plane symmetry that only one-half 
considered, and only the antisymmetrical modes are determined. For 
convenience computation, all members are considered identical. The 
boundary conditions for the antisymmetrical modes are for system 


Member (1): (0) 
w"(0) 
(0) 
Member (2): (0) 
(0) 
(0) 


for system 
Member (1) (0) 
w"(0) 
(0) 


The equations displacement compatibility and force equilibrium the cor- 
ners the two subsystems are follows: for system 


for system 


Functions selected are listed for the two subsystems separately. Since they 
are dealt with independently attempt made distinguish between the two 
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sets numbering them. For system 


For system 


Note that functions through are identical for the two systems but that 
and are added for system admit rigid body motion that 
system that cannot occur the lower frame. 

Inertia and stiffness matrices are computed for the uncoupled elements 


the two subsystems, using Eqs. 10, and 11. The ratio axial transverse 
stiffness 


1600, 


which the section radius gyration bending. These matrices are 
those the first example. The constraint matrices for the 
two subsystems are found the procedure leading Eq. and will indi- 
cated subsequently connection with the equations constraint show the 
ordering the coordinates. 
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STRUCTURAL VIBRATIONS 


MEMBER 


SUB-SYSTEM MEMBER 


SUB-SYSTEM 


MEMBER 


MEMBER 


FIG, 5.—TWO-STORY BUILDING FRAME 


FIG, MODES FRAME STRUCTURE 
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For System For System 


The generalized inertia and stiffness matrices for the two subsystems are 
computed from the equations 


using the inertia and stiffness matrices for the uncoupled elements the two 
systems. Next, the two subsystems are coupled. For this purpose, the 


deflection-force vectors must computed indicated Eq. 32. 
1,2 
These are: 
Here make the following identifications 
q2 = l vA and 42 = Pg 
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STRUCTURAL VIBRATIONS 


equating these two vectors, are led constraint matrix connecting 
the two sets coordinates Since have eight coordinates 
and six equations constraint there degrees freedom for the entire 
system. regard and free coordinates, which the first sub- 
script denotes the position the coordinate the subsystem vector and the 
second subscript denotes the subsystem the vector related. ob- 
tain the following constraining equation 


The inertia and stiffness matrices for the uncoupled subsystems are next 
determined follows: 


where the submatrices have already been computed. obtain the generalized 
inertia and stiffness matrices for the complete system operating the 
above matrices with Thus, 


Numerical values [m] and are: 


_ 
~~ 


Thus, eigenvalue problem formulated terms the vector 


j 
92,1 
93,1 
n 
K 
ge 


The corresponding normalized vectors are 


Using the constraint matrix may find the two vectors {q}, und 
for both modes. Then using matrices and for the two 
can obtain the values the coordinates for both subsystems 
and for both modes. Using these values and the functions the displace- 
ments may computed. The two natural mode shapes are shown Fig. 

may noted that, had dealt with this problem complete system 
rather than two subsystems, again would have had two degrees free- 
dom. this case would have had fifteen equations constraint the two 
junctions and seventeen coordinates, 

independent computation the first mode frequency the structure 
was made using the method transfer This frequency 


The value obtained herein exceeds this figure 3.6%. 


CONCLUSION 


method analysis has been presented for the determination natural 
modes and frequencies structural systems using the Rayleigh- 
Ritz technique. The method features procedure for obtaining mode shapes 
applicable the connected system synthesis admissible mode functions 
that are defined for the unconnected elements the system. This procedure 
eliminates the difficulty that usually attends the conventional Rayleigh analysis 
such systems. That is, the difficulty selecting appropriate mode shapes 
for the complete system. The method presented here permits the analyst 
partition the structure into subsystems that may treated independently 
the point final synthesis the generalized properties the complete sys- 
tem and solution the eigenvalue preblem. 

Two problem solutions have been presented and the results compared with 
other methods analysis undisputed accuracy. The results are considered 
quite acceptable for practical engineering analysis. 


APPENDIX. NOTATION 


Transformation matrix 
Transpose transformation matrix 
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STRUCTURAL VIBRATIONS 
Modulus rigidity tension-compression 
Modulus rigidity torsion 
Bending moments 
Mass element per unit length 
Axial force 
The j-th coordinate 


Kinetic energy the system 

Time 

Strain energy the system 
longitudinal displacement 

Shear 

Transverse displacement 
Independent variable measured along the length the member 
Torsional displacement 

Section radius gyration bending 
The j-th mode function 

Rotation due bending 


Angular frequency the vibration 
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IRROTATIONAL MOTION TWO FLUID STRATA TOWARDS LINE SINK 


SYNOPSIS 


The flow two fluids different density towards line sink located the 
bottom corner rectangular configuration with the upstream end extending 
infinity, was studied relaxation techniques. Relationships between the 
Froude numbers the two layers were determined and the critical condition, 
when the lighter fluid begins ceases flow, was estimated. 


INTRODUCTION 


Although much has been published recently problems involving stratified 
flow, few papers have dealt with discharge through outlets, orifices, towards 
sinks. Yih has made2 important contributions this respect, particu- 
larly for continuous density gradients. Bata thermal-gradient 
problem involving the flow cooling water the inlet thermal plant. 
Harleman has done* experimental work flow under sluice gates. 
Craya two-layered system with sink vertical boundary. 


Note.—Discussion open until January 1961. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Engineering Mechanics Division, Proceedings the 
American Society Civil Engineers, Vol. 86, August, 1960, 

Assoc. Prof, Mech. Engrg., McMaster Univ., Hamilton, Ontario, Canada, 

the Flow Stratified Fluid,” Yih, Proceedings, Third Natl. Con- 
gress Applied 1958, pp. 857-861. 

ASCE, Vol. 83, No. June, 1957. 

“Submerged Sluice Control Stratified Flow,” F.Harleman, Gooch, 
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“Recherches Théoriques sur Couches Superposées Fluides 


a 
x 
| 
4 
« 
by 
> 


was this last study, verified experimentally Gariel,® that formed the in- 
spiration for the present work. 


STATEMENT THE PROBLEM 


Fig. line sink located the bottom corner rectangular box 
which two fluid strata different densities extend infinity the -direc- 
tion; and are the velocity vectors infinity and are assumed 
the x-direction constant magnitude intheir respective layers; 
and and are the depths each layer far upstream and are considered 
this problem equal magnitude. The mass densities the layers are 
given symbols and shown. Each fluid considered incom- 
pressible, inviscid, and homogeneous, and the flow assumed irrotation- 


Line sink 
FIG, 1.—DEFINITION SKETCH FOR PROBLEM 


al. necessary, therefore, allow slip take place the interface. The 
pressures are common each fluid every point along the interface. 

the fluids are rest, the interface will horizontal with equal 
the strength the sink increased from zero, fluid will start flow and 


accelerate until, some critical condition, the interface will drawn down 
the sink and the lighter fluid will begin take part the flow. Further in- 
creases the strengthof the sink will cause increases the magnitude the 
total discharge andalso changes the proportions the two fluids the total 
flow. apparent that two the factors paramount importance the de- 
termination the interfacial shape and the proportions flow will inertia 
and gravity. might expected that the Froude number therefore con- 
siderable value parameter and that large Froude numbers the inertia 
effect would predominant, while low Froude numbers the gravity effects 
would particularly important. 


Under the foregoing conditions, this problem has its objective the deter- 
mination the relationships between the discharge from each layer, the den- 
sities, and the interfacial shape. 


“Recherches Expérimentales sur Couches Fluides 
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FLUID STRATA 
ANALYSIS 


Because each fluid consideredas incompressible, inviscid, homogeneous, 
and irrotational motion, the Laplace equation must satisfied throughout. 
Consequently, there must stream function each fluid field that satisfies 
the Laplace equation every point that fluid and also fulfills the boundary 
conditions imposed the problem. such stream functions can found 
each fluid field, the problem determining the velocities every point, in- 
cluding those along the interface, relatively simple. The horizontal andver- 
tical boundaries representing the top, bottom and configuration shown 
Fig. are simple boundaries, but the other boundary, the interface, more 
complex and must determined application the Bernoulli equation along 
each side the interface. order derive the interfacial equation, one 
point taken infinity, and the other some arbitrary position the in- 
terface (Fig. 1). Then, with symbols referring the upper fluid given the 
subscript and those the lower fluid the subscript the equations become 


and 


Eqs. and the pressure, the specific weight, the vertical 
distance measured from the horizontal line representing the interface in- 
the and the magnitude the velocity vector. Up- 
the equation the interface becomes 

Upon consideration Eq. boundary condition for potential-flow 
problem, the possibility rigorous mathematical determination the stream 
function each potential field becomes remote, for Eq. mixed and non- 
linear. 

Limiting and Critical Conditions.—Two extreme conditions may studied, 
critical one when the lighter fluid begins ceases flow, 
and limiting condition high discharges when gravity effects become neg- 
ligible. 

The first these conditions was studied extensively for his dif- 
ferent boundary conditions. made use the expression 
Richardson and subsequently Yih: 


gravitational acceleration, primes denote differentiation with respect and 


“Stationary Waves Water,” Richardson, Philosophical Magazine, S6, Vol. 
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G(w) function such that G1/3, and are real and fi- 
nite the interface. Craya was fortunate finding function G(w) for use 
inthis equation that satisfied approximately his boundary conditions, and through 
series transformations was able obtain numerical value for his 
limiting Froude number. the present problem the boundary conditions are 
much more complicated and corresponding value G(w) apparent. Be- 
cause the difficulties inherent this approach, mathematical analysis the 
critical condition this method was not pursued further, the hope that this 
relationship might determined with sufficient accuracy other methods. 

The limiting condition high discharges may determined rigorously 
assumed that the fluids started from rest that there vorticity 
the interface. Although this assumption was not stipulated the statement 
further analysis. Yih proved2 that under conditions high discharge 
stratified fluidin the absence vorticity andwhen gravity effects may neg- 
lected, every steady irrotational flow corresponds stratified flow the ve- 
locity whichis that forthe irrotational flow divided the square root the 
density. The streamline pattern for homogeneous fluid must have strati- 
fied flow counterpart with the interface along streamline the central one 
for the boundary conditions this problem. Since the velocity point 
the interface the homogeneous pattern must single-valued, the strati- 
fied case the following relationship must result: 


That similar relationship exists close tothe sink may shown when Eq. 
rearranged 


Multiplying through gives 


V2A V2A 

and point taken close the sink where becomes very large compared 

the other terms the equation except the limit the equation 


results, which not based any the assumptions necessary for the devel- 
opment Eq. 

the analysis close the sink continued little further relationship 
between the interfacial angle, the velocities, and the densities may derived. 
Fig. the angle that the tangent the interface the sink makes with the 
vertical and with the horizontal, arc small radius drawn, 
the discharge through the area and that through the area 
But the discharges may also written and and the 
ratio becomes since equal When the two expressions for 


the discharge ratio are equated, and Eq. substituted, the interesting relation- 
ship 
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results. 

apparent that very high discharges, Eq. accepted, the inter- 
face must make angle 45° with the axes close the sink. This fact also 
follows from the symmetry the homogeneous fluid streamline pattern 
which this limiting case stratified flow was assumed equivalent. addition, 
for any given interfacial angle the sink, the velocity and discharge ratios 
infinity are known the densities the two fluids are given. 

Dimensional functional relationship assumed follows: 


where and arethe coordinates the interface, according the Buckingham 
different functional relationship involving five dimensionless terms 
may writter 


only three fundamental units are included the relationship. When 
and are combined turn with the other variables, five dimensionless terms 


FIG, 2,—INTERFACIAL FIG, 3.—DEFINITION SKETCH FOR SIMILARITY 
ANGLE CONSIDERATIONS 
THE SINK 


the equally correct relationship 


Vom’ Ay,’ hi’ hi} 4 
results since The second two terms are Froude numbers, different 
from the customary Froude numbers because the inclusion the 
terms the denominator. Such Froude numbers are used Keulegan, 
and From this analysis, comparison Froude numbers 
the two layers for different velocity ratios andinterfaces appears useful. 
should remembered, however, that many other similar and possibly equal- 
useful relationships can found from Eq. 10. 
Similarity Considerations.—From similarity considerations two very im- 
portant results may obtained. Most important, may shown that the in- 
“Sixth Progress Report Model Laws for Density Currents: Effectiveness 


Salt Barriers Rivers,” Keulegan, Natl. Bur. Standards, No. 1700, June, 
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terfacial shape and the streamline patterns are dependent the Froude num- 
bers the two layers and not the density velocity alone. addition, the 
method presenting the results the analysis previously suggested dimen- 
sional considerations verified. 

The twodiagrams Fig. represent situations geometric and kinematic 
similarity. Therefore, the interfacial shapes and streamline patterns are sim- 
ilar each. Eq. written between corresponding interfacial points each 
diagram becomes 


and 


Eq. divided Eq. 14, the following equation results: 


Because the assumption kinematic and geometric similarity requires sim- 
ilar interfacial shapes and similar streamline patterns, and because and 
are corresponding points the interfaces, may replaced the 
length scale The right-hand side Eq. 15, which will called 
simplicity, becomes product three scale factors, that length, gravity, 
and density difference. cross-multiplication, Eq. becomes 


The corresponding potential fields Fig. must differ only constant fac- 
tor because the similar boundary conditions and the geometric and kinematic 
similarity considerations discussed previously. Therefore, from Eq. 16, the 
following relationships must hold: 


and 


From these equations, between corresponding points the similar flow pat- 
terns the following general relationship results: 


When the expression for ais substituted, Eq. becomes 
(22) 


rearrangement terms, the necessity for the equality corresponding 
Froude numbers becomes apparent: 


i 
a 
5 
ker 
ie 
| 
% 
tea 


EM4 FLUID STRATA 


The most convenient point the flow patterns choose reference 
infinity, and the similarity criterion, Eq. 23, becomes 


and ' 
g ho 42 g' hy 2 


single interfacial shape and corresponding flow pattern must result sin- 
gle ratio provided, course, that the requirements geometric sim- 
ilarity and potential flow both fluids are maintained. consequence, 
the method comparing the two Froude numbers and for various inter- 
facial shapes again suggested. interesting notice that Eq. which 
relates the interfacial angle the sink the velocity and density ratio, indi- 
cates the same conclusion, for the ratio the Froude numbers written 


for given system, the quantities Ap, and cancel out, and the ratio be- 
comes 


Systems Three More Layers.—The similarity analysis given for the 
two-layered system may extended systems three more layers, but 
the free choice density ratios found for the two-layered system not pos- 
sible. the case two interfaces, equivalent equations corresponding 
Eq. and Eq. applied each interface, turn, result two expressions 
for the velocity the intermediate layer. symbols referring the three 
layers are given subscripts 2,and where subscript applies the inter- 
mediate layer, the equations equivalent Eq. are 


and 


Since the corresponding velocities identical for any situation 
the same layer, the following relationship between the densities must hold: 


For more than three layers, additional relationships similar Eq. result. 
If, then, multilayered system some the densities are given, the others 
must satisfy relationships between the densities, such Eq. 29, order that 
certain potential-flow pattern results. The case the two-layered system 
seen unique, that for simple flow densities may 


chosen random, provided only that the Froude numbers the two layers 
kept constant. 
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SOLUTION THE PROBLEM 


The preceeding analysis indicated the desirability calculating and plotting 
the Froude numbers inthe upper and lower layers solution the problem. 
This solution requires the determination velocities and 
number interfacial shapes and potential field patterns. The difficulty ob- 
taining rigorous mathematical solution, because the complexity Eq. 
common-boundary equation, was pointed out the analysis. For this 
reason the methods approximate numerical analysis offer the most practi- 
cable way obtaining solution. The particular technique relaxation de- 


scribed10 McNown, Hsu, and Yih was followed the pres- 
ent work. 


FIG, 4.—DEFINITION SKETCH FOR RELAXATION PROCEDURE 


General Method.— The method used for the solution relaxation procedures 
necessitated, first, the assumption interfacial shape and, second, the es- 
timation the stream function throughout the two potential fields. re- 
laxation techniques the assumed values were altered until they were cor- 
rect every point for the assumed boundary conditions. Velocities were cal- 
culated number points along the interface and the interfacial equation, 
Eq. was checked determine the boundary condition was satisfied. 
the use constants the equation was possible calculate the two veloc- 
ities every point, and thus the Froude number, once the correct interfacial 
shape was obtained. 

Details the Relaxation Technique.—In Fig. the assumed values the 
stream function andvelocity are given symbols and respectively, dis- 
tinguish them from thetrue values The numerical values chosen were 
shown for two the three cases attempted. For the third case the value 
lower boundary was 19. The grid pattern used was one 0.1 unit 
squares progressively reduced 0.05 units inthe vicinity the interface, and 
finally 0.025 units close the sink. This pattern shown Fig. 


“Applications the Relaxation Technique Fluid Mechanics,” McNown, 
Hsu, and Yih, Transactions, ASCE, Vol, 120, 
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FIG, 6.—COMPARISON INTERFACIAL PROFILES 
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Each case attempted necessitated the assumption interfacial curve 
drawn that the first and second derivatives when plotted also gave smooth 
curves. The value each grid intersection was estimated and residues 
calculated using the standard For irregular star patterns near 
the interface, the residue was Taylor series expansion 
When the residues were relaxed very small value each intersection, the 
assumed interfacial shape was checked the calculation velocities each 
side the interface number points spaced along the interfacial curve, 
anda calculation todetermine these velocities satisfied Eq. the interfacial 
boundary condition. Two factors were introduced into Eq. 


and 
give 


pair points the interface was selected and values and Vo, 


found from the were substituted into Eq. along with cor- 
responding values Thedensity values used were 1.938 and 1.996. 
The resulting pair equations was solved simultaneously for and cor- 
rectly assumed interfacial curve wouldgive constant values Kand regard- 
less the pair points selected for the calculation. Kand were not found 
reasonably constant calculation involving successive pairs points, was 
necessary interfacial curve and repeat the relaxation procedures. 

was found that and were very sensitive interfacial alterations and 
was necessary make many trials, reduce the grid size and relax the res- 
idues little parts 10,000 before reasonably constant values were 
obtained. 

Three interfacial curves were determinedin this manner. These are shown 
Fig. along with the limiting profile homogeneous fluid flow. 

two quantities 


and 


are introduced into Eq. and are given their numerical values 
20, there results 


Eq. canbe solved for and between successive points without the neces- 
sity introducing numerical values and The Froude numbers can 
calculated directly from and using Eqs. and 


\ 
| 
4 
i 
; 
} 
4 
| 
{ 


come 


These expressions verify the result the similarity considerations devel- 
oped previously, which indicated that the interfacial shape was afunction the 
Froude number ratio, because their use the Froude numbers can calcu- 
lated directly from the results the relaxation procedures without the use 
numerical values for the densities. 

Results.—The final interfacial location and W-pattern for case are shown 
Fig. with streamlines indicated dashed lines. The three interfacial 
shapes are plotted Fig. along with the calculated central streamline for 
potential flow homogeneous fluid, and table coordinates included. 
Table shows the results the solution simultaneous equations for and 
for pairs points along each interface. The values and are averaged 
and the maximum deviation from the average values indicated percent. 
the upstream flat portions the interfaces, numercial difficulties resulted 
from the small size the velocity difference terms Eq. 31, and several 
points this region were omitted from the average values Kand Asa 
check the physical meaning the errors shown and the pressure 
difference between pairs points along each side the interface was calcu- 
lated shown Table The points the flat portion the interface were 
included inthis check, but the maximum errors found than 1%. Table 
2gives values Fy, Fo, and for each interfacial shape. 
The values Land Mwere calculated from the average values Kand Jgiven 
Table the use The curve shows the relationship be- 
tween the two Froude numbers. Because the intercept the 
and become zero, not apparent from this curve, another relationship, 
that the ratio Froude numbers and the sum the square roots the 
Froude numbers, shown Fig. This relationship was plotted order 
obtain abscissa proportional the total discharge but dimensionless 
terms. Fig. curve plotted against comparing the angles calcu- 


lated from points Fig. the use Eq. and those measured from the 
interfacial curves, Fig. 


DISCUSSION RESULTS 


comparison the curves and enables estimate tobe made 

the critical Froude number the lower layer which the fluid the upper 
layer ceases flow. From Fig. the value this point seen 
with reference Fig. appears reasonable. Craya value 
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TABLE 1.—VALUES AND 


Interface Interface No, 


0.0267 0.0336 
0.0259 0.0329 
0.0250 0.0337 
0.0257 0.0332 
0.0246 0.0328 
0.0271 0.0334 
0.0262 0.0334 
0.0261 0.0336 
0.0239 0.0329 
0.0276 0.0337 
0.0 


0.0325 
Average 
without 


Max. 
Error 
without 


Average 
without 
AorB 


TABLE 2,—FROUDE NUMBERS 


0.004343 
0.03637 


0.04683 
0.03791 


| 


Error Error 

face 
0.1080 1.778 0.0607 1.663 
0.5571 0.9432 0.5906 1.718 
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TABLE 3,—INTERFACIAL COORDINATES 
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TABLE DIFFERENCE ALONG INTERFACE 


Interface Difference 
Fluid 
16,272 


16.573 


15,341 
2.581 
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0.364 0.434 0.4485 
0.294 0.383 0.4038 
0.247 0.341 
0.211 0.305 0.3280 
0.1581 0.2441 0.2660 
0.13785 0.2184 
0.1205 0.1951 0.2154 
0.08105 0.1380 0.1561 
0.07105 0.1230 0.1400 
0.1009 
0.0815 
0.0735 
0.0665 
0.0604 
0.0551 
0.0505 
0.0465 
0.0399 
0.0371 
0.0345 
0.0320 
Error Error 
0.011 
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FIG, 7.—FROUDE NUMBER RELATIONSHIP 
FIG, NUMBER FIG, 9.—INTERFACIAL ANGLE VERSUS 
RELATIONSHIP FROUDE NUMBER 
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for this Froude number with his quite different boundary conditions. 
The significance ofthis for any pair fluids any geometrically 
similar potential-flow situation, the velocity infinity the lower fluid 
which the upper fluid either starts ceases participate the flow can 
calculated. 

The curves drawn verify the suggested result that large discharges 
the inertial effects predominate and the Froude numbers approach each other. 

Since, Fig. the point which both fluids are rest the origin, 
apparent that the curve drawn must start this point when the fluid be- 
gins accelerate, and continue along the x-axis until the critical Froude num- 
ber reached. Here the upper fluid begins flow, the total discharge in- 
creases slowly (Fig. 8), but the discharge the heavier fluid decreases 
that the lighter fluid increases. point reached, however, when 
changed and both discharges increase together the Froude numbers ap- 
proach equality. would expected that considerable instability would result 
the neighborhood the point minimum curve Fig. ata 
Froude number the lower fluid about 0.9. about this same Froude 
number minimum point occurs the curve Fig. Below this point the 
angle approaches its limit 45° very abruptly increases; above this 
inflection point the angle approaches the critical value 90° more gradually 
increases its critical value 2.76. The values scaled from the 
interfacial profiles, Fig. agree closely with the curve Fig. which was 
plotted from calculated values Froude number, thereby verifying the rela- 
tionship between the Froude numbers and the angles expressed Eq. 26. 

comparison the three interfacial profiles plotted Fig. and the 
central streamline calculated for homogeneous fluid suggests that the inter- 
face successively takes flatter shapes the proportion the lighter fluid 
the flow increases until, condition equivalent the inflection points 
Figs. and the trend reverses and the limiting profile approached. 


The irrotational motion two towards sink located within 
the boundary configurations chosen this problem has been shown gov- 
erned the relationship between the Froude numbers each layer ex- 
pressed the curve Fig. This relationship indicates numerical value 
for the critical Froude number the lower stratum which the upper fluid 
either begins ceases flow, and asymptote representing the equality 
Froude numbers high rates flow when inertial effects predominate. 
interesting feature this curve inflection point which, for increasing 
total flow, the rate flow the lower fluid begins increase along with that 
the upper rather than decrease. 

The shape the interface and the streamline pattern have been shown 
dependent the Froude number the upper layer, since the angle the 
interface the sink shown depend the Froude number ratio, and from 
Fig. there only one ratio possible for each value other than zero the 
Froude number the upper stratum. 


é 
ot 
x 4 ip 
= 
4. 
| 
q 
| 
Ye 


ACKNOWLEDGMENTS 


The writer gratefully acknowledges the direction and help Hunter Rouse, 
ASCE, the State University Iowa, who directed his dissertation, 
Yih, the University Michigan, under whose direction the work was started, 


and Macagno and Landweber Iowa who were generous with their 
knowledge and advice. 


| 
ik 
| 

| 

| 

ve i 
| 

1 

~ 

Br 
i 
f 


August, 1960 


Journal the 
ENGINEERING MECHANICS DIVISION 


Proceedings the American Society Civil Engineers 


UNDERGROUND STRUCTURES SUBJECT AIR OVERPRESSURE 


Eben ASCE 


SYNOPSIS 


analytical procedure presented for determining the damage un- 
derground structures induced pseudo-steady state air overpressure. 
introducing appropriate assumptions from the fields soil mechanics and 
structural dynamics, possible derive analytical expressions that re- 
late the parameters the structure, the soil, and the loading they affect 
the failure the structure. The assumptions are stated and the derivation 
the analytical procedure presented. separate section devoted the 
application the procedure for the cases analysis the vulnerability 


underground structure and the design underground structure for any de- 
sired level protection. 


INTRODUCTION 


This analytical procedure only concerned with damage induced the 
air overpressure acting the ground surface above the buried 
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structure. “damage” meant actual structural failure. Stress wave dam- 
age not considered because felt that the locations interest this 
agent, general, causes only minor structural damage. 

With sufficiently high overpressure the ground surface, stresses great 
enough fail the structure may induced underground structures through 
the surrounding soil. The structures will tend deform under the action 
these stresses, but failure cannot occur unless large masses soil can move 
with the structure order maintain the stresses. will shown that for 
given structure and given soil the overpressure required cause failure 
increases rapidly with the depth burial. 

order get movement the soil necessary that the shearing 
stresses inthe soil exceed its ultimate shearing resistance. assumedthat 
this shearing resistance afunction the normal stress onthe shear planes, 
the angle internal friction and the cohesion the soil. The shearing re- 
sistance inthe region soil structure limited the strength 
the structure. When the structure begins yield, the soil shearing re- 
sistance mobilized and increases until the ultimate shearing resistance 

further increase overpressure will extend the regionin shear 
resistance exceeded until this region reaches the ground surface. order 
tohave any movement the soil, the ultimate shearing re- 
sistance exceeded along some surface extending from the structure the 
ground surface. 

The general approach the problem select planes within the soil 
which the shear stress equal the ultimate shear resistance and then 
write the equation motion the mass soil which bounded these 
planes, the structure, and the ground surface. This particular volume soil 
will referred asthe “failure mass” and the shear planes “slip planes.” 
The solution this equation motion results expression relating the 
parameters the structure, the soil and the blast loading they effect the 
failure the structure. 

Notation.—The letter symbols used this paper are defined where they 
first appear, the text illustration, and are assembled for convenience 
reference the Appendix. 


ASSUMPTIONS 


Blast Loading.—The uniform air overpressure applied the ground sur- 
faces, that is, the pseudo-steady state overpressure (t), varies with time. 
This variation represented initially peaked, linearly decreasing over- 
pressure duration tg, given by4 


which pg(0) the initial peak overpressure for duration tg, time. 
Eq. shown graphically Fig. The loading parameters are thus (0) 


Engineering Approach Blast-Resistant Design,” Newmark, Trans- 
actions, ASCE, 121, 1956. 
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OVERPRESSURE 


The actual variation overpressure with time the 


and shown graphically Fig. The relationship between the actual and 
linearized overpressures normally maintains the same peak overpressure and 
total impulse. Thus found function to, the duration actual over- 
pressure, equating the area under the curve Fig. with that Fig. 


the actual wave form differs from that given Eq. then Eq. should 
modified accordingly. The parameters (0) and are related the weapon 
yield, height burst, and ground 

Soil Behavior.—Several fundamental assumptions are made regarding the 
behavior the soil under the type loading previously considered. These are 
based conventional static and are follows: 


The relationship betweenthe vertical and horizontal components stress 
point the soil may expressed 


which the coefficient earth pressures, the horizontal, and 
the vertical soil stress. After the soil has been deposited, whether naturally 
artificially, becomes K,, the coefficient earth pressure rest. the 
active Rankine state plastic equilibrium exists then the coefficient becomes 
Ka; the passive then the coefficient becomes K,. 

The the slip planes the soil assumed fol- 
low Coulomb’s equation 


which the ultimate shearing resistance the soil, refers the soil 
cohesion and represents the soil’s angle internal friction. 

For the active Rankine state plastic equilibrium the slip planes make 
angle with the horizontal direction. 

The compressibility the soil negligible. 

The failure mass can treated rigid body. 

Because the level the water table sufficiently far below the ground 


surface, hydrostatic pressures have considered. may reasonably 


Simple Method for Evaluating Blast Effects Buildings,” Review Edition, Ar- 
mour Research Foundation, July, 1954, 
“Effects Nuclear Weapons,” Printing Office, 1957. 
“Soil Mecahnics Engineering Practice,” Terzaghi and Peck Wiley, 
New York, Y., 1956. 
“Theoretical Soil Mechanics,” Terzaghi, Wiley, New York, Y., 1956. 
“Stress Distribution Dry and Saturated Sand Above Yielding Trap Door,” 


Terzaghi, Proceedings, Internatl. Conference Soil Mechanics, Cambridge, Mass., 
Vol, 1936, pp. 307-311, 
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ment being considered. the water table above the element may as- 
sumed that the element subjected tohydrostatic pressure, that is, the stress 
will equal the overpressure, pg(t), the soil all directions. 

Structure.—The mode failure for structural element assumed 
rigid body motion parallel its initial position. For shear failure this mode 
essentially correct while for bending failure best approximation. 
Two parameters define the structural elements. These are the ultimate re- 
sistance, and the maximum plastic displacement, 

The ultimate resistance structural elements, such walls and roofs, 
may determined standard analytical techniques. rigid plastic re- 
sistance assumed for these elements and, therefore, necessary con- 
sider only the maximum resistance under uniform load. This resistance, pm, 
expressed units pressure (pounds per square feet). 

suggested that set equal the maximum allowable displacement 
for shear mode failure and equal one-half the maximum displace- 
ment for bending failure. This assumption based the desire conserve 
the displaced volume. should pointed out that the maximum displacements 
being considered are those associated with complete failure, not merely initial 
signs failure such cracking. 


ANALYTICAL DEVELOPMENT 


Roof Failure.—The dashed lines Fig. indicate the selected slip planes. 
The failure mass this case that enclosed the dashed lines, the roof 
and the ground surface. free body diagram for horizontal slice the 


failure mass shown Fig. and the forces shown are evaluated fol- 
lows: 


the force acting the top the slice, isthe vertical stress 
soil, the length the wall roof, and structural dimension. 


the force acting the bottom the slice and the variable soil 
depth. 


the weight the slice the failure mass and the unit weight 
the soil. 


The forces and are the shearing forces the vertical faces 
the slice failure mass. assumed that the normal stress these faces 
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and 


The forces and are the normal forces the vertical faces and 
their effect resisting motion included Eqs. and 12. 


The solution the equation motion this slice the “failure mass” 
gives the relations between the various parameters that determine the failure 


the roof. Using represent the vertical motion the slice this equa- 
tion motion 


32.2 


Q+W-P-2Vp - 


placement the failure mass. expressions into Eq. yields 


which the structural parameter. Since assumed that the failure 
mass undergoes rigid body motionthe parameter independent There- 
fore Eq. becomes ordinary differential equation and Solving for 
(z,t) and imposing the boundary conditions 


and rearranging get 
which 
and 
x(t) = 
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and which isthe soil depth tothe top the adimensionless 


parameter, and the acceleration the failure mass the soil. Eq. 
the form 


which are constants. With the initial conditions (0) (0) 
where the velocity the failure mass the soil. The closed form so- 


where the pseudo-period, T', defined 


F=144p ... 


Hence substituting Eqs. into Eqs. obtain the solution Eq. 
15. When the duration infinite the overpressure constant with time and 
the peak overpressure required cause failure the roof given 


For convenience the closed form solution given Eqs. has been shown 


the initial overpressure for infinite duration. Using Eq. for and Eq. 
for then the peak overpressure the overpressure-time 


tory required fail the structurefor any duration may determined from 
Fig. 


comparison Eq. with Eq. shows that 
¥, | H r 
and 
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WITH SOIL AND STRUCTURAL PARAMETERS 
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OVERPRESSURE 


order demonstrate the effects each the parameters and 
Figs. and have been prepared. Representative values each the 
parameters were selected. These values are follows: 2000 psf; tan 

The curves Figs. and were obtained varying the parameters one 
time while holding the others fixed the selected representative values 
previously listed. 

Wall Failure.—The dashed lines Fig. indicate the selected slip planes. 
The failure mass this case that enclosed the dashed lines, the wall, 
and the ground surface. deriving the equation motion the failure mass 
divided into two parts, the upper mass and the lower mass 

The equation motion exactlythe same asthat the failure mass 
associated with roof failure, Eq. 15, with replaced tan (45 /2) and 
pt, the normal stress between failure masses and 

The equation motion determined with the following assumptions 
Fig. 9). 


The mass may treated rigid body. 

The direction motion parallel the inclined face Mg. 

The vertical component motion equal that 

There are shear forces induced the common surface between 
and 

The vertical stress and hence the horizontal stress constant through- 
out the failure mass Mg. The value the vertical stress assumed 
equal the stress acting the top surface the failure mass. 

The normal force given 


which the height the wall. 
The normal force given 


The induced shear forve parallels the direction motion and giv- 


The resistance the wall p,, acts normal the wall, which means that 
there are shear forces along the surface the wall. The force the wall 


10. The force related the force Coulomb’s equation (Eg. 7). 
The magnitude these forces are selected satisfy the equilibrium 
requirements 
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OVERPRESSURE 
the basis these assumptions equation motionfor and may 


writtenthat relatesthe resistance the wall the overpressure 
This equation 


KH 


2Km tan 


tan 
where 
LH, tan (45 /2) 
tan 
and 
KH, tan 


and which dimensionless parameter. 


For triangular overpressure-time history this equation can written 
the form given Eq. where 


208 (45-922), in (45 - 9 2) 


tan 
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which L', and are previously defined 


When the duration infinite the overpressure constant with time and the 
peak overpressure required fail the wall given by: 


p_(o) 
for the roof analysis, Fig. presents the ratio function the 


duration, tg, and pseudo-period, T', which defined Eq. 19, using 
Eqs. and for and 


order demonstrate the effects each the parameters pg(o) and 
Figs. and have been prepared. The following representative values 


MOTION 
FIG. 9.—FREE-BODY DIAGRAM LOWER FAILURE MASS, 


each the parameters were selected: 2000 psf; 203 K=1/2;y 
The curves Figs. and were obtained varying the parameters one 


time while holding the others fixed the selected representative values 
listed above. 


APPLICATIONS 


There are two major applications the analysis presented inthe preceding 
sections. One the analysis given structure, that the determination 
the overpressure-time history required tofail particular underground struc- 
ture given location. should noted that the overpressure-time history 
completely described the peak overpressure and the duration. The other 
application design underground structure withstand prescribed 
overpressure-time history. 

Analysis.—Where structure has one element tobe analyzed, for 
roof and two dissimilar walls, each element should analyzed in- 
dependently. this casethe lowest peak overpressure for given duration, 
obtained considering failure each the structural elements, considered 
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the peak overpressure that causes failure the structure. The following 
steps outline the procedure followed order determine the peak 
overpressure-duration relationship required fail particular structure. 


Select the soil properties (from Table other information not 
available). 

Determine and for each element being considered. 

Compute (Eq. Eq. 33) 

For the element with the lowest value compute (Eq. 19). 

For any to, may computed from Eq. With and from 
steps and the overpressure needed fail the structure for any duration, 
may found from Fig. 


example, consider cubic structure 10ft side, buried that the 
roof underground. The soil surrounding the structure soft moist 


TABLE 1,—APPROXIMATE VALUES SOIL 


Type Soil 


Angle internal Soil cohesion, 


degrees 


lbs 
per 


Coefficient 
earth pressure 
K 


Soft moist clay 
Medium moist 


o 


clay 0.7 
Stiff moist clay 0.7 
Hard clay 0.7 
Sandy loam 0.6 
Dry loose sand 0.5 
Dry dense sand 0.6 


Soil densities relatively narrow range. For this reason, and since 
this parameter has only limited influence the results, single value 100 per 
was 

will noted that the values for given this table are less than unity, Based 
Terzaghi’s paper, “Stress Distribution Dry and Saturated Sand Above Yield- 
ing Trap Door,” one might expect the value exceed unity for dry sand. For the 


purposes this table, values which would conservative for design have been select- 
ed, 


clay. The following steps, which parallel the previous sequence, show the 
failure analysis for this example. 


According Table the values soil parameters for soft clay are 
follows: 10°; 400 psf; 100 pcf; and 0.6. 

The ultimate resistance the roof computed from standard tech- 
niques estimated from known design information. For this example as- 
sume that for the roof 1400 psf, For the wall assume the 
following values: p,, 800 psf; and ft. 

Usingthe values determined above one obtains, from Eq. 23, for the roof 
psi and from 33, for the wall, psi. Comparing the 
value the roof with the wall, can seen that for this 
example the lower peak overpressure infinite duration that the roof. 

From 19, using the values determined previously for the roof, one 
obtains 0.183 sec. 
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Assume the duration equals sec. Then from Eq. equals 1.47 
sec and from Fig. pg(o) psi. This step can repeated often 
possible obtain other combinations duration and overpressure. 


Design.—Application the analytical procedure design probably best 
treated iterative process. Thefirst trial design maybe made 
tional methods based cost, use, etc. This structure may then subjected 
the method analysis previously outlined. the analysis shows that the 
structure does not have the desired resistance the design can modified and 
the analysis repeated until obtained that satisfiedthe requirements. 
The preceding example problem canthus serve example Based 
assumed values for and relationship between overpressure and 
duration can established. this overpressure satisfies the de- 
sign requirements then the assumed values and are satisfactory. 

One useful simplification results the loading considered in- 
finite duration. Under this condition inertial effects, given are 
eliminated, the arithmetic simplified, and the results are conservative. 
other words, structure designed for psi this method could expected 
withstand well over psi for any real condition. 

Extensions the theory have been made non-flat roof wall surfaces, 
such arches, domes and cylinder. The approach used for these shapes was 
define equivalent flat surface and proceed exactly indicated this 
paper. Such applications can considered, best, suitable for interim 
use until further research provides better solution. 


CONCLUSIONS 


The methods outlined this paper are based the assumption that the 
soil parameters are known and that they remain constant during the entire 
loading process. account taken non-uniformity soil possible 
variations soil shear parameters during the yielding process. quite 
possible that stress distributions the soil may vary yielding 
takes place. cases where any these effects are very pronounced the an- 
alytical methods given might not expected give sensible results. 

The analytical approaches this paper are dynamic modifications ac- 
cepted static Even for static conditions these approaches repre- 
sent simplifications introduced order make possible arrive so- 
lution. Certainly questions arise regarding the location failure planes, be- 
havior the failing structural elements, and behavior the soil. 

The work reported represents the only experimental study 
known the authors for the static conditions. These experiments form the 
basis for Terzaghi’s theoretical approach, which was modified for use this 
paper. Based this paper one might also conclude that the value for dry 
sand should excess unity. Application this value for questionable and 
the values recommended inthis paper (Table were selected asto con- 
servative for design. The uncertainty these values should noted however 
since the value has important influence the results. 

Therefore, the method this paper should only considered afirst ap- 
proximation for the design analysis underground structures subjected 
air overpressure. Future work, both theoretical and experimental can ex- 
pected provide improvements the analysis presented herein. 
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The following list symbols and units used the analysis. 


structural dimension, feet 

soil cohesion, pounds per square feet 

acceleration gravity, feet per square second 
soil depth top structure, feet 

coefficient earth pressure 

coefficient active earth pressure 

coefficient earth pressure rest 

coefficient passive earth pressure 

structural parameter, feet 

dimensionless parameter 

soil failure masses 


force bottom horizontal soil slice, pounds 

horizontal soil stress, pounds per square foot 

resistance structural element, pounds per square foot 

normal stress masses and Mg, inpounds 
per square foot 

time varying overpressure, pounds per square inch 

initial peak overpressure for duration, tg, pounds per 


square inch 


initial overpressure for duration, pounds per 
square inch 


2 


vertical soil stress, pounds per square foot 
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force top horizontal soil slice, pounds 
wall resistance, pounds 


dimensionless parameter 


ultimate shearing resistance soil, pounds per square 
foot 


time, seconds 

duration actual overpressure, seconds 
duration linearized overpressure, seconds 
seconds 


V2, Shearing forces the failure masses, pounds 


weight horizontal soil slice, pounds 


displacement failure mass, feet 
maximum plastic displacement structural element, feet 
velocity failure mass soil, feet per second 


acceleration failure mass soil, feet per square 
second 


variable soil depth, feet 
unit weight soil, pounds per cubic foot 
vertical stress soil, pounds per square foot 


angle internal friction soil, degrees 
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GEOMETRY MOIRE FRINGES STRAIN ANALYSIS 


Stanley August Durelli,2 and Cesar 


SYNOPSIS 


This paper shows how moire fringes can used the two-dimensional 
analysis strains. The fundamental equations the moiré method are de- 
rived and are presented the form graphs means which strains and 
rotations can obtained from simple measurements with minimum 
computation. 


INTRODUCTION 


The moire watered-silk effect optical phenomenon produced when 
two somewhat similar arrays dots lines are superimposed, resulting 
the formation alternating light and dark fringes. effect common 
everyday life, being seen mismatched window screens, the 
blind” defect television, and forth. 


Note.—Discussion open until January 1961. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Engineering Mechanics Division, Proceedings the 
American Society Civil Engineers, Vol. 86, No. August, 1960, 

Associate Engineer, Mechanics Research Div., Armour Research Foundation, 
Chicago, Illinois, 

Professor Civil Engineering, Illinois Institute Technology and Supervisor, 
Stress Analysis, Armour Research Foundation, Chicago, 
Associate Research Engineer, Institute Technology, Chicago, Illinois. 
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Because the fringes are the result displacement the elements the 
two arrays with respect each other, the effect obviously tool 
which strain measurements can made. If, for example, one begins with two 
identical sets parallel lines, fringes are produced either rota- 
tion (Fig. 1(a)) elongation (Fig. 1(b)) one set with respect the other. 


59 lines 


(a) ROTATION ALONE 


(b) DIFFERENCE 
PITCH ALONE 


COMBINATION 
ROTATION AND 
DIFFERENCE PITCH 


51 


per inch 


Simple equations govern the interpretation pure rotations and elonga- 
tions. Within the limits small deformations, these effects can considered 
linear and capable being superimposed. number papers have 
been published under these assumptions, notably Weller and 
and 


“Utilization des Reseau pour des Deformations,” Dantu, Laboratoire 
Central des Ponts Chaussees, Paris, Publication 57-6, 1957. 
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Actually, the effect not limited any way separate rotations 
and elongations nor infinitesimal deformations. Neither necessary that 
the arrays originally identical either spacing orientation. 

this paper only sets straight, parallel, nondiffracting lines will 
considered. However, moiré fringes that can analyzed could obtained 
from circular, radial, other nonparallel arrays, and diffraction gratings 
could used under suitable conditions. 

Fig. 1(c) shows the effect combined rotation and differences pitch 
straight, parallel sets lines representing homogeneous field for which 
equations will derived. The relationships obtained can also applied 
sufficiently small elements nonhomogeneous field such would 
presented general two-dimensional strain problem. 

The moire method can certainly used for the analysis large strains 
and rotations, but such application beyond the scope this paper. How- 


FIG, MOIRE FRINGES 


ever, the equations will derived the most general terms before being re- 
lated such specific quantities strain. Approximations for the case in- 
finitesimal deformations will also developed. 


BASIC PROPERTIES AND DEFINITIONS 


the subsequent derivations, fixed array straight, parallel lines called 
the “master grid” used reference both for analysis the properties 
the moire and for establishment coordinate directions. The center-to- 
center distance between the master grid lines defined the master pitch, 
and the directions perpendicular and parallel these lines are designated 

The model also bears similar array lines which the undistorted 
state are not necessarily the same pitch the master grid. The pitch 
this grid any particular point any state distortion called the model 
pitch” and designated 
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The angle defined the acute angle any point and any given state 
deformation measured from the fixed master grid lines the model grid 
lines. The angle from the fixed master grid lines fringe point meas- 
ured the same direction designated and may either acute 
obtuse. Fig. shows how obtain the inclination the fringes and 
their distances functions the master grid pitch the model grid pitch 
and the rotation both grids point. 

Let assume either the case homogeneous deformation and rotation 
(which may accompanied translation) sufficiently small element 
nonhomogeneous field. Equations relating any desired sets parameters 
for the moire phenomenon are easily derived from the simple geometry 
fringe formation. Those for thenormal fringe spacing, and the fringe angle, 
are perhaps the most fundamental. However, the information available for 
the analysis will usually the form photographs which and and 
the coordinate fringe spacings and can measured, and the purpose 
the analysis will the determination and the desired points. Obvi- 
ously, the master grid pitch, and the direction its lines must known. 

Several the expressions developed herewith can obtained other 
ways, some using the property fringes being loci points 
equal displacement. The derivations that follow, however, are very simple 
and general. representing relationship between parameters 
may prove useful practice. 

Inclination the Fringes. 


and 
therefore 
Because 
then 
sin 
Distance between Fringes.—Also from Fig. and the foregoing 
(6) 
and 
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From Eq. 


and 


Rotation Terms the Inclination the Fringes and their Spacing.— 
From Eq. 


and 


represent graphically the relationship between and convenient 
introduce the following definitions: 


(12) 


and 


Values for the entire possible range (0° 180°) and for various 
values are plotted Fig. permitting easy determination the 
model rotation point without knowledge the model pitch. 


From Eq. 11: 
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and 


+ 
+ 
_ 
«J 


Determination and Strains.—To compute the val- 
the strain using the measured quantities and graphical construc- 
tion can utilized. For this construction convenient introduce the 
following definitions: 


and 


Values for all possible values and for various values 5/p have 
been plotted Fig. permitting easy determination the model pitch, pitch 
ratio, the quantity free from the influence rotation. should 
noted that (p/p') the “true” strain perpendicular the direction 
the master grid lines and were originally equal before distortion the 
model. 

Eq. converted the form: 


obtain the conventional” “nominal” strain perpendicular the direc- 
tion the master grid lines. 

Rotation Terms Distance between Fringes along the Coordi- 
nate Directions.—In many cases may preferred measure the fringe 
spacing along the coordinate directions may desirable have sec- 
ond method for checking the model rotations. From Fig. have, for the 
two cases (a) acute, and (b) obtuse, respectively: 
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n 
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wo 


and 


Substituting Eqs. into and extracting sin trigonometric trans- 
formation, obtain: 
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Similarly, from Fig. 


avoid having give negative sign 5,, measured distance, these 


equations can combined as: 


Yr? 


CASE (b), OBTUSE ANGLE 


FIG, 5.—SKETCHES DERIVE THE FRINGE SPACING 
THE COORDINATE DIRECTIONS 


Again substituting Eq. and extracting sin 
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Following the same method used previously can define and 
follows: 


and 
Then 
(28) 
and 


Values and are plotted Figs. and respectively. Their com- 
combined usable ranges cover all possible values 

Model Terms Distances between Fringes along the Co- 
ordinate Eq. into Eq. and setting relation- 
ship similar Eq. can define: 


o 
w 


and 


Similarly, substituting Eqs. into Eq. 17: 


and 
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(DEGREES) 


6.—ROTATION FUNCTION FRINGE ANGLE AND FRINGE 
SPACING THE DIRECTION THE MASTER GRID LINES 
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FOR 


FIG. FUNCTION FRINGE ANGLE AND FRINGE SPACING 


NORMAL THE DIRECTION THE MASTER GRID LINES 


Values and are plotted and respectively. Their com- 
bined usable ranges cover all possible values 

Simplified Equations for Small Deformations and Rotations.—It apparent 
from Fig. that sufficiently close 180° (depending also the 
value the value close -1.0 This means that rota- 
tions are very small many problems, sufficient accuracy may 
result from taking and reducing Eq. to: 


€true p’ 


simple transformation 


If, furthermore, sufficiently small comparison with unity (that 
strains are very small), Eqs. and both reduce to: 


the other hand, Fig. sufficiently close 90° (depending 


the value the value close 1.0, and sufficient accuracy may 
result from reduction Eq. 


Determination the Direction the Fringe Angle.—In order determine 
the direction which the fringe angle, measured the foregoing 
equations and graphs, small amount information needed addition 
photograph the fringe pattern. For example, the direction rotation may 
often deduced from the loading conditions could always found 
means coarse set lines dots the model addition the closely- 
spaced array which produces the Once the known, the 
direction determined its definition and all ambiguity disappears. 

Alternatively, all ambiguity removed the sign the quantity 1-(p/p') 
known, that is, whether greater less than unity. This can usually 


Then: 
and 
1 8 B (33) 
| 
: 1 + 
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FIG. 8.—MODEL PITCH FUNCTION FRINGE ANGLE AND FRINGE SPACING 
PERPENDICULAR THE DIRECTION THE MASTER GRID LINES 
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FIG, 9.—MODEL PITCH FUNCTION THE FRINGE ANGLE AND FRINGE 
SPACING THE DIRECTION THE MASTER GRID LINES 
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determined from the moire pattern the actual loaded model and master 
grid are available. Differentiating Eq. leads 


Because the denominator Eq. cannot negative, the sign the 


derivaties determined the sign [(p'/p) cos this expression 
negative, 


meaning that decreases for the master grid rotated 
given direction with respect the model and the fringes rotate the 
same direction, this condition satisfied and p'/p the grids were 
originally identical, this positively indicates tensile strain. 

For cos negative the situation not quite definite because 


can satisfied number ways when large angle. If, 
not near 90° and the fringes rotate opposite the rotation the master 
grid may ordinarily considered less than indicating compres- 
sion for originally identical grids. 

Other guides the sign are sometimes available directly 
from photograph. 90° any point, from Eq. 17. Once 
fringe known represent p/p' any point, the same fringe cannot 
represent p/p' without going beyond 90°, actually through 6/2. 

Fringes parallel the master grid lines represent the condition ro- 
tation (Eq. 11). Unless they are continuations from areas where the con- 
dition can evaluated, impossible distinguish from photograph 
alone whether p/p' greater than less than unity. 


USE THE BASIC RELATIONSHIPS TWO-DIMENSIONAL 
STRAIN ANALYSIS 


order determine strains and rotations associated with two coordinate 
directions actual model, will obviously necessary have models 
with grid lines perpendicular these directions. the model not identical 
when turned 90°, two models must equipped with grids, or, alternatively, 
two sets grid lines can placed the same model and analyzed separate- 
with one-way master grid. 

Individual measurements the distance between two fringes cannot 
very accurate, will apparent from examination Figs. and 11. 
Also, the point being examined may not exactly fringe centered be- 
tween two fringes. 


p 
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much more convenient and accurate measure across series 
fringes along line including the point interest and plot the accumulated 
fringe distances from any convenient point Fig. 12. smooth curve 
drawn through the plotted points, and the slope this curve any point gives 
fringes per inch, the reciprocal which the value inches per fringe 
the case may be. Fig. 12(b) because there rotation 
along axis symmetry, the process differentiation gives the true strain 
directly. 


FIG. 10.—MOIRE FRINGES DISK HYSOL 8705 UNDER VERTICAL 
DIAMETRAL COMPRESSION WITH HORIZONTAL GRIDS 
(GRIDS 300 LINES PER INCH DISK) 


will apparent from Figs. and 11, even simplified case such 
same point with anything approaching accuracy, and some points neither 
can measured. However, can always measured and values 1-(p/p') 
and from measurements can supplement check those from measure- 

The measurements will, general, not very accurate. However, 
since one can usually make choice from among the measurements 
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and any given point interest, the sensitivity errors measure- 


ments can kept fairly low level for most fringe configurations. 


FIG, 11.—MOIRE FRINGES DISK HYSOL 8705 UNDER VERTI- 
CAL DIAMETRAL COMPRESSION WITH VERTICAL GRIDS 
(GRIDS 300 LINES PER INCH 4-IN.-DIAMETER DISK) 
The values (p/p') and found with grids the two directions 
comprise all the information needed solve the general strain problem. 
addition the normal true strains: 
with grid lines the y-direction, and 
a 
“y true ~ fi oe eee we 


with grid lines the x-direction. 
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Their nominal (conventional) strain counterparts are 


and refers the corresponding 5,, 
have also the rotations 


are 


arc sin arc sin 
with grid lines the y-direction, and 
arc tan 


arc sin arc sin 


with grid lines the x-direction. 


general, there will rigid body rotation, and signs are given 
and (say for counterclockwise), the rigid rotation will be: 


and the shear strain will be: 


APPLICATION THE METHOD ACTUAL PROBLEM 


The moiré method has been applied the familiar case circular disk 
under diametral compression. order obtain satisfactory level re- 
sponse with the 300-lines per in. grid available, sheet was made urethane 
rubber (Hysol 8705) which the grid was photoprinted using the sensitized 


resist marketed Gaco. After printing, 4-in. disk was machined from the 
1/2-in.-thick sheet. 


nom 
with grid lines the y-direction, and 

(43a) 
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The disk was loaded between flat plates, first with grid lines perpendicular 
the loading directions (Fig. 10), then with the grid lines parallel the 
loading direction (Fig. 11). Photography was the double-exposure method, 
the unloaded and loaded conditions being successively exposed the same 
sheet Kodalith film with intermediate change except for application 
the load. (In this method fringe formation the film-and-lens system must 
capable clearly resolving the individual grid lines high contrast.) 


ALONG 
TANGENCY 
| 
ALONG HORIZONTAL 
DISTANCE CURVES 
/ 
0.0 
a / 
(b) FROM FIG. 
CURVES 
/ = 


ACCUMULATED FRINGE DISTANCE 


PLOT ACCUMULATED FRINGE DISTANCE FROM ANY 
ARBITRARY BEGINNING 


Fig. plot the measured positions successive fringes along two 
representative lines where could successfully measured Figs. and 
11. Ona separate scale Fig. are also plotted the values the quantity 
determined differentiation the fringe-position curves. From the 
values (also and p/5 from separate measurements appropriate 
directions), the strains can found previously described. 
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From series measurements across the horizontal axis Fig. 10, the 
value determined number points. These values, com- 
bined with’ those for trye from Fig. and with the modulus elasticity 
(also determined the moiré method), allowed computation and 
and respectively this case), the results which are shown Fig. 13. 
Also shown are the results theoretical analysis the disk when 
subjected concentrated loads. The difference between the theoretical and 
experimental curves can attributed mainly the flattening the disk 
the flat plates used apply the load. 
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ARCH DAM ANALYSIS WITH ELECTRIC ANALOG COMPUTER 


Richard MacNeal! 


SYNOPSIS 


electrical analogy that simulates the elastic deformations arch dams 
described herein. This analogy consists resistors, transformers, and 
current generators. The results analog computer analysis the Steven- 
son Creek Dam are compared with field measurements. The results de- 
sign study the Blue Ridge Dam are summarized. 


INTRODUCTION 


The determination the stresses and deformations anarch dam has long 
been one the most difficult analysis problems civil engineer- 
ing. The difficulties that are encountered arise principally for the following 
reasons: 


Elastic deformations play animportant role the determination stres- 
ses arch dams. 

Many structural degrees freedom must includedin arealistic anal- 
ysis order several different kinds deformation. 

Each arch dam requires separate analysis due wide variations 
canyon profile, thickness, and forth. 


The most elaborate and probably the most accurate method for analyzing 
the stresses anarch dam the Trial Load this method the dam 


open until January 1961, extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE, This paper part 
the copyrighted Journal the Engineering Mechanics Division, Proceedings the 
American Society Civil Engineers, 86, August, 1960, 
Mgr., Engrg. Service Computer Engrg. Pasadena, Calif, 
“Trail Load Method Analyzing Arch Dams,” Boulder Canyon Reports, Part 
Bulletin 1938, 
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replaced intersecting gridwork horizontal arches and vertical can- 
tilevers eachof whichresists extension, bending, shear, and twisting. The so- 
lution obtained adjusting the forces interaction between arches 
tilevers asto equalize radial, tangential, and rotational displacements. Be- 
cause the large number equations required trial load analysis, com- 
putations employing desk calculators have proven extremely time con- 
Adaptations the method for solution electronic digital computers 
should greatly increase its speed and general utility. 

Since presently available metheds analyzing arch dams are either quite 
laborious over-simplified, many designers prefer the use scale models 
instrumented with strain gages. Considerable success has been achieved 
Europe with this 

The method described this paper makes use electrical analogy 
which passive electrical elements (resistors and transformers) simulate the 
elastic and load transforming properties the structure. Such analogies have 
been successfully employed for solving many other types structural problems. 
Perhaps the earliest practical application electrical analogies structural 
problems was the solution truss and framework problems Vannevar Bush 
the last decade the method has been extended tothe static and dy- 
namic analysis and shells.12 Much use has been 
made these analogies aircraft structural 

When regarded tools for engineering computations, physical analogies 
are hybrid forms that have some the characteristics scale models and 
some the characteristics conventional mathematical analysis. 
testing, the answers are obtained from physical measurements, the essential 
distinction being that the quantities measured analogy are different 
physical type than those the prototype system. The techniques measure- 
ment and the problems associated with measurement are similar. impor- 


tant difference between scale models and physical analogies that scale 


“Arch Dams: Design and Observation Arch Dams Portugal,” Rocha, Serafim, 
and Silveira, Proceedings, ASCE, Vol, 82, No. June, 1956, 
“structural Analysis Electric Circuit Analogies,” Bush, Journal Franklin 
Inst., 217, March, 1934, 289. 
“Beam Vibration Analysis with the Electric Analog Computer,” McCann and 
MacNeal, Journal Applied Mecahnics, Transactions, ASME, Vol, 72, 1950, 
“An Equivalent Circuit for Vibrating Beam which Includes Shear Motions,” 
Trent, Journal Accoustical Soc. America, Vol, 22, 1950, 
“An Improved Electrical Analogy for the Analysis Beams Bending,” 
Russell and MacNeal, Journal Applied Mechanics, Vol, 20, 1953, 349. 
ASCE, Vol, 119, 1954, 
Solution Elastic Plate Problems Electrical Analogies,” Mac- 
Neal, Journal Applied Mechanics, Transactions, ASME, 73, 1951, 59. 
“Equivalent Plate Theory for Straight Multicell Wing,” Benscoter and 
MacNeal, NACA 2786, September, 1952, 
“Electrical Analogies for Stiffened Shells with Flexible Rings,” MacNeal, 
NACA 3280, December, 1954, 
MacNeal, McCann, and Wilts, Journal Science, December, 1951, 
“Analysis Multicell Delta Wings Caltech Analog Computer,” Mac- 
Neal and Benscoter, NACA 3114, January, 1954, 
“The Integrated Use Analog and Digital Computing Machines for Aircraft Dyna- 
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model can constructedeven when knowledge the processes involved the 
prototype structure incomplete, whereas the construction analogy re- 
quires complete description the physical laws and properties the pro- 
totype. Thus, scale model, usually sufficient touse similar materials 
and produce dimensions and applied loads whose ratios are accordance 
with principles dimensional analysis. physical analogy, the other hand, 
analogy only virtue similarities its physical laws with those 
the prototype system. fact must possible describe the behavior 
the analogy set mathematical equations that essentially identical 
the set satisfied the prototype system. this respect, computation phys- 
ical analogy similar computation mathematical analysis. 

Thus, problem can solved physical analogy, can also solved 
direct mathematical analysis. For this reason analog computer more 
frequently regarded aid mathematical computation than form 
model testing. While this view expresses intrinsic truth, should real- 
ized that the operation analog computer more closely related model 
testing than the operation digital computer. 

Computation with electric analog computer has certain advantages over 
both model testing and conventional mathematical analysis. The recording 
results more rapid and usually more accurate than model testing due 
special measuring equipment built intothe computer. The which 
the analogy constructed are permanent and re-usable and can easily ad- 
justed any The main advantages over digital computation are 
the speed obtaining results and the fact thata great deal less algebraic ma- 
nipulations are required the preparation problem. 

Along with these advantages, analog computing has its own special difficul- 
ties. Unlike general purpose digital computer, given analog computer 
restrictedas the kinds problems that can used solve. The type 
computer described this paper mainly useful for determining the static 
and dynamic response structures and for the solution related mechanical 


problems. Other analog computers are designed solve heat flow problems 
control system problems. 


Accuracy limitations are often considereda disadvantage analog comput- 
ers. While true that precision comparable digital 
achieved, the degree accuracy usually sufficient for most engineering ap- 
plications. Another important disadvantage the analog computer, with respect 
digital computation, that must re-assembled for each problem that 
analyzed. Although the use removable patch boards can reduce this oper- 
ation few hours, the set-up time incomparably greater than that 
card- tape-programmed digital computer. 

The advantages claimed for should particularly signif- 
icant for the analysis arch dams. Since the numerical computations are 
lengthy and the construction and testing accurate model painstaking 
operation, the analog computer should enjoy speed advantage. Since many 
design changes are required achieve efficient design, the flexibility 


the analog computer with regard parameter changes should very useful 
feature. 
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ELECTRICAL ANALOGIES FOR THE STATIC ANALYSIS 
ELASTIC STRUCTURES 


the type electrical analogy described herein, current proportional 
force proportional displacement. Since flexibility the ratio 
displacement force, and resistance the ratio voltage current, 
follows that electrical resistance analogous elastic flexibility, 
trated Fig. 

very important component the electrical analogies elastic structures 
the ideal transformer. Although this element has many uses, commonly 
considered the analogy ideal, weightless lever, illustrated 
Fig. There fundamental relationships for anideal transformer. The 
first that, for different coils wound the same magnetic core, the voltage 
across eachcoil proportional tothe number turns the coil. The second 
that, when the current flowing each coil multiplied the number 
turns, the sum such products zero. These relationships are illustrated 
Fig. 2(b). The winding sense each coil indicated plus and minus 
signs. The positive direction current flow from the positive toward the 
negative terminal the winding that, the figure, the current the wind- 
ing with turns while the current the winding with turns 
The relationship expresses the fact that the sum currents 
flowing into the junction with voltage zero. 

From the relationships given Fig. the transformer-turns 
ratio analogous the lever ratio 

The relationships for the ideal lever may combined show that the me- 
chanical work done the external forces zero. Similarly the relationships 
for the ideal transformer may combined show that the electrical power 
consumed the transformer zero. These statements express the intrinsic 
meanings the terms ideal lever andideal transformer. Practical levers and 
practical transformers are not ideal. major portion the expense 
analog computer the type considered here involved the construction 
transformers that satisfy the relationships ideal transformers with sufficient 
accuracy. These relationships are not even approximately satisfied very 
low frequencies very high frequencies. For this reason the computer 
usually operated frequencies the medium audio range Since 
resistors are the only impedance elements, direct current concepts used 
connection with the circuits. 

addition resistors and transformers, electrical analogies for static 
structural analysis must include devices that simulate applied loads. Since 
current analogous force, the analogyfor aconstant loadis constant cur- 
rent generator. Such device may approximated very high voltage 
source feeding through very large resistance. 

Resistors, transformers, and constant current generators are the only ele- 
ments required the analogy for statically loaded structure. The manner 
which these elements are interconnected form the analogy for complete 
system illustrated Fig. good method for constructing analogy 
complete system consists the following steps: 


Subdivide the mechanical system into simple elements. 

Construct electrical analogies for these elements from standard list 
analogies. the list does not contain analogy for element, write the 
relationships satisfied the element and devise electrical circuit that sat- 
isfies these relationships. 
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FIG, BETWEEN FLEXIBILITY AND RESISTANCE 


(b) IDEAL TRANSFORMER 


FIG, 2,—ANALOGY BETWEEN IDEAL TRANSFORMER 
AND IDEAL LEVER 
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Interconnect the elemental electrical analogies satisfy continuity re- 
lationships. 


the example Fig. the elements the mechanical system are rigid 
lever, two springs and Electrical analogies for these elements have 
been previously derived. The continuity relationships satisfied are that 
the relative displacements the elements are zero the points where they 
are joined together and that the forces acting onthese joints are equilibrium. 
The displacement conditions imply the equilibrium conditions and vice versa 


Rigtd Bar 


(a) MECHANICAL SYSTEM 


Ground 
Connections 


(b) ELECTRICAL ANALOGY 


FIG, 3,—ELECTRICAL ANALOGY FOR SIMPLE MECHANICAL SYSTEM 


since work done the joints. Thus, since electrical connection con- 
sumes power, either usedas basis which intercon- 
nect the elementary electrical analogies. The displacement condition some- 
what easier apply the example Fig. since merely states that the 
voltages the electrical elements are equal points corresponding the 
joints the mechanical system. examples where the forces joint are 
not co-linear, transformers may required affect the interconnection and 
such cases the equilibrium conditions are frequently easier apply. 

The bending beam provides somewhat less trivial example. this 
example, illustrated Fig. the first step replace the continuously elas- 
tic beam setof lumped springs, concentrated loads and simple rigid bars. 
The only newelement occurring inthis idealized lumped structure the rigid 
bar for electrical analogy derived With this new elemen- 
tary analogy, the principles outlined above easily lead tothe complete electrical 
analogy shown Fig. 4(c). this circuit, mechanical symbols are used 
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(a) UNIFORM CANTILEVER BEAM WITH DISTRIBUTED LOAD 


(b) EQUIVALENT LUMPED MODEL 


(c) ELECTRICAL ANALOGY 


FIG, ANALOGY FOR CANTILEVER BEAM 
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FIG, 5,—ELECTRICAL ANALOGY FOR RIGID BAR 
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FIG, 6.—IDEALIZED MODEL STEVENSON CREEK DAM 
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designate the values electrical elements economy measure. The ac- 


tual electrical values will equal tothe mechanical values multiplied scale 
factors. 


The subsequent development the theory electrical analogies for stat- 
ically loaded elastic structures consists mainly elaboration the tech- 
niques previously described with more complex structural The el- 


ementary analogies required inthe analysis arch dams will derived sub- 
sequently. 


STRUCTURAL MODEL USED THE ANALYSIS 
CYLINDRICAL ARCH DAM 


The arch dam may characterized relatively thick shell homoge- 
neous elastic material. will assumed that displacements vary linearly 
across the thickness the dam, which assumed sufficiently thin for 
the middle surface taken the neutral surface for bending, but not 
thin that any components shearing strain may neglected. Because there 
are finite number lumped circuit elements inthe computer, finite difference 
approximations are required analysis. Furthermore, since resistors 
are linear elements, assumptions linearity are required unless expensive 
nonlinear elements are available large quantities unless iterative solu- 
tion techniques are employed. 

order simplify the presentation, will assumed that the middle 
surface the dam cylinder. This assumption was not made the anal- 
ysis the Blue Ridge Dam. 

The first the idealization the structure lay out gridwork 
horizontal and vertical lines the middle surface the dam illustrated 
Fig. The dimensions shown refer the Stevenson Creek Dam. 


The dam then replaced idealizedstructure with three kinds elements. 
These elements are: 


Horizontal arch elements, that resist axial thrust, bending about ver- 
tical axis, and shear horizontal plane. 

Vertical cantilever elements, that resist vertical thrust, bending about 
horizontal axis tangent the neutral surface, shear radial plane, and 
shear plane tangent the neutral surface. 


Torque blocks, that occupy the panels between arch elements and canti- 
lever elements, and that resist twist. 


The loads distributed over the surface the dam are replaced concen- 
trated loads acting the arch-cantilever intersections. order that such 
loads may resisted without inducing local bending stresses the arches, 
the curved archelements are replaced straight chords between intersection 
points. This substitution may justified finite-difference analysis. 

The lumped model for arch elements shown Fig. which also illus- 
trates the manner variable radius curvature may treated. The 
bending stiffness the arch represented rotational springs. The 
flexibility each spring computed the formula 
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FIG, MODEL FOR ARCH ELEMENTS 
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The limits integration are from the center one arch element the center 
the next arch element. the moment inertia the arch cross section 
and equal 1/12 Ah, where the thickness the dam and the 
vertical dimension the arch. The thrust and shear flexibilities respectively 
arch element are given the formulas 


shear correction factor, approximately equal 0.8, and the cross 
sectional area the arch. 

With regard vertical cantilever elements, the usual assumption the trial 
load method, that the cantilevers are inextensional for vertical thrust, will 
This assumption has the effect decoupling the effects vertical and 
horizontal loads that only horizontal loads need included the analysis 
considered here. The assumption also has the effect making shear stresses 
the neutral surface functions only the relative tangential motions ad- 
jacent Thus, the resistance tangential shear deformation may 
represented simple shear springs connected between adjacent arch-canti- 
lever intersections shown Fig. 

The resistance bending and shear radial plane provided 
lumped beam model similar that shown Fig. 4(b). 

Neither the arch elements nor the cantilever elements resist twist. Re- 
sistance twist provided third type ofelement, referred torque 
block, which occupies the panel defined arch and cantilever centerlines. 
Twisting the dam may regarded warping the plane connecting the 
four corners sucha panel caused concentrated loads acting the corners. 
These corner loads may replaced couples acting the top and bottom 
faces the torque block shown Fig. The block then replaced 
torsional spring with flexibility 


This spring joined adjacent arch elements their centers. un- 
necessary add additional springs connecting adjacent cantilevers since the 
twisting deformation the dam completely defined the relative rotation 
adjacent arches. Although the method treating twist described here cor- 
rectly accounts for the strain energy twist and leads simple electrical 
analogy, introduces ficticious components radial shear into the arch and 
cantilever elements that should subtracted fromthe measured 
the analog computer. 

The structural model that has been described more-or-less equivalent 
that employedin conventional trial load analysis. The principal differences 
are inthe manner treating twist and inthe substitution concentrated loads 
for loads continuously distributed alongthe arches. viewpoint shell 
theory, the main objections the model are the thrust deforma- 
tion the cantilevers and the omission Poisson’s ratio coupling. Both 
these effects canbe the expense additionalelements inthe model. 
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Since the number electrical elements available analog computer 
limited, the inclusion these relatively unimportant effects would decrease 
the number grid points that could represented. 

The concentrated loads applied arch-cantilever intersections require spe- 
cial attention. Due the uniform character the loads sustained arch 
dam, small percentage change the load given point may produce 
disproportionately large change the local bending moment. For most load ef- 
fects the lateral load applied given point may regarded the product 
pressure and surface area. The value surface area given arch- 
cantilever intersection should chosen such that, anidealizedcircular arch 
model consisting chord members, uniform radial pressure gives uniform 


(a) ELEMENT WITH CORNER (b) ELEMENT WITH COUPLES 
LOADS HORIZONTAL FACES 


FIG, BLOCK 


stress equal r/t the absence rib shortening. This condition satis- 
fied point the model Fig. the radial load equal to: 


must equal for all points the same arch. 

The effects due atemperature rise, AT, uniform across 
the thickness the dam, may simulated equivalent radial and tangental 
loads. The basis this equivalence fact that the total strain the cir- 
cumferential direction, may written 


6 
| 
1 | 
Rae 
ag 
4 
a 
| 
| 


DAM ANALYSIS 


139 


where the actual tensile stress and the coefficient thermal expansion. 
thickness, constant radius arch, this applied stress may produced ra- 
dial pressure equal more general cases the equivalent ther- 
mal loads applied point Fig. are: 


After the deformations and stresses due this loading are measured, 
computed, the quantity must subtracted from the resulting tensile 
stress the circumferential direction. 

The effects uniform temperature gradient through the thickness the 
dam may simulated applied moments vertical axis which may, 


turn, replaced distribution point loads with predominating tangential 


ELECTRICAL ANALOGY FOR CYLINDRICAL ARCH DAM 


The remaining steps the construction analogy are draw 
electrical circuit analogies for these elements and then toconnect the circuits 
together satisfy the conditions continuity. 

The elements the arch dam and their corresponding electrical analogies 
are shown Fig. 10. The analogy for arch chord element contains re- 
sistor simulating the thrust flexibility and separate circuit simulating the 
shear and rotational properties the element. This latter circuit similar 
the circuit for rigid bar shown Fig. except for the addition re- 
sistor simulating the shear flexibility given Eq. 

The special circuit Fig. 10(c) required satisfy the relationships 
existing joint between arch elements such joint Fig. Since 
work done the joint, either the displacement conditions the equilibrium 
conditions may usedin the derivation the analogy. The displacement con- 
ditions will Transformed displacements are defined follows: 
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(a) ARCH CHORD ELEMENT 
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FIG, FOR THE ELEMENTS ARCH DAM 
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and, 


These relationships are satisfied the circuit Fig. 10(c), provided that the 


transformer turns ratios are set the following values for transformer 


The scale transformation Eq. must also applied the arch chord ele- 
ments. 

The current generators simulating radial and tangential loads are connected 
nodes where the voltages are proportional and The circuits for 


the other elements Fig. make use previously derived analogies. 

The manner which the elements are connected together illustrated 
Fig. 11, which shows the complete analogy for very simple model consisting 
two arches and one cantilever. The foundation assumed infinetely 
rigid that all components displacement are zero the abutments. This 
condition simulated grounding all electrical elements points the 
abutments. The circuit shown four separate diagrams, one diagram for 
each type coordinate represented the analysis. Expressed inelectrical 
terms, the conditions displacement continuity are that, points where ele- 
ments are connected together, voltages corresponding the same type co- 
ordinate are equal. These conditions are satisfied solid electrical connec- 
tions between corresponding nodes. 

The electrical analogy for arch dam model containing many arches and 
many cantilevers similar form Fig. 11, but contains many more elec- 
trical elements. Special problems encountered practical analyses will 
discussed subsequently. 


ANALYSIS THE STEVENSON CREEK DAM 


The Stevenson Creek Dam small test dam built 1926 specifically for 
the purpose obtaining field measurements the stresses and deformations 
arch dam. complete account the construction and testing given 
The dimensions the dam are shown Fig. where may 
seen that the abutment profile almost exactly symmetrical 45° The 
dam was instrumented the downstream face with gages that measured the 
horizontal and vertical components strain. The small reservoir could 
emptied filled few hours that fairly reliable measurements the 


“Arch Dam Investigation,” Report ASCE Committee, November, 
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strains due water load were obtained. The dam was filled progressively 
higher levels the tests proceeded. Large vertical cracks appeared the 
top and bottom the crown section when the water level reached ft, the 
cracks extending about ftup from the bottom and down from the crest. 
The cracks did not appear grow whenthe water level was raised the crest 
the dam (60 ft). 

The Stevenson Creek Dam was selected for investigation onthe analog com- 
puter because the availability carefully recorded field measurements, the 
simplicity the cylindrical structure, and the general the dam 
test structure. The crown section was considered plane symmetry 
sothat arepresentation only one half the dam was the analysis. 
shown Fig. the idealized model the half-dam contained six arches 
and three vertical cantilevers. The number cantilevers was made purposely 
conditions where these effects would large. 

complete description the analog computer investigation given else- 
The electrical analogy contained total transformers and 140 
resistors. large number cases were analyzed corresponding different 
water levels, different temperature conditions, different elastic modulifor the 
foundation material, and different conditions cracking along the abutments 
and along the crown section. 

The effects uniform temperature change were examined the equiva- 
lent load method discussed previously. The presence the foundation was 
simulated shear, bending and thrust springs computed the method 
These springs were replaced resistors which were added 
series with corresponding adjacent elements the dam. 

The effects cracks could only guessed. the case the crack the 
top the crown section, was that this crack did not close until 
the dam was nearly filled with water. For this reason the thrust stiffness 
the top arch (B) the centerline was made thrust stiffness the 
centerline was given reducedvalue. Arches and were given 
zero bending stiffness the centerline. 

Results the cracked dam simulation with full head (60 ft) water and 
temperature change are shown Figs. 14, where they are com- 
pared with stresses derived from field measurements. For this comparison 
was assumed that the elastic modulus the concrete used the dam 
2.9 106 psi, that the elastic modulus the foundation, E,is 1.6 106 psi and 
that Poisson’s ratio, the dam 0.15. rather than stresses 
were measured the field, values andy are required for the reduction 
stresses from field measurements. 

Fig. will observed that the peak horizontal stress the down- 
stream face obtained the analog computer about 100 psi less than the peak 
stress obtained from field measurements. There reasonto believe that part 
this discrepancy due tofinite difference approximation inthe arch bending 
model which increases twisting moments the expense arch bending mo- 


ments. This approximation canbe greatly reduced arefinement the arch 
bending model. 


“analysis Stevenson Creek Arch Dam,” MacNeal and Hill, Com- 
puter Engrg. Assocs, report project REE 148, March, 1959, 


“Uber die Berechnung der Fundament deformation,” Vogt, Anhandlinger det 
Norske Videnskaps-Academi, Oslo, 1925, 
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field measurements were made shear strains. The analog computer 
results given Fig. show shear stresses along the abut- 
ments. These shear stresses measured horizontal-vertical coordinate 
system indicate tendency the dam bend about axes that were coincident 
with the abutments. 

The contour plots stresses were obtained from analog computer results 
the following steps. Electrical currents were converted into internal loads 
and moments the idealized model. These loads and moments were converted 
into stress components the faces the dam specific points. Contribu- 
tions stress component from different sources, such horizontal stress 
due arch thrust and arch bending, were interpolated the same points and 
added, Finally contour lines were drawn with these values bench marks. 


DESIGN STUDY FOR THE BLUE RIDGE DAM 


Sketches the Blue Ridge Dam are shown Fig. 15. The canyon profile 
distinctly unsymmetrical with afairly steep slope the right abutment and 
relatively gentle slope the left abutment terminating thrust block. 
The radius curvature varies from 223 the crest less than 100 
stream bed. The arch barrel exceptionally thin with ratio base width 
height less than 10%. order eliminate extensive undercutting due 
change radius with elevation, the crown section inclined downstream 
average angle about 15°. 

The distortions the geometrical shape the dam from true cylinder 
create problems analysis not previously mentioned this paper. The first 
problem the matter defining the centerlines cantilevers the idealized 
model, This problem was solved passing planes through vertical refer- 
ence axis arbitrarily located near the “average” center curvature the dam 
shown Fig. 15. Intersections these planes with the arch centerlines 
define the arch-cantilever intersection points. elevation diagram was then 
prepared projecting these points cylinder with the reference axis 
center, shown Fig. 15. this diagram arches and cantilevers are hori- 
zontal and vertical lines. may seen, total six arches and nine can- 
tilevers were used the analysis. 

The second problem the determination the principal planes bending 
have infinite bending stiffness about axis perpendicular the face the dam 
and infinite stiffness about axis tangent the face the dam. non- 
cylindrical dam the direction the normal changes along each cantilever 
that the principal planes bending are rotated the elevation changed. 
This problem was solved changing the directions the principal planes 
arch-cantilever intersection points, that each cantilever element remained 
untwisted. The tangential and radial displacements cantilever element 
were then related the tangential and radial displacements the adjacent 
arches the electrical analogy means additional taps the transformers 
Fig. 10(c). 

Perhaps the most important problem the resolution applied loads into 
components. Since the cantilevers are inclined they cannot resist vertical loads 
without bending. The crown section is, fact, statically unstable under the 
action gravity. This problem was handled resolving applied loads into 
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horizontal components and components tangent the axes the cantilevers. 
The horizontal components defined were assumed resisted the lat- 
eral load-carrying capacity the dam whole, including arch bending, 
cantilever bending, and on. 

The components tangent the axes cantilevers were assumed car- 
ried cantilever thrust without bending. These assumptions appear reason- 
able long the radius curvature vertical plane large compared 
the radius curvature horizontal plane. 

The electrical circuit analogy for the dam was similar form that 
Fig. and employed approximately 140 transformers, 250 resistors, and 
current generators. inthe analysis the Stevenson Creek Dam foundation 
flexibility was simulated the method indicated Vogt.18 The foundation 
modulus for thrust was made less than the foundation modulus for bending 
order approximate the effect plastic deformation. 

The primary objective the study find asatisfactory design which 
allowable stresses were not Subsidiary objectives were investi- 
gate the mechanisms load distribution and analyze the final design con- 
figuration for variety loading conditions. Complete sets stresses and 
deflections were obtained for total cases. 

the first seven cases the elastic model the dam was built adding 
additional class elements each successive case. Thus arch thrust 
flexibility, cantilever bending stiffness, tangential shear stiffness, twisting 
stiffness, and foundation flexibility were successively added. Besides pro- 
viding simplified models that could easily checked, these cases gave valu- 
able insight into the importance each type element and the mechanisms 
load transfer generally. 

The next twenty cases were usedin the investigation design modifications. 
the initial design, extremely high cantilever bending stresses were observed 
the base the dam. was soon realized that these high stresses could not 
eliminated either reasonable changes the thickness the dam 
minor changes its shape. Major improvements were obtained introducing 
slip plane midway between arches and (Fig. 17) and hinging the dam 
along the abutments. After these changes were made, the region maximum 
cantilever bending stresses was shifted the neighborhood the intersection 
arch andcantilever No. were further reduced chang- 
ing the thickness distribution the top half the dam. Satisfactory stress 
levels were finally achieved eliminating the effects temperature drop. 
The progress made reducing cantilever bending stresses indicated the 
bar chart Fig. 16. 

The last nine cases were used for the investigation minor design alter- 
natives and for the analysis additional loading conditions such water loads 
for partly filled reservoir and equivalent earthquake loads. The principal 
stresses the upstream and downstream faces the final dam design are 
shown Fig. 17. The radial deflections for this case are shown Fig. 18. 
study the Blue Ridge Dam occupied the analog computer for 
total time ten days. 
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FIG, 18,—RADIAL DEFLECTIONS FOR THE FINAL DESIGN 
THE BLUE RIDGE DAM 


ACKNOWLEDGMENTS 
The arch dam design study for the Blue Ridge Dam was performed for the 


firm Leeds, Hill and Jewett Los Angeles, Calif. complete report the 
results this study given elsewhere. 


Study the Blue Ridge Dam,” MacNeal and Hill, Computer 
Engrg. Assocs. report project ESS 143, March, 1959, 


| | | 
| | | 
} | 
| | 
| | 
| | 
| | =i 
L 
+ 
} 
é 
be 


Bell 
il 
: 
4 


2580 August, 1960 


Journal the 
ENGINEERING MECHANICS DIVISION 


Proceedings the American Society Civil Engineers 


ELASTIC MODEL DESIGN THE B-58 AIRPLANE 


SYNOPSIS 


This paper outlines the extensive elastic model test program used the 
design highly redundant structure the B-58 airplane. Considerations 
the design the various models the program are listed and several unique 
methods obtaining and presenting test data are shown. Scale factors for 
both structure and loadings are developed detail the Appendix. 


INTRODUCTION 


The elastic model program used the stress analysis the B-58 bomber 
Convair-Fort Worth, Division General Dynamics Corporation, repre- 
sents one the most extensive programs its kind ever attempted. was 
started time when other available methods for analyzing highly redundant 
structures had failed completely producing the degree accuracy required 
this radical, long range weapons system with its minimum structural weight 
requirements (Figs. and 2). 

The latest model this program now being used the evaluating 
number design improvements the production B-58 airplane. This model 
represents the last four distinct model configurations used, since 1952, 
evaluating more than 300 full airplane loading conditions (Fig. 3). 


Note.—Discussion open until January 1961. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Engineering Mechanics Division, Proceedings the 
American Society Civil Engineers, Vol. 86, No. August, 

Sr. Struct. Engr., Convair-Ft. Worth, Division General Dynamics Corporation. 
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DEFINITION ELASTIC MODEL 


elastic model representation that structurally simulates its proto- 
type one more important respects, and has the necessary provisions in- 
cluded for simulating predicted loadings, usually hydraulic rams lead 
weights. Additional provisions are commonly made the form strain and 
dial gages for determining stresses and deflections significant points during 
test loadings. 

The size the model, degree simulation, extent instrumentation and 
number conditions tested are strictly function the accuracy needed 
and time and money available for the program. Scale models range size 
from small 1/100 for small plastic models full size the case 
for static test airplanes. 


OTHER TYPES SCALE MODELS 


number other types models are used the developmental phases 
airplane. Several these are the force model, pressure model, and flutter 
model. Each these true aerodynamic model and represents the outside 
configuration its prototype high degree. While the force model (Fig. 
determines only value total net airload for given attitude, the pressure 
model yields refined pressure distribution over its surface. The flutter 
model simulates the external airplane contour, its elastic properties, and its 


internal mass distribution the evaluation certain aeroelastic phenomena 
(Fig. 5). 


NEED FOR THE ELASTIC MODEL PROGRAM 


the initiation the elastic model program 1952, the most advanced 
analytical method available for the analysis redundant structures consisted 
attempted hand solutions structural matrices near order 30. Later, the 
IBM-701 high speed digital computer became available for use this project. 
The largest passive analog computer available was the SES Analog installation 
the California Institute Technology. However, all attempts using both 
these techniques over-allanalysis the wing proved futile. Where cer- 
tain simplifying and overlapping assumptions would have been possible for 
different type weapons system, the B-58 with its maximum range concept 
demanded absolute, unyielding loyalty tothe task eliminating all unneces- 
structural weight. This could only attained devising more rigor- 
ous and extensive analysis technique for redundant structures than had previ- 
ously been attempted. 

Because the success achieved using experimental methods other prob- 
lems which had defied analytical solution, was decided attempt experi- 
mental technique this program. Therefore, feasibility study was started 
outline the requirements the initial model configuration. 

This experimental model technique has proved successful that the first 
crudely designed 1/4 scale model has been succeeded three others, each 
incorporating extensive improvements degree simulation 
instrumentation. The last model, the 3/8 scale, has provisions for reading 
strains 1,826 locations and deflections 170 locations during given test. 
Virtually allextensive data reduction, such principal stress determinations 
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and internal station balances, are coded for solution the computer installa- 
tion. 

Because this latest model incorporates all the refinements its predeces- 
sors, the balance this paper will concerned with explaining only the de- 
tails the 3/8 scale elastic model program. 


DESIGN AND INSTRUMENTATION 


Design 3/8 Scale Elastic Model.—The unusual 3/8 scale this final 
having scales 1/4 and 1/2 respectively. After two the earlier 1/4 scale 
models had been built became apparent that they were too small provide 
space for refined representation small fittings and installation necessary 
strain gages. Because the 1/2 scale model, built represent only the outer 
portion the wing, had proved more than adequate size and costly 
construction was decided proceed with 3/8 scale representation. After 
more than three years testing this model, its scale selection has been 
fully justified. One consideration keeping model’s scale small pos- 
sible the reduction the amount simulated load that must applied for 
each test condition. should noted that keeping all other considerations 
equal, increase model scale factor 2.0 increases the amount 
applied load for given condition factor 4.0. The necessity applying 
additional load this magnitude during test canbe very time consuming and 
costly. Therefore, any increase model scale should fully justified before 
being approved. 

Model Scaling Factors.—Requirements and assumptions that must con- 
sidered any general scaled prototype representation and those considered 
the 3/8 scale elastic model program are included detail the Appendix. 

Structure Simulated The Model.—Because the elastic model program 
was started primarily for use tool wing stress analysis, great de- 
gree refinement was simulation other airplane components. 
Photographs the earlier models show nacelle vertical tail represen- 
tation and show only crude box structure used simulate fuselage stiffness 
(Figs. and 7). 

However, became evident early the program that some refined simu- 
lation other components would justified because their elastic properties 
affected wing internal loads and the model presented economical method 
obtaining their design stresses well those for the wing. Again, study 
the four models the program indicates successive increase the de- 
gree component simulations each. 

The 3/8 scale model shows the greatest development this concept all 
component representation. found model that simulates virtually all 
the wing, fuselage, and nacelle structure. Even the vertical tail simulated 
near its root insure proper load introductions the aft fuselage structure 
(Figs. and 10). The only remaining components not simulated are the 
landing gear and elevous. However, the computed loads these components 
are introduced concentrated loads their simulated points attachment 
the model. 

Although the earlier models simulated the stiffness several structural 
members only one model member, this lumping has been almost completely 
eliminated the latest design. Although lumping stiffnesses members and 
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their external loadings presents difficulty, the unlumping the resulting 
data determine individual prototype member loads next impossible. 

Load fittings are provided the model accommodate every major point 
load introduction the prototype such pod hooks, landing gear fittings, 
and others. addition, load points are provided virtually every structural 
intersection for the application needed loads and moments. This extensive 
load introduction system the model permits application scaled prototype 
distributed loadings with minimum loss accuracy due iumping. 

Representation Unconventional Prototype Structure.— number ad- 
vanced manufacturing and design methods used the prototype could not 
incorporated the model because their unavailability the time model 
design the excessive cost tooling for scaled representation. One 
these was the corrugated spar and bulkhead webs the wing substructure. 
These corrugations were designed replace the conventional vertical stiff- 
eners previously used the prevention web buckling. Because their exact 
representation was not feasible, was decided use conventional flat web 
stiffener combinations the model (Fig. 11). However, the longer length 
curved web per inch spar onthe airplane had considered flat- 
web design such that correct representation web shear deformation per 
inch spar the model would result. 

This was done substituting model webs having thickness factor in- 
versely proportional the web length “growth factor” the prototype. This 
resulted thinner web shorter length operating higher stress than 
the longer airplane counterpart. The correct amount web shear deformation 
per unit length spar was obtained. Because this one cases where 
model stresses are not directly applicable those the airplane, the data 
obtained are presented terms equivalent airplane load pounds. This 
especially useful the few cases where airplane webs have been replaced 
diagonal trusses after the model design was completed. 

decisionto simulate the structural sandwich the airplane wing surfaces 
would have been entirely impossible when model was designed, 
because bonded sandwich panels were still the primitive stages develop- 
ment. Even later direct scaled representation the model would have been 
prohibitive both cost and time. Therefore, these surfaces were simulated 
conventional flat skins and slugs (Fig. 11). The centroids the skins 
the located closely possible tothe comparable surface panel 
centroid the airplane. The simulation the internal concentrated slug 
bending material inthe airplane panels was included spar and bulkhead caps 
and external slugs where model cap areas became excessive. 

Where feasible, exact airplane methods construction are incorporated 
the model. the case the outboard pylon plate, which provides attachment 
the wing and nacelle, welded steel structure provided simulate 
similar steel forging used the prototype (Fig. 12). 

Instrumentation Installed The order provide internal load 
distributions for its prototype, the model includes one the most extensive 
installations strain gages ever attempted (Fig. 13). These gages provide 
measurements axial, shear, and combinations axial strain various lo- 
cations the model depending the type gage and its installation. spar 
and bulkhead flanges and fittings requiring only uniaxial stress determinations, 
approximately 550 axial gages are installed. spar and bulkhead webs and 
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certain surface material for which only shearing stresses are re- 
quired, 496 V-setts gages have been installed. The most refined strain deter- 
minations are obtained model surface panels where 260 rosette gages are 
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mounted. The latter gages measure three values strain 45° each other. 
These values strain point are used determine the complete state 
stress using equations graphically represented Mohr’s Circle. 
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Wherever possible, each gage installed the model “backed” and wired 
into the recording instrument eliminate the effect local plate buckling 
from the recorded strains. 

One the most serious problems that has existed throughout the model 
program the evaluation long time drift due temperature and other ef- 
fects the measured strain data. The temperature effects are notone heat 
being applied part the test, but rather case the model changing 
temperature due surroundings during the eight more hours between the 
initial and final data recordings for test. 

Unfortunately, the unaided strain gage reflects total strain experienced 
the part mounted and has noway differentiating between strain 
produced loading and that produced thermal expansion. avoid this, 
each active gage the model wired similar dummy gage mounted 
unloaded piece identical material nearby. The wiring such that the 
strain the dummy gage (which wholly thermal nature) subtracted 
from the total reading the active gage (which both thermal and loaded 
nature) yielding only the strain produced load. assumed that the ma- 
terial the model and the material for the dummy gage will change tempera- 
ture uniformly during the test. 

The following additional provisions are made either minimize, help 
evaluating these thermal effects: 


The model mounted air conditioned test laboratory where air 
temperature held constant possible. 

wall has been built around area protect from 
the blasts hot cold air occasioned opening doors permit entrance 
equipment personnel other areas the test lab during model test. 

Both the initial and final readings for each gage are recorded with only 
tare load applied. Any difference between these readings indicates some sort 
long time drift the strain recording setup. Most this probably caused 
temperature change the model structure between these readings. 


Data Recording Equipment Model.—The strains and deflections the 
elastic model are recorded number different types instruments de- 
pending their locations and the accuracy required. Deflections near the 
root the wing and along the centerline the fuselage are measured dial 
gages that canbe read the nearest 0.001 in. these gages provide maxi- 
mum travel only in., deflections inthe outer portion the wing are meas- 
ured making observations with surveyor’s level toanaccuracy 0.01 in. 
(Fig. 13). 

Indicated strains are measured both manually read Brown Indicators and 
SADIC automatic recorders. Various switching techniques and patter box ar- 
rangements have been used expand the capacity the indicators where 
approximately 1700 channels instrumentation may recorded using only 
two 100 channel SADIC’s and Brown recorders. Although the model has 
1,826 channels installed instrumentation, condition yet tested has indi- 
cated pre-test severity justifying the reading more than 1500 channels. 
This elimination strain readings apparently non-critical areas has con- 
tributed greatly the over-all economy the program. 

The validity each recording instrument during test checked read- 
ing least one channel strain data each instrument from calibrated 
load bar. the end the test, each recorded value strain from the load 
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bar compared against the known value check the accuracy each 
instrument. number cases, faulty data from certain instruments have 
been corrected after the test and retests avoided using correction factors 
derived from this check information. 


LOAD PREPARATION AND TESTING 


Load condition preparation for the model involves selecting the most criti- 
cal conditions the airplane required sustain, determining the complete 
airplane loading for each, and converting these loads model magnitude, 
keeping mind both scale factors and model loading capabilities. 

The B-58 bomber designed operate allspeeds twice the speed 
sound altitudes above 50,000 under multiplicity gross weight and 
weapons configurations. This array possible conditions extensive that 
criteria group set determine the airplane’s capability and select 
the structurally significant conditions within the criteria approved the Air 
Force. After its selection, test condition outlined listing such items 
its Mach number, gross weight, altitude, angular and scalar accelerations, 
attitude, control surface deflections, and weapons configuration. 

This condition description used another group, charged with the re- 
sponsibility determining distributed loadings, guide combining 
unit data for the condition. This unit data consists predicted pressure dis- 
tribution over the airplane for certain types situations such unit change 
angle attack, elevator setting, yaw angle, and many others. After all 
necessary unit data have been combined dictated the condition descrip- 
tion, resulting pressure distribution for the entire airplane results. This 
then combined with inertia distribution consisting each item weight 
times its local acceleration, yielding perfectly balanced airplane loading. 
For all conditions, total inertia loads exactly balance total applied loads. For 
flight conditions the applied loads result from pressure distributions, whereas 
those for ground conditions originate either from landing gear drag chute 
loads. 

After the airplane distributed loading for the condition has been obtained, 
lumped model load introduction panel points and reduced magnitude, 
through the use model scale factors. This tabulation “panel point loads” 
released the elastic model test group. 

Pre-Test Condition Evaluation.—Because design stresses are urgently 
needed during the design stages airplane, efforts are made restrict 
consideration the more critical conditions for test the model saving 
both time and money. Therefore, soon test condition has been re- 
ceived, its external load summaries are rigorously compared against those 
for previously tested critical conditions. Unless appears comparatively 
critical, the condition either omitted from testing relegated the latter 
stages the test program. This severe scrutiny external loads for con- 
dition comparison with others from the same airplane sometimes points 
peculiarities that might indicate errors some stage the condition’s prepa- 
ration. course, similar checks both severity and validity are performed 
both the criteria and loads preparation groups while the condition being 
prepared. 

Preparation Pre-Test the evaluation ofa condition proves 
both valid and significantly critical, released the engineering test 
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lab for testing the model. that time, test lab personnel are notified 
expected critical areas the condition and its relative severity with respect 
other conditions currently scheduled for testing. This information used 
determining what instrumented areas will monitored during the test 
the limited strain recording equipment. 

While the actual rigging load points, wiring required strain gage chan- 
nels, and pretest operations are continuing the test lab, certain test facili- 
tation data are being prepared the elastic model group. The most important 
these are the “load sequence cards” which indicate the order loading panel 
points during given test. Because both personnel and equipment limita- 
tions impractical load all 340 model points simultaneously. These load 
sequence cards provide recommended order loading that maintains safe 
balance between and down loads all times. was learned early the 
program that any random loading model panel points would soon result 
disaster. addition maintaining approximate vertical balance, the load 
sequence cards also keep side, roll, yaw, and pitch unbalances within safe 
limits. 

Model Support Points.—The model supported only three points, which 
are instrumented indicate any additional them during testing. 
Center line supports are located the fuselage near the leading and trailing 
edges the wing. They include universal joint arrangement that prevents 
them from resisting either roll local pitching moment. The remaining 
support resists vertical load near the wing-leading edge airplane span sta- 
tion 56.5. Because this point provides the model’s only resistance unbal- 
anced roll, its load carefully monitored throughout each test. buildup 
load assures the test engineer that only correct, properly sequenced loads are 
being applied during any load increment. 

Preparation For Testing Model Area.—After the panel point loads for 
condition have been received the test group, approximately four days are 
required ready the model for testing. This time consumed checking 
electrical instrumentation and making provisions load each panel point with 
its correct direction and magnitude load. 

The model mounted the inverted position (Figs. 14, 15, and 16) solely 
reduce the number panel points that must rigged for upload application. 
For virtually all airplane flight situations, most the airplane area has net 
up-loads applied. This would require up-rigging the majority model panel 
points unless the inverted. However, inverting the model, most 
points are simply down-loaded applying weights their load pans. This 
leaves fewer number up-points that must rigged through overhead 
truss and pulley arrangement. 

For all ground handling conditions where the ratio up-to-down loads 
the airplane reverses, all the loads are applied the reverse direction the 
model and the signs resulting strain and deflection data are changed during 
the manual data reduction operation. Because the low stresses maintained 
the model, the effects compressive buckling may disregarded. 

After allrigging and other model setup completed, load card attached 
each panel point load pan indicating its correct incremental load for each 
phase the planned test. These cards are used both the mechanic loading 
the point and the engineer checking the point during the test. 

Preload Model Condition.—In the early phases the program was 
found that retest generally produced better, more linear data than the original 
test. After some study, was decided that this result was caused the ear- 
lier testing having removed the slack “preloaded” the modelin the direction 
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load for the given test. This prevented sudden unexplained jumps 
during the latter test. After number tests were conducted verify this 
assumption, simplified preload was established preliminary the test- 
ing each condition. Whereas the actual condition seldom has loading above 
35% equivalent airplane-magnitude, the preload generally carried only 
25%. This believed provide sufficient “stretch” and leaves margin 
safety that only fewcritical gages are monitored during preload. addi- 
tion, this loading serves indicate possible points trouble such load ca- 
bles binding deflected model structure and others that could cause ex- 
pensive retest unless corrected prior actual testing. 


FIG, 


Test Model Condition.—The test elastic model condition anti- 
climatic when compared with the time involved both pre-test and post-test 
work for condition, requires only one day out the six weeks allocated 
for the complete processing condition. 

Testing begins applying tare load and recording all strain and de- 
flection readings which, where possible, are then set zero. The test con- 
tinues with additional 10% loadings being applied and instrumentation again 
being read each stage. Normally, after attaining 35% load and recording, 
the load removed and again raised 5%. This final tare reading then 
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recorded and compared against the initial reading. Although these tare read- 
ings should theoretically equal and all the intermediate readings perfectly 
linear, much the comments data reduction concerned with why they are 
not. The unloading zero load and return for final readings consist- 
ently provides more linear data than simple reduction from 35% 5%. 


DATA REDUCTION 


The raw strain and deflectiondata measured during test generally include 
five values (5, 15, 25, and 5%) for each significant point. Manual data- 
reduction consists establishing some logical sequence for these values 
that extrapolation 100% may obtained. all cases other than rosette 
data, these extrapolated values are converted airplane magnitude stresses, 
loads, deflections, and are presented form for use stress 
analysis. the case rosette gages, most the data reduction accom- 
plished the use high speed digital computers. 

Computer Programmed Data Reduction.— Most the extensive repetitive 
computations involved elastic model data reduction are coded for solution 
the digital computer. The first these the solution the strains meas- 
ured the 45° rectangular rosettes the model. Using the equations repre- 
sented Mohr’s Circle, these strains are reduced that the complete state 
stress the point obtained. number cases, this procedure pro- 
duces many eleven values stresses point for use stress analy- 
sis the area. 

part this data reduction program provides the final measure the va- 
lidity the test data obtaining comparison internally measured versus 
externally applied loads number “free bodies” taken strategic points 
the model (Fig. 17). performing this comparison, the program uses 
virtually all the strain data from the model for condition and the model 
section properties the region each gage. most cases, these compari- 
sons indicate less than 10% error locations such complete wing isola- 
tion outboard airplane span-station 80.0. practically all cases this com- 
parison internal external loads indicates that the model stresses are 
slightly conservative. Later measurements the B-58 static test program 
have verified this. 

Presentation Model Data.—The entire elastic model program can 
justified only its results have made possible simplified the stress analy- 
sis the B-58 airplane. For this reason, all final model data from each test 
are presented manner that convenient for use the engineer perform- 
ing stress analyses each the airplane components. most cases, stress 
results are shown exact location scaled planforms airplane structure 
(Figs. and 19). This facilitates the stress analyst’s becoming familiar with 
the over-all stress distribution for each condition and enables him inter- 
polate between instrumented points with ease. Extensive data for the wing in- 
clude separate planforms for rosette principal stresses, spar and bulkhead 
flange axial stresses, and spar and bulkhead web shear loads. Similar, less 
extensive planforms are provided for fuselage and vertical tail structure with 
the remaining stress data being tabulated. 

Most the deflection data from the model presented chordwise plots 
vertical deflection number span stations from the airplane center line 
the wing tip. The deflection most points the wing fuselage can 
obtained interpolation within these basic curves. 


a 
oy 
i 
4 
4 
| 
| 
& 


August, 1960 


170 


[\ \ 
\ \ 
\ » 
/ \ 
\ WNYILXI OL TWNUILNI 40 SNOILVLS 
\ \ | (sjaudg 0204195) NIVHLS 
| 
| \ 
| s 
| 
| 
| | | | 


Pal 


F 
| 
= 4 
i. 
f | 
if 
\ 


(SIS6) 


(96 

8 i+ 


+ 


(692¢1) 


(18%) 


x” 


+ 


£292) + 


(0242) 6£2S) | (pees) (gece) 


Hell + 


(ded) 


9082i+ 
(8196) 


40221 + 


(v980!) (8888) 


226! 


9086! + 
(9281!) 


(2212!) 
2+ 


+ 


(826%!) 
OlLe2+ 


(4€98) 
£605! + 


($2191) 
- 


(6162!) 
+ 


882-7 


171 i 
$ 
| | 2 2 
~ 
Z 
| | pre 
© 
2 | 
8 f - | = 
| ++ = + 
in 
| 
| © 
| = | | 
= 
5 
Pa Vs 4 
N ive 
€ 
~ + 
* & 7 
” 
| 


August, 1960 


172 


+ 


9091- 
26% + 


a! 
a! 


| 
' 
I 
4 
3 
> 
q 
q 
= 
n ~ 
2 
+ 
| 
| 
3 


MODEL DESIGN 173 


the end each major phase model testing, composite wing planforms 
are prepared indicate only the most significant value obtained each point 
and the condition responsible. Because the B-58 airplane often experiences 
surface temperatures above 200°F due aerodynamic heating, some consider- 
ation reduced these elevated temperatures neces- 
sary. Because the stress analyst must consider both stress magnitude and 
structural temperature, separate maximum surface stress planforms are pre- 
pared for “standard” and “elevated” temperature conditions. This simplifica- 
tion cannot used understructure, therefore only one maximum web shear 


load planform presented. complete list maximum wing planforms 
follows: 


Rosette Stress Data.— 
Maximum Principal Stress (Hot and Cold Conditions Separately) 
Minimum Principal Stress (Hot and Cold Conditions Separately) 
Maximum Shearing Stress (Hot and Cold Conditions Separately) 

Axial Stress Data.— 
Maximum Spar and Bulkhead Flange Tension Stress (Hot and Cold Conditions 
Separately) 
Maximum Spar and Bulkhead Flange Compression Stress (Hot and Cold Condi- 
tions Separately) 

Web V-sette 
Maximum Spar and Bulkhead Web Shear Loads (Irrespective Sign With Hot 
and Cold Conditions Combined) 


Remaining Direct Model Uses.—The 3/8 scale model now fully occupied 
with determining stress distributions for the B-58 certain final conditions 
that will enhance its use tactical weapon the Air Force. 

result fatigue type failures found number airplanes vary- 
ing types its inventory, the Air Force has placed new emphasis the de- 
termination the service life and fatigue characteristics all airplanes 
scheduled for future procurement. 

This opens entirely new field search for critical conditions because 
they now must evaluated for both load magnitude and number applications 
during the airplane service life. For static considerations, only magnitude 
loading considered. the cyclic fatigue certification program for the B-58 
now progress, the elastic model scheduled tohave major responsibility 
refined stress determination for number conditions. 

addition totesting new conditions, the catalog more than 100 conditions 
already tested will used the less refined stress predictions needed the 
early phases this program. Significant savings are realized since this pre- 
liminary work minimizes the total number conditions tested both the 
model and the full scale cyclic fatigue airplane. 

Related Uses The Model.— Because the 3/8 scale elastic one 
the largest, most extensively instrumented and consistently loaded redundant 
structures ever built, natural source information for the verification 
proposed analytical methods analyzing highly redundant structures. 

Three distinct methods analysis are under study. These are (1) Direct 
Analogy Electric Analog Computer (DAEAC) solution, (2) Consistent Distor- 
tion Solution (Q-07) and (3) Stiffness Matrix Analysis Structures (Q-01). 
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Although all three methods simulated either the B-58 hypothetical struc- 
ture their early stages, they now simulate the 3/8 scale elastic model. 
Where possible, they apply distributions previously tested the model 
and use its stress and deflection data evaluating the accuracy their solu- 
tions. 

Elastic Model Contributions.—Through the use elastic models the de- 
sign the B-58 airplane, the time required for the final qualifications its 
structure has been significantly shortened. Model stresses used the design 
are being closely substantiated the static test program how progress. 
The almost complete absence major structural failures during this static 
testing has validated the merits this program and justified the many man- 
hours involved. This means that inventory airplanes now production will 
not plagued with excessive retrofits and delays the date their combat 
readiness. 


APPENDIX. MODEL SCALING FACTORS 


The derivation complete set scaling factors for loads, material sizes 
and material substitutions must the beginning any model program. 
These factors for given model are dictated such items cost, ease 
fabrication and versatility the finished structure and its data. proper 
adjustment scale factors the early phases model program, certain 
forms output data are assured. most cases, structural models are de- 
signed have their strains equal their prototype. Only this type model will 
considered the following examination. 

Scaling factors will derived three stages. The simple “true scaled 
model,” the adjustments required for dissimilar material substitutions, and 
special scaling factors necessitated desired model deviations will out- 
lines, with special emphasis being shown for those used the 3/8 scale elas- 
tic model. 

True Scaled Models.—A true scaled model one built exact scale with 
direct substitution prototype material and methods for construction. Al- 
though such model design generally impractical, the derivation its scal- 
ing factors demonstrates the basic considerations required model design. 

The following equation expresses the general scaling factor required for 
conversion parameters from prototype model magnitude. 


(Model Magnitude) (Prototype Magnitude) (Model (1) 


which (n) represents the number times the model scale factor (S) will 
involved the particular parameter. For example, the area (A) involves the 
square linear dimension; therefore, the case loading (P) 
which involves area maintaining equal strains, the case moment 
inertia, involving both area and the square linear distance, 

Using these basic rules, the following geometric relationships between the 
model (subscript and its prototype (subscript can stated. 


(2b) 
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Since the area (A) function width (b) times thickness (t), 
(3a) 


Similarly the moment inertia (I) and polar moment inertia (J) are 
functions area times the square distance, therefore, 


but 


Applied loads (P) must adjusted maintain equal strains between the 
model and its prototype. Starting with this equality strains, 


and 


Because moment (M) and torsion (T) involve load times distance, 


These equations are equally true for external internal loads. 
The scale factor for deflection (5) equivalence may verified using can- 
tilever beam length and constant with vertical load its free end. 
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(16b) 


the same manner can shown that torsion applied the free end 
yields, 


Tm lm 
and 
(18b) 


Scaled Models Dissimilar Materials.—In many cases during the develop- 
scaled model program becomes apparent some material 
other than that used the prototype either convenient absolutely neces- 
sary. This concept easily incorporated applying certain adjustments 
the true scale model factors. 

Adjustments will only required where scale factor involves consid- 
eration thickness for either web surface material. Geometric dimen- 
sions such width, length and depth will not affected. The general equa- 
tion for equivalent thickness substitutions for dissimilar model materials is, 


Because the equations expressing area (A), moment inertia (I), and polar 
moment inertia (J) are functions material thickness, their model magni- 
tude may obtained multiplying their true scaled model values the ratio 
the modulus elasticity the prototype material that the model, or; 
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The consistent application these scale factors yields model and proto- 
type with equal strains. convenient reducing strain data determine 


prototype stress directly using only model strain and prototype moduli elas- 
ticity rigidity. 


sirable incorporate special features that are apparently inconsistent with 


TABLE 1.—MODEL CONVERSION FACTORS FOR SIMILAR MATERIALS 
THE PROTOTYPE AND MODEL 


Subject 
Indepen- 3/8 Model 
Conversion dent Conversion 
Subject (S)n Magnitude Factor Factor 


Geometry 


Width 
Length 

Depth 

Material Thick- 
ness 


Material Properties 


Area 

Moment Inertia 

Polar Moment 
Inertia 


S2A 
s4 A d2 


Load 
Moment 
Surface Strain 
Deflection 
Twist 


Apply “K” factor correction the true scale model conversion factor, 


The factor vanishes from these equations because appears once both the 
numerator and denominator. 


both its prototype and true scaled model. Two these special features used 
the 3/8 scale elastic model design were the substitution flat webs the 
model for comparable corrugated prototype webs and the substitution equal 
thickness magnesium skins and webs for prototype aluminum structures (Fig. 
11). both cases, valid model data canbe obtained proper adjustments are 
made scaling factors for both the model’s structure and external loads. 

The true scaled simulation rather thick prototype sandwich panel with 
single flat skin the model would yield model surface which would obvi- 
ously much more criticalin compressive buckling than the airplane. How- 
ever, the combined thickness prototype aluminum skins simulated with 
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equal thickness model magnesium skin, standard sheet stock may used 
and the buckling problem decreased. All that remains are adjustments 
the affected scaling factors. 

The substitution magnesium flat webs the model for equal thickness 
corrugated aluminum webs the prototype requires scale factor adjust- 
ment other than the “growth factor” reduction explained the paper. How- 
ever, the original substitution magnesium skins and webs directly for pro- 
totype aluminum counterparts requires the following scaling factor adjustments. 

adjustment (K) must derived for the external loading (P) make the second 
equation valid. 


(22) 
(23) 
(24) 
Therefore, 


the specific case the 3/8 scale model where magnesium versus alumi- 

num substitutions are made, this becomes 
6.5 

Because the “K” factor function the model load (P) all items that are 
function this load will require the same correction. 

The scaling factors for the 3/8 scale model and true scale model, using 
similar materials, are listed Table this Appendix. 


Model Scale Ratio 
Distance (inches) 


Material Thickness (inches) 


“ 


Material Area 

Moment Inertia 

Polar Moment Inertia 
Modulus Elasticity (pounds/inch2) 
Modulus Rigidity (pounds/inch2) 


F 
are 


MODEL DESIGN 
Load (pounds) 


Moment (inch-pounds) 


Torsion (inch-pounds) 

Axial Strain (inches/inch) 
Deflection (inches) 

Twist (radians) 

Subscript, pertaining prototype 


Subscript, pertaining model 
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SHRINKAGE, SWELLING AND CREEP 


Closure Hrennikoff 


HRENNIKOFF, ASCE.—The writer wishes acknowledge the private 
communications Mr. Powers, Research Counsellor the Portland 
Cement Association, the subject this paper. 

Mr. Powers calls attention the author’s statement, “The presence film 
(adsorbed) water cement was demonstrated several other investigators 
without, however, attributing the mechanical action (shrinkage-swelling) 
brought this paper.” This statement incorrect since Mr. 
clearly designated the adsorbed water cause shrinkage-swelling phenom- 
enon, 

Although recognizing the Mr. Powersin realizing the role played 
the adsorbed water, the writer considers only fair point out that Mr. 
Powers’ article, unlike the present paper, contains detailed explanation 
the mechanical action film water. Furthermore, alongside the adsorbed 
water theory, Mr. Powers devotes comparable space the exposition the 
capillary action another cause shrinkage, assertion with which the 
writer, explains some length his paper, strongly disagrees. 

Mr. Powers also states that according the recent observations with the 
electron microscope, the units cement gel are not spherical, reported 
Bogue, but have the shape crumpled sheets and ribbons. 

Leaving the final word this subject the specialists electron micro- 
photography, may said that the shape the gel unit not vital the in- 
tegrity the author’s theory, since the spherical shape vouched for the time 
acompetent authority, was used the author only for the sake illustration. 
The significant remaining fact, irrespective the gel elements being round, 
flat elongated, that they are distinct units surrounded water, and the 
author’s reasoning, leading the conclusion shrinkage-swelling and result- 
ing from the attraction the adsorbed water, way affected their 
shape. Also unaltered the conclusion that shrinkage produced drying 
much smaller than the accompanying decrease thickness water films. 
The disparity between the loss water and the extent shrinkage dis- 
cussed also both Messrs. Neville and Keith Jones. Mr. Neville pre- 
sents some figures referring mortar cylinders illustrating this point. 

correctly states that the rigid framework the particles sand and un- 
hydrated cement clinker accounts for the smallness shrinkage compar- 
ison with the loss water. However, the writer wishes add this the point 


Civ. Engrg., Univ. British Columbia, Vancouver, Canada, 
Bulletin 22, Portland Cement Assn., 582, 


“Chemistry Portland Cement,” Bogue, 2nd Edition, Reinhold Publishing 
Corp., 1955, 
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made him the basis his theory. Every time the water film around the 
elements cement gel decreases, the attractionbetween the latter brings them 
closer together, but the extent this mutual approach, the nature the laws 
attraction betweenthe water andcement, muchsmaller than the lost thick- 
ness the water film. 

Mr. Jones suggests the unequal size gel units possible cause for the 
disparity between the amount shrinkage and the corresponding water loss. 
visualizes the gel units belonging two categories; the ones contact 
with eachother forming framework giving rigidity the whole mass, and the 
others occupying, more less loosely, the empty spaces between the units 
the first kind. The above mentioned information supplied Mr. Powers with 
regard the shape the gel units suggests irregularity even more far 
reaching than the one intimated Mr. Jones. diversity conditions with 
respect the available moisture the grain contacts then becomes strong 
probability. the same time, the writer doubtful the presence float- 
ing and loosely connected grains the mass hydrated cement. feels that 
the units cement gel either form situ out the grains unhydrated ce- 
ment, they form precipitation from the water solution. both cases the 
units gel, being under the action strong chemical and physical forces, 
would form ina close contact with their neighbours and not ina loose aggrega- 

The writer does not know what inconsistency between the author’s theory 
and his experimental results Mr. Jones refers his discussion. means 
the disparity between the loss water and the resultant shrinkage, then the 
writer considers this point not inconsistency but strong factor 
favor his theory. 

The writer grateful all discussers for their constructive contributions 
the subject the paper. regrets that the private communications Mr. 
Powers had not been made part the public discussion, but considers 
inappropriate present here both Mr. Powers’ and his views number 
interesting points brought their correspondence. 
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COMPARATIVE STUDY SEGMENTAL ARCH 


VELTROP®, ASCE.—Mr. Zienkiewicz has presented interesting 
and valuable discussion three the most frequently used systems shell 
equations for arch-dam design comparing moments, thrusts, shears, and de- 
flections along the arches for four values the central angle: 90°, 120°, 180°, 
and 360°. The last two values are little interest arch-dam designers, but 
the last one serves the purpose illustrating the limitations Donnell’s equa- 
tions. 

The comparison with the classical theory and especially with the corrected 
classical theory most interesting, particularly the explanation the phys- 
ical meaning the differences between Timoshenko’s and Flugge’s equations. 

Since the purpose the paper primarily compare various shell theo- 
ries, the author has limited himself constant radius and constant thickness 
arches. Numerical values are further limited one set arch dimensions 
with R/t 

Mr. Zienkiewicz has also limited himself arches with fixed abutments and 
with Poisson’s ratio equal zero. should noted that the author’s conclu- 
sions are not limited the specific numerical cases presented, but hold, gen- 
erally, the common range R/t ratios and central angles. 

The present discussion expands the original paper somewhat and adds few 
observations which should particular interest designers. The 
computations covered this discussion consist of: 


Calculations for moments, thrusts, and deflections abutment and crown 
fora number R/t ratios: 0.025, 0.050, 0.075, 0.100, 0.150, and 0.200. These 
are limited toa 120°, which representative for arch-dam de- 
sign. These computations are performed first all for three the methods 
discussed the author, namely Donnell, classical, and classical corrected 
addition, the following methods have been added: classical 
and classical corrected with 0.2, data from USBR tables which include abut- 
ment deformations, data from Parme’s curves, and Donnell’s equation with 
abutment deformations. 


Like Mr. Zienkiewicz, the writer has limited these computations constant 
thickness circular arches, otherwise the amount work would become pro- 
hibitive without adding the basic purpose the paper the discussion. 
For clarity the ensuring discussion the writer has numbered and compiled 
the various methods and their limitations, far they are relevent this 
discussion, Table Needless say, all these methods include, course, 


Head, Arch Dam Sect., Harza Engrg. Co., Chicago, 
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the so-called ribshortening effect and the usual assumption the theory elas- 
ticity. 

The results the writer’s computations for methods through are pre- 
sented Table together with the additional data computed for the author’s 
method The author’s and writer’s results both cover the case half central 
angle 60° These values aretherefore tabulated separately 
Table for direct comparison all ten methods computation. 

All formulas and numerical results have pressure applied 
the middle plane the arch. 

Explanation order comparison between the var- 
ious methods used the writer and the classical equation presented Mr. 
Zienkiewicz, the various equations are written down and assumptions briefly 
discussed illustrate pertinent differences. 

Methods and 6.—The classical theory arch rings includes deformations 
due axial forces, moments, and shears, while abutments are considered 


TABLE 2,—COMPARISON METHODS 


Method 


Abutment 
Deformation 


Shear De- 
flection 


Effect 
Axial Force 
Rotation 


Poisson’s 
Ratio 


Donnell 
Timoshenko 
Fligge Yes 
Classical 
Classical Corrected Yes 
Classical 0.2 
Classical Corrected Yes 


rigid. Using standard methods, one arrives the following expressions (using 
the author’s notation, for ease comparison): 


crown 
and 


abutment 


| 
ak 
sin 
a 


Donnel 
Class, 
Class, corr, 
Class, 
Class, 
USBR Tables 
Parme’s Curves 
Donnell Corr. 


Donnell 

Class, 
corr, 
Class, 
Class, corr, 

USBR Tables 
Parme’s Curves 
Donnell Corr, 


Donnell 

Class, 

corr, 
Class, 

Class, corr, 
USBR Tables 

Parme’s Curves 
10, Donnell corr. 


Donnell 

Class 

Class, corr, 
Class, 

Class, corr, 
USBR Tables 

Parme’s Curves 
10, Donnell corr, 


Donnell 

Class, 

Class, corr, 
Class, 

Class, corr, 
USBR Tables 

Parme’s Curves 
10, Donnell corr, 
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.0006690 


.002432 
.002282 


.0007116 


Moment 


.001414 
.001278 
-00131 


.0003252 
.0003239 


.005377 


Crown 


.003201 
.002845 
.002822 
-00291 


Thrust Abutment 


Thrust Crown 


1,8909 
1,8908 


.011058 
.009330 
.010 
.008300 

.011111 


.005529 


1,8410 
1.84 

1,8270 


.01051 


.9344 
.9730 
-976 


.9399 
-939 


¢€ . 
t 
co 


-0167 


1,6670 
1,7084 
1,7001 
1.71 


TABLE 
1 - | - - ,02417|- .04090 
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with 
and where 


first approximation. 
Eqs. and are identical Eqs. and 9b, when made equal zero. 
The contribution shear included the coefficient with the term E/G. 


TABLE 


Method 


Parme’s Curves -0.931 0.00494 1.84 
Donnell Corr, 0,00556 1,8270 


From the theory elasticity have E/G 2(1 that the following nu- 
merical values are possible: 


writer: 0.2, k=1.25 3.0 (method 


The use the values 2.5 and 3.0, therefore, illustrates the effect Poisson’s 
ratio. should noted that the author’s value equal the writer. 

Methods and 7.—The classical corrected method includes all effects 
classical methods and but corrected for rotation due axial forces and 
The equations for moments, thrusts, and deflection are: 

crown 


| 
a 


and 


abutment 
wi 


where again Eq. applicable. 

Several commentsare inorder. First, should noted that Eq. for con- 
stant derived directly from Eq. for constant going through 
the classical derivations for both the energy method can that these 
differences are the direct result correction term for the work done the 
axial force rotation. However, expressions for moments and thrusts (Eqs. 
11, 12, 14, and 15) can derived directly from Eqs. 18, 19, 21, and re- 
placing h/R with zero. actual numerical calculations appears that the 


term always very close however, the term h/R cannot neg- 


lected the end Eqs. and 21. This verifies the author’s statement that 
the main effect the correction the classical method just the additional 
term h/R. 


Comparison crown deflection (Eqs. and 20) shows that Eq. has 
extra term besides few multipliers The numericalimportance will 


elaborated upon subsequently discussing the results. 

Method 8.—The method used the USBR (25) (26) essentially the classi- 
cal method, without the “correction” methods and (25a). The impor- 
tant difference with methods and inclusion abutment deformations 
the USBR method. The resulting equations are rather cumbersome and are not 
therefore worked out into simple expression. 

Method 9.—Parme’s curves are based the classical method for 0.25, 
and for rigid abutments. Details the derivation and equations for moments, 
thrusts, and crown deflection can found (27). 

Method 10.—Donnell’s equations were corrected for abutment deformations, 
since these expressions are somewhat easier handle than Timoshenko’s and 
The writer derived moment, thrust, and deflection equations crown 
and abutment including the effect abutment rotations. This effect was intro- 
using the following boundary conditions: 


tained (25b). 
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The equations for moment, thrust, and crown deflectionare the sameas Eqs. 
through 7d, except that the coefficient changes 


with 


Discussion Results.— convenience the subsequent materialis divided 
into two categories. The first (A) compares moments, thrusts, and deflections 
due different methods (No. through and 9), while the second (B) com- 
pares the effects abutment deformation and axial work rotation again 
for moments, thrusts, and deflections. 

Comparison Methods. 

already indicated the author his conclusion (a), the deflections 
computed from the three shell theories show very good agreement. Table 
shows that this also holds true for those computed the classical methods 
(both with and 0.2). The small differences between Donnell’s re- 
sults and the classical methods increase for thicker arches, shown Table 
This due the effect shear deflection which included the classi- 
cal methods but not Donnell’s equations. 

Substantially the same can said about the thrust abutments and 
crown, again Donnell’s values are somewhat smaller. Moments show rela- 
tively larger differences than thrusts and deflections because inclusion 
the higher order terms Flugge’s equations and the work effect axial force 
rotation the classical corrected method, further discussed under 

The differences between values obtained from Parme’s curves (method 
and methods through are generally somewhat larger than differences 
among methods through This may due the inaccuracies involved 
the rather rough graphical interpolations from Parme’s curves. 

Conclusions. 

generally give satisfactory results but lead substantial 
errors, which increase magnitude when applied arches, with t/R 0.1, 
approximately, indicated Donnell’s values Table 

For the range variables covered this analysis there practically 
effect due Poisson’s ratio. 

The suggestion the author, made his conclusion (d), use the sim- 
ple crown deflection formula derived from Donnell, useful one, but the 
error involved may larger than suggests, can seen from Table 
comparison with method 

Comparison Basic Assumptions. 

The greatest differences moments, and deflections but not thrusts 
are found between the USBR method which includes abutment deformations and 
the other methods, which not include this effect. This follows directly 
comparing the classical method with the USBR method the USBR being 
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based exactly the same derivation throughout, except for the abutment de- 
formation. This deformation decreases abutment moments considerably while 
hardly affecting crown moments. From Table can seen that changes 
abutment moments appear more consistent than those crown moments, 
even though numerical values were computed with two more decimals than 
shown for methods and Perhaps the limited number digits pre- 
sented the USBR tables accounts for the lack consistency data for 
method Crown deflections also increase substantially. 

This difference occurs for all values t/R, but can seen from Table 
that the significance these differences increases for thicker arches. Table 
also illustrates the large differences crown deflection between the USBR 
method and all other methods which overshadow the relatively small differ- 
ences among the others. This leads the conclusion that inclusion abut- 
ment deformations inthe archcomputation important for arch-dam analysis 
especially when cantilever deflections are equated with those the 
arches. 

The effect abutment rotation was also added method the results 
are shown method 10. comparing Eq. with the expression the author 
for the case fixed abutments canbe seen that Eq. 
should normally lead smaller values and, therefore, larger values 
deflection, thrust, and moment, can seen also from Table comparing 
results method with those method 10. This obviously contrary ex- 
pectations, surely abutment moments should decrease, confirmed the 
classical methods and USBR data. The error seems inherent Donnell’s 
equations for moment and thrust (see author’s Eq. and 7d, respectively). 
distinction made there between abutment and crown thrust, and, further- 
more, abutment moments are always twice crown moments, regardless 
dimensions. Both the expressions for moment and thrust are apparently too 
oversimplified reflect such refinements abutment deformations. Conse- 
quently inclusion abutment deformation does not lead the correct effects. 

Inclusion the work axial forces rotation inthe classical corrected 
method has substantial effect moments, but not thrusts and crown de- 
flection. One would expect this the correction does not, first instance 
least, alter the static equilibrium, and that what determines the order 
magnitude the thrust. (Note, for example, that for simple ring stresses 
from the cylinder formula, when ribshortening due thrust neglected, the 
coefficient for thrust crown and abutment would 1.00). Then also 
follows that the geometric pattern radial deflections not basically changed 
illustrated the small differences coefficients for clas- 
Sical and classical corrected methods. 

One can see directly that with the effect axial forces rotation the mo- 
ment the crown needed guarantee continuity now smaller, the coef- 
ficient shows. Subsequently, static equilibrium half the arch will show that 
the abutment moment must then increase. 

compare the effects abutment deformation with those the work 
axial forces rotation reference made tothe coefficients Table for 
methods and For thrusts and deflections inclusion abutment 
deformations significant, whereas the “correction” for work not. For 
moments the results are quite different. Abutment deformation considerably 
decreases abutment moments, but the “correction” causes almost equal 
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increase. This effect would seem contradict the USBR statement (25a) and 
needs further investigation. Abutment deformation changes the crown moment 
only slightly, but the “correction” causes considerable decrease. 

Conclusions. 

USBR data show that inclusion abutment deformations arch compu- 
tations important, especially for crown deflections and abutment moments. 

has been shown that Donnell’s equations are not accurate enough 
properly reflect the effect abutment deformations. 

The inclusion work axial forces rotation, has significant effect 
moments abutment and crown, but effects thrust and deflection are 
negligible. This conclusion appears contradict USBR practice neglecting 
this effect. 

writer would like thank the several engineers 
the arch dam section who have checked the numerical calculations, especially 
Wengler who checked all derivations and reviewed the entire discussion. 
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WHITE NOISE REPRESENTATION 


describes interesting electrical experiment which artificial earthquake 
produced random (white noise) function generator used excite struc- 
tures. Fig. the author compares the results his experiment with av- 
erage earthquake spectra and shows that there close agreement except 
the low period end where the earthquake spectra have characteristic hump. 
The largest deviation for the zero-damped spectrum curve, and here the 
author’s points are, presumably, not too reliable the zero-damped case 
particularly difficult handle accurately and larger samples would need 
averaged give accurate results. 

interest compare the author’s spectrum points with the results 
simple approximate analysis based energy balance 
the excitation oscillator has random characteristic, then, the aver- 
age, the the system increases until the average rate energy 
input equal the average rate energy dissipation. For harmonic motion 
single degree freedom viscous damped system, the average rate 
energy loss and, therefore, the average power input 


The solution this equation 


ten terms the relative velocity spectrum V(A, 


Combining Figs. and yields 


California Inst. Tech., Pasadena, Calif, 

Engrg., Calif. Inst. Tech., Pasadena, Calif, 

“Response Spectrum Techniques Engineering Seismology,” Hudson, 
Proceedings, World Conf, Earthquake Engrg., San Francisco, Calif., June, 
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Pty 

and 
choosing such way that the best representation over the whole range 
parameters obtained. Eq. 12, the time which the average 
computed, however, the author’s case the value always chosen 
give the largest value V(A, T). The time will, thus, less than the 
sec duration the input. The value will also influenced the fact 
that the input limited finite band width. 

putting 6.5 sec Eq. theoretical spectrum curves are obtained 
shown Fig. 17, which aiso shows the author’s data points. The agreement 
this simple analysis with the author’s tests very good. All the author’s 
points seem bit high the low period end and this presumably indica- 
tion that the input was not true white noise. would, course, possible 
accentuate this effect the use appropriate filters modify the input 
that produces spectra virtually identical with earthquake spectra. 


Eq. can fitted the scale used Fig. setting 


FT/SEC 


% 02 04 O06 08 1.0 1.2 1.4 1.6 18 20 22 24 26 28 3.0 
PERIOD (SECONDS) 


FIG, 17.—THEORETICAL AND EXPERIMENTAL SPECTRA 


actual earthquake record not Gaussian and, course, the author’s 
input could not truly Gaussian. would interest see some 
the frequency distribution the experimentally determined values. 
the range these V(A, values were inconsistent with the range for actual 
earthquakes the input could modified correct for this. Another possibility 
would construct artificial earthquake whose average spectra 
agrees exactly with those real earthquakes. Earthquake studies could then 
made along the same lines that random shaking tests are used for missile 
reliability test studies. 

The maximum yield displacements computed Mr. Bycroft appear quite 
reasonable. may that some the differences between the white noise 


“Properties Strong Ground Motion Earthquakes,” Housner, Bulletin, 
Seismol, Soc. Amer., Vol. 44, July, 
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yield displacements and those computed the previous paper are accounted 
for inaccuracies the earlier method computation. 

difficult know exactly what interpretation place the computed 
maximum yield displacements since ideal-plastic force-displacement re- 
lation was assumed for the oscillator. made for elasto- 
plastic also showed sizeable permanent displacements. actual 
building, course, the force-displacement diagram would not have flat plas- 
tic yielding but the diagram would curved with slope greater than zero and 
result the maximum yield displacements would presumable smaller. 

Studies along the lines the author’s work are much needed, because one 
the key questions regarding the behavior structures during earthquakes 
how much energy can absorbed during the vibration without having the 
structure collapse. 


“The Analysis Structural Response Earthquakes,” Berg, 291, 
Univ., Michigan, May, 1958. 
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STRENGTH AND EFFICIENCY ASPECTS PLATE 


Discussion George Pincus and George Winter 


GEORGE A.M. ASCE and GEORGE ASCE.—The 
first part Mr. Gerard’s paper re-formulates well established relationships 
the post-buckling strength compressed plates. The supporting evidence 
consists numerous tests the field largely non-ferrous aeronautical 
structures. gratifying note that these relations are strikingly close 
agreement with similar information developed some fifteen years ago one 
the present writers from tests light-gage steel structural members 
the civil engineering type.17,18,19 The design which has been force 
nationally since 1946 for this widespread type construction, incorporates 
this information explicit design formulation.22 For this reason may 


some interest compare and correlate Mr. Gerard’s presentation with this 
earlier information. 


Th. von Karmanwas the first23 present semi-intuitive equation for the 
post-buckling strength plates, formulated interms ofthe effective equiva- 
lent width rather than terms Mr. Gerard’s average stress failure, 
For rectangular plates with simply supported edges, this equation reads 


April, 1960, George Gerard, 
Graduate Asst., School Civ. Engrg., Cornell Univ., Ithaca, 

and Head, Dept. Structural Engrg., Cornell Univ., Ithaca, 

“Strength Thin Steel Compression Flange,” George Winter, Transactions, 
ASCE, 112, 527, 1947, reprinted with added appendix Cornell University 
Engrg. Experiment Sta, Reprint No, 32, 

“Performance Thin Steel Compression Flanges,” George Winter, Third Con- 
gress Int. for Bridge and Structural Engrg., Liege, Publ., 317, 1948; 
reprinted Four Papers the Performance Thin Walled Steel Structures, Cornell 
Univ, Engrg. Experiment Sta. Reprint No, 33, 1950. 

“Performance Compression Plates Parts Structural Members,” George 
Winter, Research Structures Suppl., 49, London, 1949, 

“Specification for the Design Light Gage Cold-Formed Steel Structural Mem- 
bers,” Iron and Steel Inst., first issued 1946, latest edition 1960. 

“Cold-Formed, Light Gage Steel Construction,” George Winter, Proceedings, 
ASCE, Vol. 85, No. 151, 
“Commentary the 1956 Edition the Light Gage Cold-Formed Steel Design 
Manual,” George Winter, Amer. Iron and Steel 1958, 

“The Strength Thin Plates Compression,” Th. van Karman, Sech- 
ler, and Donnell, Transactions, ASME, 54, 1932, 53, 
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Using Eq. with this can re-written 


this can also written 


which identical with Mr. Gerard’s Eq. with 1.0 0.5. 
its alternate form the same relation can also written (see Eq. 37) 


which identical with Mr. Gerard’s Eq. with 1.9, 

Sechler24 was the first note that the constant 1.9 von Karman’s 
equation needs adjusting match test results. Further explorations post- 
buckling strength were carried out during World War II. that time security 
restrictions resulted independent investigations aeronautical engineering 


the one hand, and civil engineering the other. the basis his 
tests, gave approximate lower bound, that for simply sup- 


ported edges, 
1/2 1/2 


which, include other edge conditions, can written more generally 


b — Se . Je - 


seen that the terms the square bracketts Eqs. and represent 
the experimental corrections von Karman’s original Eqs. and 41. 

and are valid not only failure, but for any edge stress lar- 
ger than that causing initial buckling. The collapse load was found 
reached when this edge stress attained the yield point, finding which 
Mr.Gerard states tobe agreement with the later analyses Stowell, Mayer 
and Budiansky. 

important note that Eqs. and are formulated terms the 
effective width rather than the average collapse stress This permits 
conventional strength materials methods tobe applied tothe design mem- 
bers (such beams) which the particular compression plate merely one 
several components (such the compression flange). this case, the usual 
section properties apply, provided they are based the effective, rather than 
the actual width. The fact that the effective width concept applies not only 


“The Ultimate Strength Thin Flat Sheet Compression,” Sechler, No. 
27, Guggenheim Aeron, Lab., Calif. Inst. Tech., 
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Since, the definition the effective width, 
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the edge stress but also edge stresses smaller value, 
provides more general tool and design. permits one analyze 
not only the case when the limiting load the member governed the col- 
lapse the compression plate, but also the frequent situations when de- 
termined prior yielding the tension flange, excessive deflection 
design load.22 For these reasons believed that the use effective 
width preferable that average stress collapse the com- 
pressed plate. 

interesting compare numerically, the results Eq. with those 
Eq. (Fig. 10). The former was obtained primarily from tests individual 
plates non-ferrous metals while the latter was developed from tests en- 
tire thin-walled steel members, which the pertinent plate was merely one 
component. Correspondingly, these latter tests some rotational edge re- 
straint was always present, its amount depending the configuration the 
section. However, many sections this restraint was small, and since Eq. 

1.0 


FIG, COLLAPSE CONDITIONS FOR SIMPLY SUPPORTED 
COMPRESSED STEEL 


was fitted closely the lower bound the scattering range, represents 
substantially the conditions simply supported edges. 

For these conditions Mr. Gerard’s Eq. with 0.85 and 1.42 from 
Eq. and using for mild structural steel ksi and 30,000 ksi, be- 
comes 

0.85 


the other hand, for the same material’s constants, Eq. becomes 


These two expressions have been plotted Fig. 10. 


; 
4 
0.4 
Eq. 43 
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seen that very close agreement obtained over the entire relevant 
range. slightly more conservative, which consistent with its char- 
acter lower bound expression. Only the region low b/t-ratios, when 
the plate fails average stress excess about 85% the yield strength, 
the two expression deviate more strongly. this range, however, buckling 
occurs plastic domain and the problem post-buckling strength becomes 
largely meaningless. Also, performance this range depends strongly the 
detailed inelastic stress-strain characteristics the particular material 
that single expression cannot expected apply equally well all situ- 
ations. 

view the close agreement for this special case, there little reason 
doubt that would exist betweenthe more general expres- 
sions, Eq. and the writers’ Eq. 43. 

summarize, gratifying note that the information post-buckling 
strength compressed plates, which for many years has been widespread 
use light-gage steel structures the civil engineering type, very satis- 
factorily confirmed Mr. Gerard’s recent data, obtained the field non- 
ferrous aeronautical structures. 
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129. the beginning the second full paragraph, change “The decrease” 
“The decreases” 

131, line 17. Change 0.55R 0.055R 

134, line Change Gto A+C 
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DIVISION ACTIVITIES 
ENGINEERING MECHANICS DIVISION 


Proceedings the American Society Civil Engineers 


NEWS 


August, 1960 


UNDERGRADUATE PROGRAMS ENGINEERING MECHANICS 


Dr. Pletta, Chairman Engineering Mechanics Division ASCE and 
Head Dept. Engr. Mech. V.P.I., submitted the following data the 
growth undergraduate engineering mechanics programs the United 
States. The first was established the Mich. 1928. was followed 
another Purdue about 1940 which 1954 was renamed engineering sci- 
ence. 1957, one was established and 1959, another Columbia. 
All these are now approved the ECPD. addition Michigan, North- 
western and Penn State were awarded accreditation engineering science 
1959. engineering physics, Cornell (1951), Harvard (1951), Kansas (1949), 
Maine (1949), Michigan (1954), Oklahoma (1953), Toledo (1956), and the 
Naval Post Graduate School (1955) were accredited ECPD. The growth 
undergraduate curricula these disciplines has been quite healthy, especially 
recently, and should the future have considerable impact upon our division. 
There are the moment other schools such Manhattan College, University 
Wisconsin, that are just starting programs but these are not yet eligible 
for accreditation because the fact that they have not yet graduated any 
students. 

The editor has noted that many regard the development these separate 
programs sound evolution engineering education, while others regard 
the apparent trend with apprehension. The second group prefers regard 
engineering science and mechanics important part the spectrum be- 
longing the professional depts. rather than separate discipline crowd- 
between the professions and the sciences. Perhaps Purdue’s recent action 
combine the Engr. Science Dept. with professional dept. (Aero) reflects 
this philosophy. 

Concise reader comments are welcome and will considered for presen- 
tation here. 


Note.—No. 1960-24 part the copyrighted Journal the Engineering Mechanics 
Division, Proceedings the American Society Civil Engineers, Vol. 86, No. 
August, 1960. 


Copyright 1960 the American Society Civil Engineers. 
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SECOND CONFERENCE ELECTRONIC COMPUTATION 


was announced previously, this conference will held September 
and 1960, the Pittsburgh Hilton Hotel. planned mail information 


the program 18,000 ASCE members. list the papers and authors 
follows: 


Thursday morning, Sept. 


Impact Electronic Computation Civil Engineers (Keynote Address), 
Richardson. 

Computer Program Exchange: Myth and Reality, Chang. 
Computer Design Structural Steel for Buildings, Zar and Beck. 
Computer Solutions Linear Buckling Problems, Sylvester. 
Error Analysis for Eigenvalue Problems, McCalley, Jr. 

Optimum Design Transmission Towers, Anaston. 

Optimum Design Reinforced Concrete Buildings, Graham. 

Vision Our Automatic Future (Luncheon Address), Dean. 


Thursday afternoon, Sept. 


Computer Analysis Structures, Borges. 
Structural Idealization for Digital Computer Analysis, Archer and 
Samson, Jr. 
Structural Design Systematic Synthesis, Schmit, Jr. 
Simultaneous Equations Solved Over-Relaxation, Lehman. 
The Electronic Computer Tool Moment Distribution, Soehrens. 
Multi-Story Frame Analysis Digital Computer, Rubinstein. 
Stiffness Method Rigid Frame Analysis, Pei. 
Computing Maximums Due Moving Loads, Lingeman. 
Some Fundamental Concepts and Their Use Matrix Structural Analysis, 
Berman. 
The Elements Matrix Structural Analysis, Shore. 
Dynamic Analysis Circular Arches, Eppink, and Veletsos. 
Numerical Procedure for the Analysis Continuous Plates, 


Ang, and Newmark. 


Friday morning, Sept. 


Sloping Surcharge Retaining Wall Design, Wadsworth. 
Vessel Foundation Designs Using Digital Computer, Horton, and 
Eichmann. 
Stress Distribution Patterns and Settlement Characteristics Structural 
Pile Foundations, Rausa. 
Numerical Analysis Laterally Loaded Piles, Reese, and Mat- 
lock. 
Computational Technique for Three-Dimensional Pin-Jointed Structures, 
Pearson. 
The Finite Element Method Plane Stress Analysis, Clough. 
Numerical Analysis Applied Beam Vibrations, Wilson. 
Inversion Band Matrices, Asplund. 
The Use Computers Election Forecasting (Luncheon Address), 
Woodbury. 
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Friday afternoon, Sept. 


Computer Analysis Twin Box Culverts, Hoffman. 

Computer Analysis Pier with Unlimited Shape, Barkocy. 

Computer Design Prestressed Concrete Beams for Deflection Control, 
Kulka. 

Computer Methods for Dynamic Structural Response, Neubert. 

Matrix Analysis Non-Linear Structures, Wilson. 

New Methods Matrix Structural Analysis, Klein, and Chirico. 


NATIONAL AMERICAN WELDING SOCIETY MEETING 


According preliminary program printed the June 1960 issue 
WELDING JOURNAL, ASCE will participate the sponsorship four ses- 
sions the AWS Fall Meeting held September 26-29, the Hotel Penn- 
Sheraton, Pittsburgh, Pa. Titles and authors for these sessions are fol- 
lows: 


Monday afternoon, Sept. 


Welding High-strength Low-alloy Structural Steels for Bridges, 
Beaton, and Jonas. 

Behavior Welded Built-up Beams Under Repeated Loads, Stall- 
meyer, Hall, and Munse. 

The Flexural Fatigue Strength Stiffened Beams, Stallmeyer, and 
Munse. 


Tuesday morning, Sept. 


Tests Welded Plate Girders, Yen. 

Design Columns Multi-story Frames, Ellis. 

Welded Interior Beam-to-column Connections, Graham, Sher- 
bourne. Khabbaz, and Jensen. 


Tuesday afternoon, Sept. 


Design Features for Welded Joints Space Trusses for Retractable-roof 
Auditorium, Cohen. 

Plastic Design Fixed-base, Gabled Frames, Prawel, Jr. and 
Ketter. 

Driscoll. 


Wednesday afternoon, Sept. 


Strength Welded Box Beams, Brungraber. 
Research and Development Continuous Welded Rail, Magee. 
Painting Welds, Keane. 


planned publish more detailed program the September 1960 
issue WELDING JOURNAL. 
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WEST COAST COMMITTEE MINUTES 


The following are excerpts from committee meeting held April 1960, 
the University California, Berkeley: 


The ASME Applied Mechanics Executive Committee has approved 
principle that the national summer conference Applied Mechanics 
should held west the Rockies every fourth year and those years 
would amalgamated with the West Coast Conference. The ASCE 
committee was complete agreement with this proposal and expressed 
the desire still participant these particular meetings, although 
they would primarily run the ASME National Committee Ap- 
plied Mechanics. 

Because larger number papers, the ASME National committee 
like schedule more their papers for presentation the 
West Coast Conference. This means either longer meetings with more 
sessions, simultaneous sessions, fewer ASCE papers. 


This committee feels there objection simultaneous sections. 
(There were several the 1959 West Coast conference.) There are 
definite advantages retaining ASCE participation this conference. 
particular, this committee felt there are, least some years, al- 
ready too many separate conferences and that applied mechanics 
binding area joint,interest both ASCE and ASME members. 

The 1962 National Congress Applied Mechanics scheduled 
for Berkeley 1962. The Committee agreed that probably desir 
able not have separate 1962 West Coast Conference. 

desirable that the West Coast Committees both ASCE and ASME 
have representative the committee from the school area which 
sponsoring meeting. 

ASCE participation the West Coast conference should continue 
grow. The West Coast Committee should request budget for these 
conferences help defray expenses, previously carried solely 
ASME. This could used for printing and mailing costs, losses 
meals, Honorariums for General Lecturers, etc. 

ASCE should investigate the possibility furnishing preprints ASCE 
papers the meeting. There was also some discussion length 
ASCE papers Engineering Mechanics, and was felt they were too 
long. possible that the Division should have shorter limit, 
due the nature the papers, than the standard ASCE policy. 

The West Coast Committee needs able mail advance notices 
West Coast Conferences the EMD membership. This necessary 
for publicity, well insure ASCE sponsored papers the con- 
ferences. 

The West Coast Committee does not wish review papers. wishes 
only screen papers that have been accepted for publication the 
Journal and that may available for presentation the West Coast 
conference. The following steps are necessary for this screening 
process: 


Select from the Publications Committee list accepted papers, 
number that appear attractive for the West Coast conference. For 
this purpose, necessary that the chairman and secretary the 
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West Coast Committee receive all reports the Publications 
Committee indicating the status all papers. 

Obtain one more copies papers selected, screened 
members the West Coast Committee. may necessary 
contact the program chairman, also, determine arrangements 
can made for presentation this paper the West Coast meet- 
ing. 

Have papers screened members the West Coast Committee and 
select those suitable for inclusion the program. 

Clear with Program Committee and write author inviting presen- 
tation the West Coast conference. 


PUBLICATIONS COMMITTEE MINUTES 


The following are the Minutes for meeting held May 1960, Purdue 


University: 


Report Secretary Wenk: 


Number manuscripts submitted EMD has increased from about 
year 1955-1956 about year. About per cent these have 
been rejected withdrawn. 

the thirty papers submitted during the past six months, twenty-seven 
were about evenly divided among the fields Elasticity, Mathematical 
Methods, Plasticity and Structural Dynamics. Two papers were con- 
cerned with Experimental Analysis and Analogs, one with Fluid Dyna- 
mics. 

About two thirds the papers submitted are reviewed Technical 
Committees within two months, few take more than four months. 

The present backlog unpublished accepted papers such that after 
acceptance paper must wait from nine twelve months before publi- 
cation. This means that the time interval between submission 
manuscript and its publication the Journal usually 15-18 months 
with minimum months. The change from four issues the EMD 
Journal per year six that planned the Society Publications De- 
partment should reduce the backlog and may allow reduction the in- 
terval between submission and publication low six months. 


Reviewing Standards. The Division Secretary has authority question 
and reject reviews manuscripts that not meet Division standards. 
was suggested that the form sent out the Secretary requesting com- 
mittee review manuscripts reworded remove any vague inter- 
mediate category acceptability, e.g. “publish space available” the 
equivalent. Reviewers should take responsibility terms flat yes 
no. Doubtful cases should rejected. 

Short Technical Notes. This category manuscript was originally pro- 
posed means reducing processing time for short papers. view 
the short reviewing time for normal manuscripts, there appears need 
for it. 

Manuals Preparation. The Commentary Plastic Design, the last in- 
stallment which appeared the April Journal, will revised and ready 
for possible issue manual practice September, 1960. 
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The Protective Construction Manual now being edited the committee 
and the manuscript should completed within few weeks. Haltiwanger was 
requested find out its preliminary publication (for criticism Society 
Membership) can single issue the Journal order avoid the 
separate printing with private funds originally planned. 


NEW EMD POLICY TECHNICAL PROGRAMS 


memorandum dated May 31, 1960, White, Division Secretary, 
informed the Program Committees revision the present EMD policy 
technical programs. Until now technical sessions have been limited almost 
entirely the annual meeting the Society and the programs co-sponsored 
the West Coast Committee. Three factors leading revision policy 
are: (1) Enough papers are received fill ten technical sessions per 
year, (2) the division cannot count more than six sessions the annual 
meetings the Society, (3) has been shown that given the right conditions 
engineering mechanics meetings can successful away from the annual 
meeting. 

After discussion these factors, was agreed the May EMD Ex- 
ecutive Committee meeting that the division should the future plan tech- 
nical sessions meetings additional the Annual Meeting. These might 
(a) connection with regional ASCE conventions, (b) jointly sponsored with 
some other society, (c) pure EMD meetings preferably with help from 
local ASCE group. Certain things were considered necessary for success: 

(a) least two three sessions having related themes, (b) reasonable sup- 
ply potential attendance the surrounding areas, (c) enthusiastic planning. 


ITEMS FOR OCTOBER LETTER 


News items and comments, considered for publication Newsletter 
for circulation with the Oct. 1960 EMD Journal, must your editor’s desk 
Wed., Sept. 28, 1960. Please note the change address. 


COMBINED INDEX ASCE PUBLICATIONS 


For complete coverage the Society’s 1959 year print, there now 
Combined Index covering the Division Journals, Transactions, and Civil Engi- 
neering. Also included are reprints the Proceedings Abstracts that are 
published each month Civil Engineering. The price the Combined Index 
(ASCE publication 1960-10) $2.00 with the usual 50% discount members. 
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The coupon herewith will make ordering easy: 


CUT HERE 


American Society Civil Engineers 
West 39th St. 
New York 18, 


Please send me... copy(s) ASCE 1960-10. Enclosed remittance 


Print Name 


Address 


Zone State 
1960-10 


Donald Dean, Newsletter Editor, EMD 
Department Civil Engineering 
University Delaware 

Newark, Delaware 
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